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A METHOD OF DETERMINING THE OPTIMAL PULSE RESPONSE FUNCTION 
FOR ONE CLASS OF DISTURBANCES 


P. S. Matveev 
(Moscow ) 


A method is presented for the determination of the optimum pulse response 
function for a servo system in the case when the system's input consists of given 
harmonic and exponential functions of time plus a stationary random function, In 
solving the problem, use is made of the relationship between the correlation func- 
tion and Green's function. Two examples illustrate the method. 


In a number of works (for example, [1-3]), there are obtained theoretical results related to the determina- 
tion of optimal characteristics of dynamics systems which are subjected to random disturbances. These results 
are useful in those cases when the functions of time are given as polynomials of degree r, the solutions being 
obtained synthetically, Although reference is made in [2,4] to the possibility of using the aforementioned results 
when the disturbances are given in the form of sums of harmonic and exponential functions, no developed method 
has appeared in the literature for determining the pulse response function for these latter cases. The problem 
which the present paper proposes to solve is, using the mathematical relationship between the Greens function 
and the correlation function for a self-adjoint differential system (cf., for example, [5,6]), to give a method for 
determining the optimal pulse response function, given the requirements that the least possible mean-square 
error be attained by a system for which the correlation function of the signal and noise, and the duration of the 
transient response, are given. It is assumed that the signal may be represented in the time domain as the sum 
of given harmonic and exponential components. We shall consider the solution of the problem to be the obtain- 
ing of the reproduction operators: 


H (p) = >) 


i=o 
1. Posing of the Problem 


We shall assume that the system under consideration is acted upon by a controlling stimulus y(t) and a 
disturbing stimulus n(t), Let the controlling stimulus be the sum of two components, y(t) = g(t) + m(t), where 
g(t) is a given function of time in the form of a sum of harmonic and exponential functions and certain of their 
_ derivatives, and m(t) is a stationary random function with a given correlation function R,,(T) or a given spectral 
density S,,(w). 


The disturbing stimulus (noise) is also assumed to be a stationary random function whose correlation func- 
tion R,(T) or spectral density S_(w) is given, For simplicity we shall assume that the functions m(t) and n(t) 
are uncorrelated and have mean values of zero. We shall assume that the quantity to be reproduced is the function 


h(t) =H (p)y a) 


—o 


| 


where 


x(t) = \ H (jo) do. 


—o 


(3) 


We shall further assume that the correlation function R(t) of the sum of the processes m(t) and n(t) has 
2q continuous derivatives and, consequently, that the process g(t) + m(t) + n(t) can be differentiated q times. 
Then, the process x(t) at the system's output may be written (cf. [6]) in the form 


where the Wj(r) are functions of bounded variation which characterize the optimal system. We shall assume 
that they are continuous within the interval of integration (0, T) and have at most first-order discontinuities at 


t = 0 and t = T, such that the Stieltjes integrals in (4) exist. Integration of expression (4) by parts with zero ini- 
tial conditions leads to . 


z(t) + m(t— + (2) dz 


gi 

j=0 

Stee) 
j=o 


where k(r) is the pulse response function of the system for 0 < t < T. 
The system's error of reproduction may be written in the form 


=H (p)y(t)—2(t) = | +. de — 


T 
0 


(6) 
ai 
j=0 
q 
j=0 
We require here that the average value of reproduction error be identically zero, i.e., 
T 
—co 0 (1) 


— — mo — 8 (t—*) = 0. 
j=0 


We assume that in identity (7) g(t) = dor (t). Then, (7) can be rewritten in the form 
v=0 


ayP, (t — t) x (t),dt= 


v==0) 


(t—t) k(t) dt 
v= 0 (8) 


q j 
+ + (t—2) Dy har. 
j=0 v=() j=0 v=0 


We consider further the class of harmonic and exponential functions (sums of sines, cosines and exponen- 
tials) Py (t) which possess properties which permit the representation of Py (t— rT ) in the form 


| 
7 
| 
| 
| 
| 
| 
| 
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= > by (2), (9) 
u=t 


where $j (t) and bi(r) (u =1,2,...,1) are linearly independent functions. By taking (9) into account, we 
may write identity (8) in the following way: 


r 


r T 


v=0 v=0 v=) 


Identity Pr defines (n + 1)r constraints, imposed on the pulse response function by the transforming 
operator H(p), These constraints, as is easily seen, have the form: 


(x) k + 


(t) C5 + 


(11) 
(v=0, i, 2,...,%, 14, 
j=0 


and are constants. 


For simplicity in the argument which follows, we let 


T 
(t —t) + m(t—*) +n (t —*)] (t) de = (t) + (t) + (t) (12) 


and 
+ m(t— 2) + de = 
g(t —T) + (t —T) + oD (t—T). (13) 


Taking (12) and (13) into account, we rewrite (10) and (11) as 


v=0 (14) 
T a 
0 j=0 


Thus, for example, if g(t) = ag cos wet, H(p) = Hy + Hyp and q < 0 then, based on identity (14), one may 
write 


\ COS Wy T COS Wot + Ap SiN Wp t SiN Wot] [H (t) + dt = 


H COS Wot COS Wot do Siz Wo t SiN Wot) + 
HH, [— COS Wot Sin Wot + SiN Wot COS WoT} = (16) 


| 


— 
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(continued) 


T (16) 
=] wor 608 wot + dq SiN Wot SiN Wot] dt = 
0 
T T 
ay 008 wot \ cos wet k de + a sin wet | sin age k (2) de 
0 0 
or 
T T 
(t)dt, Hyw. = sin (t) de. (17) 


In the example given we obtained two constraining conditions, such that n + 1 = 1 and r = 2 and, con- 
sequently, (n + 1)r = 2. If, in this case, we set Hy = 0, the constraining conditions then take the form: 
T T | 
(z)de=0, (e)de = H yup. (18) 


As a second example, we consider the case when g(t) = apsin wot + ay cos 2wyt, H(p) = Hp and q < 0. For 
this example, identity (14) is written as follows: 


( Hg [dy Sin Wot COS Wot — ay COS Wot SiN Wot a COS 2,t cos + 


T 
a, sin sin 3 de a9 [sin cos — (19) 
0 


COS Wot Sin w,t] k (t)dt + a, 
or 0 
T 


T 
H,= (t) dt, = (sin dt, 


[cos 2,t cos + sin 2w,¢ sin 2@,t] k(t) dt 


(20) 
Hy = cos (x) de, 0 = sin 2wyck de. 


In this case, we obtain four constraints, since (yn + 1) = 2,r = 2, and (m + 1)r = 4, 
We now consider a third example, wherein g(t) = e~™', H(p) = Hy, and q < 0. For this example we may 
write 


T 
\% [e—*! e**] 8 (t) dt = e** k(t) dr, 


and, consequently, we get one constraining condition [since here (n + 1) = 1 and r = 1): 


T 
H,= k (x) dt. 


Finally, we consider the example where 


g (t) = a, cosa, H(p)=H, & q=0. 
For this example, one may write 


a 
| 
a 
| 
| 
0 0 
4 


[@q COS Wot COS Wot + SiN Sin wot] b(t) dt = 


T 
[2,008 cos + a, 8in sin wr] (2) + CoB (2) T) = 
0 


T 
ay sin wet sin (2) + C,8(t) + Dd (t—T)) de + 


T 


T 
H,= [Ak (x) + €,8(«) + (« —T)] de, 
0 = (singe (<) + + Do de. 


Here, (n + 1) =1,r = 2 and + 1)r = 2. 


2. Conditions for a Minimum Mean-Square Error 


We write the expression for the error as 


T 


(t) = H (p) y(t) (¢— +) + m(t—) + x 


j=0 j=0 
Then, if identity (14) is satisfied, we have 


T 
(t) = H (p) m(t)—{ {m(t—*) x 
0 


j=0 j=0 


If we square the right and left members of expression (21), we obtain 
e*(t) = \ \ m (t — 0) m(t (8) (x) d0 de — 


oT 
m (t [m(t —2)] (0) x 
x + + [m (t — 6) +n (t—0)]x 
q 


x[m(t—+) +n (0) + C8 (0) + @—7)] x 
j=0 j=0 


q q 
x (<) + | 


a 
5 


and, consequently, 
M [e2(t)] = \ (= — 0) «(0) — 


—o 


+ + + x ton 
00 j=0 


q 
x [A (8) + C; 8 (6) 4 — d= 


where R(t 6) )+R,(r @). 


7 | The problem is now to find a pulse response function k(t) which simultaneously minimizes M[€?(t)] and | 
: satisfies conditions (15). To solve this problem, as is well known, it is necessary to set up the expression 


v=0 


(24) 


and then to add to k(t) its variation Ak(t). 
For the result, we get 


\ Rm — 0) (e) (8) de dd + 

—co —co 

0 j=0 j j=0 


=0 


@ T 
+ + + + (25) 
j=0 0 j=0 j=0 


q qa 
+ A [k (0) + C;8 (6) + (6— dd — 2P de 


O<t<T 


Aly = M [e? (t)] + 2AE, + 


Ey =\ + + («—7)] de x 
r a 
{\Re—9 (0) + (0) + (0—7)| 
0 j=0 j=0 
| 
0 j=0 
T 
| (0) + 8) (6) 4+ SD; 8 (9 — dé, (28) 
0 j=0 j=0 


P(t)= \ Rm (t —*) x(t) dt + >) 


v=9 


| 

q - 

| 

| 5 (26) 


By differentiating (26) with respect to A and then setting A = 0, we obtain the necessary and sufficient 
condition that Al, be a minimum in the form E, = 0, or 


T 
+ (a) + dt — P(t) = 0(<t<7), (29) 
j=0 j=0 


0 
where q = j— k—1, and j and k are defined by the expressions for the spectral densities. 
3. Solution of the Integral Equation 


The solution of integral Eq. (29) will be obtained for the class of stationary processes whose spectral densi- 
ties can be expressed as rational fractions: 


j + “L(@) L*(@) - 


It is possible to show (cf. (5]) that the correlation function of a stationary process which has a spectral 
density of the form of (30) is related to the Greens function of the self-conjugate boundary problem 


L(p)L*(p) y(t) =9, 
(6) lim y™ (¢) (n= 0, 1,...,j—1) 


(30) 


(31) 


by the relationship 
R(t —+) = M (p) M*(p) H(t (32) 


where H(t — r ) is the Greens function of the self-conjugate differential system of (31),i,e,, the solution of the 
equation 


L (p)L*(p) H (t —t) = 8(t (33) 


and the operators M(p) and L(p) equal the M(w) and L(w) of (30) for p = jw, and may be easily determined 
from the expressions for the spectral densities, Taking (32) into account, we rewrite integral Eq. (29) in the form 


T 


q q 
M (p) (p) H(t (2) + 89 (x) + («—7)] de = 
j=0 j=0 


= + Bm (¢ — 2) x(x) de 


M (p) M"(p) + (2) + 7] art = 
j=o 


> (t) + Rm (¢ x(x) de, 


v=0 #=1 


The left member of the linear nonhomogenevus differential equation obtained has constant coefficients 
and is of order 2k. If we denote the integral in the left member of (34) by N(t), we get 


M (p)M*(p)N(t)= + | Rm (35) 
v=0 p=1 
On the one hand, N(t) is determined as the general solution of the differential equation which has the form: 
2k a f 
N(t) = + (p)M* + | 
n=1 


p= 1 


T 


(40) 


and, on the other hand, 


j=0 i-0 


where the Agp are the roots of the equation 


M(QA) M* (\) = 0. (38) 
If we apply the operator L(p)L* (p) to both members of (37), we obtain, on the basis of (33), 


j=o j=0 


0 
=L(p)L'(p)N(t) O<t<7), 


k(t) = L(p) (p) N(t) (0<t< T) 


k(t) = Bx L (0) + 


n=1 


+ |  O<t<7). 
In expression (40), terms of the form 


L(p) L*(p) M>(p) = 


v=0 
v=0 
v=0 


— 
M(p) M*(p) | M (je) ° 


L (p) L*(p) jot — LL (je) jot 
M(p)M*(p) 


then, with account takea of the discontinuities at the points t = 0 and t = T, we can write the final expression 
for k(t) in the closed form: 


k(t) = +3 (t) + 


n=1 v=0 p=1 
m(t —*) x(t) dt 
+ L(p)L*(p)M R (t + 


q 
+ (t) + —T) Gri <7). 
j=o j=o 


i.e., 
} (39) 
| In expanded form, 
| 
| and 
| 
| 
| 
| 
8 


The arbitrary constants B, and Ane as well as C; and Dj, are determined from integral Eq. (29) and from 
conditions (15). 


Example 1, Let g(t) = cos wet, H(p) = Hy, m(t) = 6, and Sp(w) = 2a0/(a* + u*). 


In this case, 2k = 0, 2j = 2,q = 0, M(p)M®* (p) = 2a, L(p)L* (p) = a? — p*, n +1 =1,4r = 2, and P(t—Tr) = 
= COS Wot COS WeT + Sin Wot sin WoT. 


On the basis of expressions (35) and (41) we may write 


2aN (t) = 72 cos Wot + wot, 


N (t)= +> + sin 


(0) = (a — p*) cos + 22 sin aye] = 
= A, cost + Aysinat O<t<7), 


By substituting expression (43) in identity (14), we get 
H, (cos wot Cos + Sin Wot sin 6 (t) dt = | (cos wor cos agt + 
0 


-+sin @,t8in [A, cos @,t + A, sin + C,8 (t) + —T)] de. 


T 
H,= «yt cos «gt + A, Sin gt + C,8(t) + D8 (x — dt = 


= Ag[ + Sin 2aye] sin |" 4+- Cy COS + 
+ Dy cos = A, + {sin day?) Sin* + C, + Dy cos mT 


T 
0 = [sin [4p cos + Ay sin + Cod + Dod —T)] dt = 
0 


= sin* + Ai — 8in Cy SiN Wot + 
+ = sin? + — sin 20,7) + D, sin 


By substituting the value of k(t) from (43) into integral Eq. (29), we obtain 


T . 
[A, cos + A, sin + C,8 () — T)) dt = 
0 


= 712. COS Wot + 72 Sin Wel. 
With the change of variable t— r = 6, we rewrite integral Eq, (46) in the form 


0 t 
[4,008 0, (¢ — 0) + A, sin — 6)] A, c08 — 6) + 


t—T 
4+ Aysin (t —0)}d0 Cye-*! + = 


and 
Whence 
(47) 
9 


(continued) = [410 — Aga + Cy (2? + + (A, cos Wot + A, sin wot)+ 


a? + «? 4 (47) 
[Ay sin og? — Aye cos — 
— Aya sin oT — Ato, 008 D(a? = co + 
From this we obtain the following equations: 
— Ao& + Cy (a* + at 
a? + 
SiN — COS — Aza Sin — COS + D, (a* + (48) 
a? + w? 


The joint solution of Eqs. (44), (45), (47), and (48) determines all the unknowns which enter into the pulse 
response function. Moreover, it follows from Eq. (47) that 


2a.A 2a.A 


It is obvious that, account being taken of the discontinuities at the points t = 0 and t = T, the pulse re- 
sponse function can be presented as 


q q 
k(t) =k + (t) + (50) 
j=0 j=0 


and, consequently, in the example under consideration, 
= Ay cos wot + A, sin wot + Cod (t) + DoS (t —T) (O<t<7). 
By using Eq. (50) we can find the expression for the mean-square error for m(t) = 0. With m(t) = 0, in- 
tegral Eq. (29) goes over into the form 


v=0 p=1 


qa q T 
M = \ (=) + (x) + de] Ra (t — 8) x 
0 j=0 j=0 0 


q a (52) 
x [ (0) + (6) + — 7) ] 
j=0 j=0 
Whence, taking (51) into account, we obtain 
v=0 p=1 
(0<t<7). 
In our example, _— 
M (t)] = = (54) 


Example 2, Let g(t) = cos wt, m(t) = 0, S,(w) = N*, H(p) = 1, and q < 0, It is obvious that 


k(t) = Aycoswst + Aysinwyt (0<t<7). (55) 


> 
} 

| 
10 


(COS Wp COS Wot SiN Wot sin Wot) 8 (t) dt = 


T 
= (COS Wot COS Wot Sin Wot Sin wt) (t) dx 
0 


T 
wp Ap 608 + Ay sin wot] de = + Sin 2 + 
0 


+ sin* = 1 


By solving, simultaneously, the equations in (56), we get 


The transient response may be written in the form 
t 


x (t) = \ cos wy (t (x) dt = cos aot + (24 44) (O<t<7). 


Taking Formulas (55), (57), and (58) into account, we consider two cases, 
First Case, Let T = 1 and w, = 0.25. Then Ay = 3.98, Ay = —23.66, and 
k(t) = 3.98 cos 0,25 t — 23.66 sin wot , (59) 


x(t) = 1.99cos 0.25t + sin 0.25t (7.96 ~ 11,83t). (60) 
Figure 1 gives the results computed from Formulas (59) and (60), and also the value of cos 0,25t, 


Fig. 1 


Second Case, Let T = 1 and w) = 5. Then A; = 2.19, A, = —0.38, and 
k(t) = 2.19 cos 5t— 0.38 sin St, 


x(t) = 1,095cos 5t + (0,22 — 0,19 t)sin5t, 


Then, 
or 
(56) 
and T As a 
Sin Wot [Ap COS Wot + A, SiN Wot] dt = sin*w 7 +- T — sin = 0. 
and 
tft) | ft) 
cos Q25t tt) 
10 
2 
06 
a2 
0+ 0 ik 4 I / / sec 
-06 
and 
(62) 


Fig. 2 


Figure 2 gives the results computed from Formulas (61) and (62), and also the value of cos 5t, 


SUMMARY 


1. This paper presented a generalization of the method for obtaining the optimal pulse response function, 
as presented in [1], to the case when it is required to obtain the minimum mean-square error, given the correla- 
tion functions of the signal and the noise, the time functions in the form of sums of sines, cosines, and exponen- 
tials, and the duration of the transient response. 


In [1-4], the question of choosing the time of observation T remained open, Here, as in [3], we have 
taken the time of observation to be the duration of the system's transient response to the signal g(t). 


2, Application of this method may, in many cases, provide an essential simplification in the determina- 
tion of the pulse response function, as compared to the method of work [1], since the approximations of sines 
and cosines may require polynomials of high degree. This fact could also lead to an increase in the mean- 
square error, 


There is practical significance to such a problem in the case when it is necessary to filter signals which 
are represented on a finite interval by means of Fourier series with finite numbers of terms, 


3. The examples herein provided show that there is no particular difficulty in applying the method pre- 
sented here for the determination of pulse response functions. 
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APPROXIMATE DETERMINATION OF THE AUTO-OSCILLATIONS 
IN THE AUTOMATIC CONTROL SYSTEM FOR A SYNCHRONOUS MOTOR 


D. P. Petelin 
(Moscow ) 


The auto-oscillations are investigated in the automatic control system for a 
a motor, where the nonlinear portion is characterized by the functions 
(sino pe, (sind p*e , and (sino p*e. 


The automatic control of a synchronous motor (SM),working on a load,has for its purpose the maintenance 
of the best working regimens, from the standpoints of voltage oscillations and energy loss in the circuit, both of 


the SM itself and of the energy system supplying it. 


As was shown in [1], underlying the best law for controlling the SM excitation is the principle of the con- 
stancy of reactive power output. It was shown in [2] that to implement such a control law as a function of the 
angular shift of the rotating axial field of the stator and rotor — the angle @ — it is necessary to introduce a non- 
linear element — a squaring block — into the control system, 


For the development of automated synchronous electrodrives for such machines as rolling mills, planers, 
broaching machines, grooving machines, power compressors, etc., for which the load varies from zero up to one, 
or one-and-a-half times the nominal moment, it is necessary to take the nonlinearity of the motors into account, 


In our investigation, we shall take the following nonlinearities into account: the angular moment non- 
linearity, of a sinusoidal nature, of the SM; the nonlinearity induced by the reaction of a synchronous machine's 
armature, and the nonlinearity of the squaring block introduced into the system to provide control of the SM 
under a constant output of reactive power. 


An investigation of the periodic regimens of systems with nonlinear elements having sinusoidal character- 
istics was made by E, L. Urman [3] for the investigation of the synchronization conditions for synchronous 
machines, 


The basis for our investigation is the principle of harmonic equilibrium of N, M, Krylov and N. N, Bogo- 
lyubov [4]. 


The block schematic for the system of automatic control of an SM is shown in Fig. 1, We will now set 
up the equations characterizing the processes in the nonlinear control system, The controlled object, the SM, 
is described by the following equatious. 


a) Equation of the SM's rotor motion: 


(1) 


where the first term is the electromagnetic moment of the SM, m is the number of phases, U is the voltage in 
the supply network, Eg is the idling emf, ws is -ne synchronized angular velocity, Xq is the reactance along the 
axis, J is the moment of inertia applied to the SM'‘s shaft, p is the number of pairs of poles of the SM,and Mg is 
the !oad moment on the shaft. 


mUE, J 


> 


—o 
= SM 
= 
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EMA = 
+ 30 00 
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Fig. 1 
b) The equations for the transient response in the excitation winding of the SM, neglecting magnetic 
saturation, and for the electromagnetic transient responses in the stator circuit: 


X4—X, 
Eae == Ea+Ta,pEa, Ea = Fa cos, (2) 


where Ege is the steady-state value of the idling emf, Eq, is the idling emf during the transient response, EQ is 
the emf of the SM during the reactance transients, Tg, is the idling time constant of the SM's excitation winding, 
and X is the transient reactance along the axis of the SM. 


In correspondence with (2), we give the equation for the SM's emf in the following form: 
= (1 + Tap) a—T4,—y — (Ta = Ta, Xa/Xa). (3) 
The steady-state value of the emf, Ege, is related to the voltage at the input of the SM's excitation wind- 
ing, Ujn, by the relationship 
Ede = (4) 


where ky, = Xad/ tex; Xad is the reactance of the reaction armature along the axis of the SM, and re, is the 
ohmic impedance of the excitation winding. 


The equation of the measuring element is written in the form: 
Ug = (5) 
where Ug is the output quantity and kg is the gain, 
The equation for the squaring block, connected in series with a linear amplifier, has the form: 


Ue = kaksql (6) 


where U, is the magnitude at the output of the electronic amplifier, being simultaneously the input quantity to 
the EMA, and k,, is the gain of the squaring block. 


With complete compensation for the armature reaction, the equation for the EMA will be: 
(1 + TowP) (1 + Tsp) Vin = keU,, (7) 


where T,y is the time constant of the control winding, T, is the time constant of the short-circuited network, 
and k, is the voltage gain. 


Solving Eqs. (1), (3) and (4)-(7) simultaneously, we find 
af? sin & + B, (sin 6?) p 6 + 8, (sin 6)? p* 6 + B, (sin 6)? p* 6 + 
J Ts PA + Town) (1 + 7s + 
P... ~ 2 
Pm knkeka kg k ak sq 
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where 


_ mu m(Xq—Xq) 
@,Xq’ ©, X4 kikek aksq 
m(X4—X4) 
Ps ®, Knkekaksq ké (Tow Ms Pal 


We can put Eq. (8) into the following form: 
Toy p) (A+ 7, p) (1+ 
(1+ + Ty p) (1+ 
kaksq 


Yr + Bi Ye + Bo Ys + Bs + 


where 
yy = @sin6, y,=(sin6)?p0, = (sin 6)? p?9, Y4 = (sin 6)? p* 4. 


In seeking the auto-oscillations of the nonlinear system under consideration, we shall look for solutions 
in the form @ =@» + Op =O + asinuwt, where @» is the steady-state value of the angle and a and w are, respec- 
tively, the amplitude and frequency of the auto-oscillation being sought. With this, p@ = awcoswt, p'e= 
= sin wt and = —aw* cos wt. 


Letting wt = U, we have 
= % +asinU, p§=awcosU, p?4 = —aw*sinU, —aw* cosU. (10) 


By using the method of harmonic linearization [4] for the nonlinear functions entering into Eq. (9) (cf. 
the Appendix), we can separate out from (9) the equations for the periodic and the constant components. 


The equation for the periodic components has the form: 


J P(t + Top) (1+ Ts + Ta P) 
+, Prot Faken = 0. 


The equation for the constant sacieiines has the form 


a (a, Oo) + By (a, @) + Be q, (a, @) + Bs (a, fo, w) +- 
J TowP) (1+ 7, p)(1+T% P) (1 + Tewp) (1 + T's p) (1 + p) 
Pm knteka ksq Inkeka ksq 


Equations (11) and (12) provide the possibility of obtaining the desired parameters of the auto-oscillation. 
However, in our case, the joint analytical solution of Eqs. (11) and (12) does not appear possible. Therefore , to 
determine the desired parameters of the auto-oscillation we employ, for some previously given value of @», a 
graphicoanalytic approximate method [5], using the Mikhailov stability criterion. 


+ 


The characteristic equation corresponding to the equation for the periodic components has the form: 


| 


8, 9a(@, 9%, @) p? +p, p+ E 9, (a, 94) +B, 9; (a, 


— @ @ 
(a, 05, @) 
+ Bs ———— | p + [xq (a, + 81.9" (a, %, @) + Bs (a, + (13) 
J p) +7) p) 
= 


We substitute p = jw in Eq. (13) and then, for each of several numerical values of w, we construct the 
complex plot as a varies from zero to infinity. By using the method of successive approximations, we find those 
values of the frequency w for which the origin of coordinates lies between two neighboring curves. Then, by 
constructing lines of identical amplitude, we isolate the region within which the origin of coordinates lies. Then, 
by the use of interpolation, we can determine the desired auto-oscillation parameters for the control system, It 


is obvious that the smaller the region inside of which the origin of coordinates lies, the more accurately are the 
auto-oscillation parameters determined. 


The parameters thus obtained for the assumed auto-oscillation are investigated for stability by the use of 
the analytic relationship 


Oy Ox 


characterizing stability according to the Mikhailov curves [5,6]. Here, x and y are, respectively, the real and 


imaginary parts of the Mikhailov curve, and ap and wp are, respectively, the amplitude and frequency of the 
auto-oscillation sought. 


Since the oscillatory process in our control system is characterized by a fifth-degree equation (13), we 
must require,as an additional stability condition on the periodic solution found by us,that the Mikhailov curve, 
constructed in accordance with (13), successively traverse five quadrants in the counterclockwise sense. 


If the automatic control system of the synchronous motor has a sufficiently large gain then, as follows 
from (8), the values of the coefficients 8,, 83, and Bs will be small, and the influence of the nonlinearity on 
the reaction of the synchronous machine's armature will be insignificant. With this, the characteristic equation, 
(13), can be expressed in the following form: 


, q, %) (1+ Ty p) (1 + 
a, 9.) + p+a = 0, (14) 
where 


In this case, for given values of 0, it is simpler to use the frequency method of L.S, Goldfarb [7] to deter- 
mine the parameters of the auto-oscillation and to investigate their stability. The Gol*dfarb method makes it 
comparatively simple to gauge the effect of the flywheel mass on the synchronous motor's shaft and the control 
system's gain on the auto-oscillation. 


In accordance with the method presented above, a trial computation was carried out to determine the 
auto-oscillation in the automatic control system for a synchronous motor characterized t) .~e following data. 


The power of the synchronous motor is P = 3.6 kilowatts, kr = 3.76, Tq, = 0.53 seconds, and Tg = 0.163 
seconds, 


The power of the electromachine amplifier is P = 4.5 kilowatts, ke = 10, Tew = 0.04 seconds, and Ts = 
0.1 second, 


The total gain is kaksgkg =1. The flywheel moment of the synchronous electrodrive attached to the 
motor shaft was taken to be, in one case, 0.95 kilograms” and, in the other case, 300 kilograms*, The steady- 
state values of the angle 6, were taken equal to 15°, 25°, and 35°, 
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Fig. 2 


As sli the computations, there was no auto-oscillation in the automatic control system for GD* = 
= 0.95 kilograms’. Fig. 2 shows no intersections between the amplitude characteristics of the nonlinear part 
(curves 1, 2 and 3) and the reciprocal amplitude-phase characteristic of the linear part (curve 4) and, consequent- 
ly, there are no periodic regimens. For GD* = 300 kilograms*, identical values for the auto-oscillation param- 
eters were obtained using the methods of E, P, Popov and L, S. Goldfarb, On Fig. 2, the periodic solutions 
characterized by points My, Mg and Mg correspond to stable auto-oscillations in the control system, but the 
periodic solutions of points Ny, Nz and Ns characterize unstable oscillations. It is clear from Fig. 2 that with 
decreasing angle 6, the amplitude characteristic of the nonlinear part contracts, and approximates to the axis 

of abscissas. In the limiting case when 0, = 0, the amplitude characteristic of the nonlinear part coincides with 
the axis of abscissas from —oo to +o0 as a varies from zero to infinity. In this case, there is always, for Eq. (14), 
an intersection of the amplitude characteristic of the nonlinear part (the axis of abscissas) with the reciprocal 
amplitude-phase characteristic of the linear part. Based on the previous investigation, the following conclusions 
may be drawn. 


1. For small values of the angle 6, the automatic control system for a synchronous motor is subject to 
auto-oscillations, 


2. A significantly large flywheel mass on the shaft of the SM is capable of inducing auto-osciilations in 
the system, so that in developing automation for synchronous drives it is necessary to decrease the flywheel mass 
as much as possible, This step is very useful in decreasing the torsional oscillations of the motor shaft when the 
applied load varies periodically, and in avoiding damage to the shaft. 


3. As the gain of the control system is increased, the auto-oscillations engendered by the nonlinearities - 
considered here die out. 


Appendix 


Using the method of N. M, Krylov and N. N. Bogolyubov, the harmonic linearization of the nonlinear 
functions = sin, = (sine) ys = (sin) and y, = (sine) p*@ is carried out in the following manner: 


thy 
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a (0, 09) + 4° (a, 04), 


(4, 


Ya = (4, Oo» ++ 9” Op, + p0 


Pp’ 


Oy, @) + q2 (a. @) 
P% 


ya = (a, @) + 
ua == (a, + (a, 99, @) + eax 
The coefficients a of the linearized nonlinear functions are constant components and the coefficients of 


the Fourier series when these functions are expanded in Fourier series of the first harmonic. 


By taking into account the fact that the functions sin(asin U) and cos(asin U) may be expressed in the 
’ form of rapidly converging Neyman series [8] 

sin (a sin U) = 2 J; (a) sin U + 2 J3(a)sin 3U + 2J5(a)sin5U +..., 

cos (a sin U) = J, (a) + (a) sin2U + 2/,(a)sin4dU +..., 
where the coefficients 


a” a? 
= 2” nl (: 2x4 (2n + 2) (2a +4) 


are Bessel functions of the n-th order, first type, and by also neglecting the higher harmonics in the Neyman 
series, since the control system is actually a low-frequency filter, we may present the nonlinear functions y,, 
Y2» Ys and y, in the following way: 


vi = sin [J, (a) + 2J2 (a) sin + 02 cosO, 2J; (a) sin U + 20, sin @ [Jy (a) + 


+ 2J3(a) sin 2U] sin U +- 409 cos OoaJ; (a) sin? U + sin 09a? [Jo (a) + 2Ja(a) sin 207) + 
+ 2 cos 0, a* J; (a) sin® U, 


y2 = [sin? 0, a J? (a) cos U + 2sin® 0, a J, (a) Jz (a) (sin U + sin 3U) + 
+ sin? 0, a J? (a) (2 cos U — cos 3U — cos 5U) + sin 20, a Jy (a) Ji (a) sin 2U + 
+ 2 sin 20, J; (a) (a) sin? 2U + cos* 6, a (a) (sin 3U — sin U)jo, 
Ys = (sin? 0, a J? (a) sin U + 2 sin® 0, (a) J2 (a) (cos U — cos 3U) + 
+ sin? (a) (2 sin U 4- sin — sin 5U) + 2 sin 20, a J, (a) J; (a) sin? U + 
+ sin 20, a J; (a) Jz (a) (2 sin'2U — sin 4U) + 4 cos? 0, a (a) sin? 
Ys = — [sin? 0, a (a) cosU + 2 sin? 6, a Jy (a) J2(a) (sin U + sin 3U) + 
+- sin? 0, a J3 (a) (2 cos U -- cos 3U —- cos 5U) + sin 2 0, a Jo (a) Ji (a) sin 2U + 
+ 2sin 20, 4 J; (a) J; (2) sin? 2U + cos? 6,4 (a) (sin 3U -- sin U)' 


With this, the coefficients of the linearized nonlinear functions, defined as 
or on 


on 
a, (a, 9, @) = | ydu, q? i (a, 0,, @) = y sin udu, qi (a, = y cos udu, 


q,, (a, 9) = 02 sin Oy J, (a) + 20, cos 0, a Jy (a) + 5 sin 0, a? J, (a), 


q’ (a, 0.) == 202 cos Ji (a) / a 4- 20, sin Oy (a) +- 0080) J; (a), 
(a, 99) = 209 sin 0» J, (a), 


| 
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4, (@, = sin 2 0p @ Ji (a) Ja (a) @, 

q" (2, Oo, @) = 2 sin® Jy (a) Ja (a) « — cos* J? (a) w, 

9; (@, 0, = sin® J? (a) + 2 sin® J? (a) w, 
Oo, @) = — sin 209 Jo (a) J (a) aw, 

qa (a, 80, = — [sin® Oo (a) + 2 sin® O J3 (a) w® + 2 sin 209 Jo (a) Ji (a) + 

+ 3 cos? J? (a) 

(a, == — 2 sin? 0) Jo (a) J2 (a) 
(a, 99, @) = — sin 269 a J; (2) Je (a) 

(a, Oo, = cos? Oo J? (a) — 2 0, Jo (a) Js (a) w?, 

(a, Oo, @) = — sin® J? (a) — 2 sin® J3 (a) 


The coefficients a thus found are constant for constant values of a and w. In making the transition to 
oscillatory processes with variable a and w, the values of the coefficients also change. 
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TRANSIENT RESPONSE IN A DC CIRCUIT 
CONSISTING OF AN OHMIC RESISTANCE AND A THERMISTOR 


A. G. Shashkov 
(Minsk) 


For a circuit consisting of a thermistor and an ordinary resistor, the block 
schematic is given and the transfer function and stability conditions are found. A 
method is presented for constructing the transient response in this circuit. The in- 
vestigation is based on a common idealization, the averaging of the thermistor 
temperature over its conductance. 


An electrical circuit consisting of an ohmic resistance and a thermistor reacts not only to variations of the 
voltage on its terminals and to changes in the magnitude of the ohmic resistance included in the network with it, 
but also to variations of the thermal relationship of the thermistor with the ambient medium, 


Many works have recently appeared which are devoted to the study of .ue physical properties of thermistors, 
and to their applications [1-5]. In them, the properties of thermistors are predominantly described and studied 
in stationary regimens, The few works in which the nonstationary flow of current through thermistors has been 
investigated considered primarily foreign types of thermistors and, moreover, did not treat the questions fully. 
The present paper is devoted to the study of a dc Ry, R* circuit and to the construction of its transient response, 
The analysis and computation of the possible states of such an elementary circuit are the first necessary steps in 
the study of more coraplicated thermistor circuits, The investigation is based on a common idealization -the 
averaging of the thermistor temperature over its conductance. The temperature of the thermistor is considered 
to be a temperature T, uniform for the entire thermistor, the magnitude of which, when substituted in the formula 
Ry = Rg exp B/T, will give the value of the resistance Ry which is actually measured. 


Circuit Equations and Transient Response 


Any cause whichengenders achange in the conditions of heat transfer between the thermistor and its en- 
vironment will lead to a redistribution of voltage and power dissipation between 
the elements R and Ry (Fig. 1). This process continues until the following 
equations are satisfied: 


Po =Po—Pr, 
Px = 1,0 =-aF(T,— 69), 
Fig. 1 


Here, I, and Ut» are, respectively, the current flowing through the thermistor and the voltage drop across 
it in the steady state, a is the coefficient of heat transmission of the thermistor with surface area F through the sur- 
rounding medium. If there is a stable equilibrium state then, at the termination of the transient response, the 


* We will thus denote a circuit consisting of a resistor and thermistor in series. 
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thermistor is heated to temperature Ty. In the process of temperature variation, the power dissipated by the 
thermistor is partly dissipated to the surrounding medium (power P,) and partly acts to vary the thermal content 
of the thermistor itself (power P,). Thus, 


P = Py + Pp, 
he. 
P=IR,, 
P, =aF(T — 6), 
P,= ve, 
R (1) 
R, = Ree’ , 
Ve 


Here, c is the thermal capacity and V the volume of the thermistor, and t is the time, 
From the system (1) of equations describing the transient response in the circuit, we find that 


aT 
c T 


It is clear from Eq. (2) that a circuit containing a thermistor is sensitive to changes in any quantity which 
enters into the right member of the equation. With changes of circuit parameters or in the thermic relationship 
of the thermistor with its environment, the circuit may become unstable, leading to a rapid increase or decrease 
in the current flowing through it. In order to derive stability conditions for the circuit, we consider system (1) in 
the form of linear approximations to deviations from the equilibrium state, introducing the symbol A to denote 


a deviation: 
AP = AP, + AP», 
AP == 2U,,AI +- PAR,, 
AP, = a)F(AT — Af) + F(T’, — %)Ac , 


dA (3) 
AP, = cV 
1 U, 
Al = Tit Tap + AR). 


Here, By = B/ n. 


The quantities A@, Aa, AR and AU, are mutually independent. Therefore, we may take any of them as 
an input quantity depending functionally on the circuit values, As the output quantity, we may take any one 
(or several) of the remaining “unknowns,” since they are mutually dependent. 


We now consider two basic cases, the reaction of the Ry, R 
circuit to a change in the temperature of the medium surrounding 


4i |v the thermistor (the input quantity is A@), and the reaction to a 
W (Rah change in the medium or in its state (the input quantity is Aa), 
In both cases, we take as the output quantity the change in 
temperature of the thermistor. The first case occurs when the 
= thermistor is used as a sensing element of the temperature of a 
Fig. 2 stagnant fluid or gaseous medium, With this, the conditions U, = 


= const and R = const usually hold, The coefficient a depends on 
the difference T—@. This difference is defined by the accuracy with which the temperature is to be measured, 
and ordinarily does not exceed a fraction of a degree, We may therefore set « = const (6]. With these condi- 
tions holding, we can use Eqs.(3) to set up the block schematic (Fig. 2) and to obtain the circuit's transfer function: 


| 


AT(p) _ Ko 
AO(p)  Hp+T’ 


where 


(4) 


The block schematic, showing the causal chains, was set up by taking into account the deviations from 
the steady-state value of the temperature of the medium surrounding the thermistor, said variations giving rise 
to variations in the temperature and resistance of the thermistor, This leads to a corresponding power drain by 
the thermistor as a function of the current passed, and to a further variation of its temperature. The impact of 
the feedback loop is determined by the error which can occur due to superheating of the thermistor by the 
measuring current. 


The high accuracy requirements on the measurements determine the position of the working portion on 
the ascending arm of the thermistor's volt-ampere characteristic (here the difference T—@ is small, and so 
Ohm's law is valid), Therefore, the dynamic properties of the circuit being considered can be characterized by 
the time and amplitude-phase characteristics: 


AT (ja) x 
AT = K,A8(1—e 4), AO(ja) +i» 


where w is the oscillatory frequency of the temperature 6, 


It is known that with temperature variations of the surrounding medium the volt-ampere characteristic 
is deformed and, under definite conditions, the R7,R circuit becomes unstable, leading to a rapid increase or 
decrease of current. It is clear from (4) that the condition of static stability of the circuit has the form: 


U 
ar he (Ur —Ur) Io <1. (5) 


The points of the volt-ampere characteristic for which Ut) < Up 
or, what amounts to the same thing, Rt») < R, always correspond to stable 
circuits. 


In the first quadrant of Fig. 3 are shown the volt-ampere character- 
istic of the thermistor for two different ambient temperatures, 0; and 03, 
and the lines 1-3 which characterize the ohmic resistance in series with 
the thermistor, The points labeled with two numbers are the limits of 
stability. The solid lines in the second, third, and fourth quadrants define 
the sets of stable points which can be obtained by a continuous variation 
of the ambient temperature 9 from 6, to 6,. The unstable points lie on 
the dashed lines, In the temperature interval 6; < 6 < 42, the thermistor 
may be used as a proportional transducer only on the stable segments 
1-1,2 and 2,3-3, However, it is inefficient to work on the segment 2,3-3 
since the static sensitivity of the thermistor, 8) = B/T} , is lower on this 
segment than on the segment 1-1,2 for the same ambient temperature, 
Fig. 3 and the time constant is greater, Moreover, the superheating of the 
thermistor, T— 0, on the segment 2,3-3 can reach a significant level. 


The presence of unstable equilibrium states (points 1,2 and 2,3) of the Ry,R circuit allows its use in circuits 
for thermal protection of machines and for automatic transformer blow-out [5], in which are used processes 
leading from point 1,2 to point 3, and conversely (Fig. 3). 


* Experimental determinations of the quantity cV/aF and subsequent calculations show that, for example, for a 
type KMT~1 thermistor placed in transformer oil, with w = 0,01 second’, the phase shift between AT and A@ 
is about 16,5°, and the amplitude deviation is about 4%, 
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It is not possible to study these processes in linear ap- 
proximations, since this gives infinitely large thermistor 
temperature variations for arbitrarily small increases of ambi- 
ent temperature above 62 or decreases below 6;. Therefore, 
we do not limit ourselves to linearity in considering the pro- 
cess of moving from point 1,2 to point 3, 


Let the ambient temperature 6 increase and, at a cer- 
tain moment, become equal to @,. Then the imbalance of 
the powers drained and dissipated by the thermistor, AP+ = 
= Py(T)— where Pp = — IR) and = aF(T—@6), induces a variation in the thermistor tempera- 
ture, governed by the equation cV(dT/ dt) = AP. 


Fig. 4 


Since there is no analytic expression for AP; (T), the solution of this equation may be cbtained graphically 
(cf. Appendix 2), For this, it is necessary to know the dependence of APy/cV on the temperature, T = 
B/ In [U, — IR)/IRg)], which is easily constructed from the given static volt-ampere characteristics of the ther- 
mistor.* 


Let us consider the second case. This case occurs when thermistors are used as indicators (or signallers) { 
of the medium and its state, characterized by the speed, pressure, concentration, etc. Such a use of a thermistor 
is based on the dependence of the thermistor's heat-transfer coefficient on the ambient medium and its state. 
As before, Eq. (2) describes the transient response of the circuit when the heat-transfer coefficient varies. In this 
case, U, = const, R = const, @ = const, and the system of equations in (3) may be represented by a block schematic 
+ of the causal chains as in Fig. 4. The corresponding transfer function has the form: 


AP 
Aa (p) Hp+i 
where 


We may see, in analyzing Eq. (6), that the working portion of the volt-ampere characteristic must be 
located at the maximum voltage for which the gain Ky is sufficiently large but the influence of the variations 
in ambient temperature 6 on the thermistor temperature is small, 


P, mw 


Py 2U 


P, 


ts Fig. 5 Fig. 6 


We shall consider separately the reaction of the R7,R circuit to a step-wise variation in the thermistor's 
heat-transfer coefficient. This case occurs in level control or medium inspection. 


*The method for determining the magnitude of cV is given in Appendix 1. 
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Let a thermistor be found in a medium which, up to a certain moment, is characterized by the heat- 
transfer coefficient @,; the current through the thermistor, and the voltage drop across it, are defined by point 1 
(Fig. 5). At time t = 0, let there be an instantaneous change in the properties or state of the ambient medium, 
There will then be a corresponding step-wise jump in the heat-transfer coefficient by the amount Aa =a;— G,, 
where a, is the thermistor's heat-transfer coefficient after the moment when the properties or state of the ambi- 
ent medium changed. A variation of the thermistor's heat-transfer by the amount F(T,—0 )Ac induces a varia- 
tion of its temperature in accordance with the equation: 


cV = AP, (T) sign Ac. 


’ + 1 for Aa > 0, 
for An <0. 

The instantaneous value of the power consumed by the thermistor in the transient response is determined 
by the expression Py = I1(U, — IR). The power dissipated in steady-state heat exchange with the surrounding 
medium, as characterized by the heat-transfer coefficients a; and O,, is given by the formulas, respectively, 
Pry = ) = and Py, = — ) = | depends on the sign of Aa and is defined by 
the relationship: 
for Ax>0, 
| Py— Pyq|for Aan<0. 


Figure 6 gives the dependence of Py, Py, and Py: on the temperature. 


The magnitude of Py is given as a function of the instantaneous thermistor temperature during the tran- 
sient response, T = B/In(RT/ Rg), where RT = (Uc — IR)/I. 


The magnitude of Py, is given as a function of any steady-state temperature of a thermistor immersed in 
a medium with heat-transfer coefficient a: 


B 
="). 
* 


The magnitude of Py, is a function of the steady-state temperature of a thermistor in a medium with 
heat-transfer coefficient a,: 


B U 
(Fre =) 
With these functions known, it is simple to construct the transient response, Of practical interest is the 
fall of current flowing through the thermistor, or the voltage drop across it, as functions of 9 or a, These func- 


_ tions are found by the corresponding computations. The transient response induced in the circuit when there is 


a step-wise jump of the series ohmic resistance, or in the voltage on the terminals, can be investigated and con- 
structed in the analogous manner. 


SUMMARY 


1. The block schematics of the causal chains (Figs. 2 and 4) were drawn, the transfer functions, (4) and 
(6) were obtained, as was the condition for circuit stability, (5). 


2, The processes arising in the R7.R circuit given a relay effect, or a step-wise variation in the heat- 
transfer coefficient, may be considered on the basis of the equation cV(dT/dt)= APy. The static volt-ampere 
characteristic and the numerical value of cV are required for the graphic solution of this equation. 


3. The block schematics of the causal chains and the transfer functions of the Ry,R circuit may be used 
for its analysis and design, 
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Appendix 1 


For the analysis and computation of the transient response in a thermistor circuit, it is necessary to know 
the numerical value of the product of the thermistor's thermal capacity and volume (cV). Given oscillograms 


of the transient response, for example, from point 
‘ 1,2 to point 3 (Fig. 7),this quantity may be deter- 
mined from the equation cV (dT/dt) = APy. In 
fact, the right member of this equation can be com- 
puted from the given static volt-ampere character- 
istic, and the left member can be determined from 
the oscillogram of the process I = I(t), constructed 
in temperature-time coordinates (Fig. 8). This 
curve has a point of inflection which corresponds to 
the maximum value of APy = APy max. To deter- 
mine the magnitude of cV, it is convenient to use 
the linear segment of the curve in the neighborhood 
of the point where dT and dt can be replaced by, 
respectively, AT and At, 


7 
420 
580 
3 
2 Ab at, 
At, be 
IOS 
t,sec. ? 7 0 


Fig. 8 


Consequently, cV = AtAPT max/ AT. 


The magnitude of cV, determined in the foregoing manner for a type KMT~-11 thermistor with a nominal 
resistance of Rey = 1690 kilohm, B = 4600°K, Rg = 0.281 ohm, was equal to 1.4 milliwatt second/ degree. 


Appendix 2 


Below is given the construction of the process of current increase in a circuit containing a type KMT-11 
thermistor immersed in transformer oil at a temper- 


Z TABLE ature of 20.7°C and a linear resistance of R = 10 
—- oo kilohm, The increase of current engendered a dis- 
'y ruption of the circuit's static stability due to the 
kohmj 7,*K | mw 
* increase of the oil temperature above 20.7°C. 

0.2 | 4047.5 |305.0| 41.8 0 To obtain the weneene process from point 1,2 
0.3 694.2 | 313.0] 62.5 9.5 to point 3 of the thermistor's static volt- ampere 
characteristic (Fig. 7), we construct the graph (the 

5 414.5 7324, 5. 

1.0 | 199.5 |342.0]199.5| 44.5 right portion of Fig. 8) of APy /cV as a function of 

2.0 93.5 | 362.0)374.0) 128.0 the instantaneous thermistor temperature , calculated 

from the formula T = where Ry = 

22.8 | 407.01822.0] 348. 

8.0 14.0 | 427.0} 896.0} 336.0 = (U, — IR)/1, and Ug = 180 volts, The data neces- 
10.0 8.7 | 446.0) 870.0} 250.0 f tion of the of /eV=f(T) 
12.0 5.2 |468.0|748.0| 60.0 cor. graph of APr/eV= fi 
12.5 A.5 | 476.0 0 
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are given in the Table, One then draws the line ab, parallel to the axis of ordinates, at some distance which is 
convenient for the construction (in the given case, equal to 100) and n equidistant lines [n = (T,— Ty 2/AT, 
AT = 20°] parallel to the axis of abscissas. The distances between the points of intersection of the lines 0-1, 
0,-2, . . « » On-g-m with the line ab and the corresponding equidistant lines give the time increments Aty, At,, 

- ++,» A4ty. During each of these increments of time, the thermistor temperature will change by the amount AT, 
Measuring off the At to the left of the axis of ordinates, we obtain the curve T = T(t). On Fig. 8, the dotted 
line gives the oscillogram of the process T = T(t). Knowing the graph of Ry = Ry(T) (from the Table), we may 
construct the curves showing the dependence on time of the variations of thermistor resistance and current I = 

= U,/(Rr(t) +R). This latter curve is given by the dotted line on the left side of Fig. 7. 
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TRANSIENT RESPONSES IN THE MAGNETIC CIRCUITS 
OF ELECTROMAGNETIC CLUTCHES 


G.M. Flidlider 


(Moscow) 


A method is proposed for the engineering calculations of transient responses 
in the magnetic circuits of electromagnetic clutches consisting of several solid 
members, account being taken of the influence of eddy currents. 


The equations for the transient responses are given, as are the formulas to 
be used in the computations, and an example of the use of these formulas, together 
with a comparison of computed and empirical data, is provided. A more precise expo- 
sition of the influence of eddy currents on rotor motion is included. 


1. Posing of the Problem 


Electromagnetic clutches (EMC) are widely used in machine construction and in automatic devices. In 
the majority of cases in which EMC are used — in high-speed executive mechanisms, electromagnetic brakes, 
relay systems of automatic control, etc, — the essential factor is the behavior of the clutch during transient 
responses determined by its own electromagnetic transient responses when its excitation circuit is switched on 
and off, 


The majority of clutch designs have solid magnetic circuits.* In such clutches the eddy currents have 
decisive importance in determining the course of the transient response when the circuit is switched off and, to 
a lesser extent, when it is switched on. The influence of the eddy currents is particularly significant during 
switching off when, as is generally the case, the external circuit is open and the course of the transient response 
is completely determined by the parameters of the magnetic circuit. However, the existing methods for calcu- 
lating electromagnetic clutches [1-4] do not take into account the effect of eddy currents and, in practice, do 
not allow computation of the transient response when the clutch is released and, moreover, lead to errors in com- 
puting the effect of forced switching-in. The absence of a theoretical solution and of a method for engineering 
calculations, taking into account the inertia of a clutch considered as an electromagnetic system with solid mag- 
netic circuits, renders difficult both the choice of clutch parameters in the design phase and the analysis of con- 
trol systems which contain EMCs, 


The present paper presents the fundamental results of a theoretical and experimental investigation into 
the transient responses in electromagnetic clutches with the effect of eddy currents taken into account. The 
problems to be considered are the finding of general equations for the electromagnetic transient responses and 
the development of methods for their engineering computation. 


The investigation was carried out so as to be applicable to single-disk electromagnetic friction clutches, 
We have also considered the special features involved in computing clutches with ro. normally switched in, 
as well as powdered- metal and magneto-emulsion clutches. The method of investigat.on and our results may be 
used for solving the analogous problems for other types of electromagnetic clutches, as weil as for various electro- 
Magnetic systems with small armature travels (for example, relays). 


*We have excluded certain types of clutches where the executive organ employs powdered metal or emulsions. 
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2. The Method 


Of prime importance in the problem under consideration is the analytic expression for the effect of eddy 
currents on the course of the transient response, To obtain such expressions for magnetic circuits consisting of 
several solid members, one might use the classical methods of field ware of equivalent short-circuited wind- 
ings [5], or reluctance operators [6]. 


We will make use here of the latter method, thus allowing ourselves to employ the well-known advantages 
of the operator calculus. 


This method was previously employed in three cases of which we know for the determination of the inertia 
of solid cores [6-8]. 


In the works cited there was considered the relatively simple question of the effect of one solid member 
of the magnetic circuit on the transient response of the system, Therefore, their manner of employing the 
general method [6] and, in particular, their method for engineering calculations, was applied by these authors 
only to particular cases and cannot be used by us in view of the complex configurations of the magnetic circuits 
of EMCs and the increased requiremnts as to simplicity and accuracy. In our case it is necessary to find a more 
general and, at the same time, simpler method of carrying out the calculations. 


We note, in accordance with the method chosen, the basic statements which determine the sequence of 
the steps in the solution. 


1) The magnetic circuit is rationally decomposed into a series of members, for each of which there is 
defined a magnetic admittance operator Gj(p) and a reluctance operator 


Rj (p) = (1) 
2) The total reluctance operator is computed 
i=k 
R(p) = 2 Ri(p) (2) 
=1 


as well as the resulting induction operator 


W2.10-8 
L(p) = (3) 


3) Using the known L(p) and R(p), we find the ampere-turns operator Wi(p) from the equation 


= U (p) 
+ PL (p) (4) 


and the operator for magnetic flux on switching in, (p): 


= (5) 


We now make the transition to the original functions, i.e., to functions of time. 
4) The analogous problem for switching-off of the clutch is solved. 


5) An analysis of the fundamental equations of the processes is made, and a whole set of questions related 
to the engineering computations of the transient responses is considered. 


Assumptions. Naturally, the analytic solution of the problem posed cannot be carried out with all the 
actual physical conditions taken into account. We therefore make the following assumptions: 1) the magnetic 
Circuit is not saturated, i.e., 1 = const; 2) flux leakage is small; 3) the induction vector B for each portion of 


* While it is theoretically possible to solve a system of field equations in partial derivatives, it is tremendously 


difficult in practice, 
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the magnetic circuit is parallel to the bounding surface of 
the portion; and 4) in the gaps between the solid portions, 
i.e., in the zones in which the vector B changes its orienta- 
tion, the magnetic circuit is ideal, i.e., p 


Veo 


The first two assumptions are usually made, and are 
YLOEII_. sufficiently well satisfied for single-disk clutches. In the 
Yt; | case of parallel bounding surfaces, the third assumption co- 
5 incides with the second. The fourth assumption is valid if 

ts Ny the reluctance of the orientation zones is small in compari- 
son with the reluctance of the magnetic circuit where, of 
-+— — + + course, the air gaps are neglected, and this condition is al- 
ways satisfied. 


itz It should also be mentioned that, subsequently, after 
= ROY the basic relationships are obtained, the first two assumptions 


yy); = er and the corresponding factors be taken into 


Turning directly now to the solution of the problem, 
we decompose the clutch's magnetic circuit into the follow- 
ing portions (Fig. 1): 1 is the rotor disk, 2 is an air gap, 
Fig. 1. Sketch of the clutch's magnetic 3 is the outside cylindrical portion, 4 is the faceplate disk of 
circuit. the magnetic circuit, 5 is the inside cylindrical portion, and 
6 is again an air gap. 


/) 


Leaving aside the air gaps, we find four portions making up the clutch's magnetic circuit, portions which 
are pairwise identical in form — two disks with apertures in their centers and two hollow cylinders. 


The next step is to determine the parameters cf the clutch's magnetic circuit in the transient response — 
the reluctance operator R(p) and the induction operator L(p), 


3. Parameters of the Clutch's Magnetic Circuit in the Transient Response 


For the determination of R(p) and L(p), in accordance with (2) and (3), it is necessary to compute the 
reluctance (or admittance) operator of a disk and a hollow cylinder. 


By solving Maxwell's equations for a portion, assuming an excitation of one ampere-turn [6], we obtain 
the transient magnetic admittance of the disk, for a radial flow of flux, in the form 


Bnd 
Gqd(p) = 4 (n= 1, 3, 5,000) 
(6) 


Figure 2 shows the frequency spectrum of the transient magnetic admittance of the disk. The transient 
magnetic admittance of a hollow cylinder for a flux directed parallel to a generator has the form: 


G (p) = 1+ PT m 


Formula (7) was obtained by solving Maxwell's equations and then employing an asymptotic expression for the 
Bessel functions in the computations for and An. 


ra) [V ra — (— 1)" 
+ re) (Vig + (— 1)" 


Bn = 


where n = 1,2,3,4,..., 


* Omission in original Russian 
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Fig. 2, Frequency spectrum of the disk's field Fig. 3. Frequency spectrum of the hollow cylinder's 

for excitation by one ampere-turn. field for excitation by one ampere~-turn. 
working 


tg is the outside radius of the cylinder in centimeters, r, is the inside radius of 
the cylinder in centimeters, and 


Y = = ra i 


Figure 3 shows the frequency spectrum of (p). 
Fig. 4. Dependence of 
Bic/Bag on Tp. 


a 


It follows from Fig. 4 that, for r;/r,;= 4/2 we may, with an error not 
greater than 6%, take 


(8) 
Bac = Brd 


where n=1,3,5,... 
Thus, expressions (6) and (7) are identical. Their original function in the time domain is given on Fig. 5. 


With the identity n = 2k—1,n=1,3,5,...,k=1,2,3,..., and with the use of the relationship (from 
page 50 of [10]) 


we may write expressions (6) and (7) in closed form 


1 thz (10) 
G(p) = — 
= 
where z and Ty = Typ forn=1. 
By setting tanhz /z = 1/ f (p), we obtain the following expression for the disk's reluctance operator 


Rg (p) = Rud fa(p). (11) 
For the hollow cylinder we find, analogously, 
Re (Pp) = Ruc fc (p). (12) 


We now sum the reluctance operators of all k portions of the clutch's magnetic circuit (k is the number 
of portions having different time constants Ty n): 
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ers 


of 


R(p) =2 Rufi(p) (3) 
=1 


gwar? Going now to relative quantities, we write the 
one term of the bes clutch's reluctance operator in the form 


k 
using two terms of the series — R(p) =R(0) D aft (p). (14) 
im] 


> In expression (14) 

Fig. 5. Transient magnetic admittance of the R(0)= DW Rui , (15) 
disk and the hollow cylinder. ‘ 

Rui 
x Rui (16) 


We note that, in magnetic circuits similar to the one under consideration, it is reasonable to assume that 
the phenomenon of flux displacement occurs also in the air gaps [7]. In this case, the reluctance of the air gap 
is taken into account by introducing an equivalent magnetic permeability, 


k 
(17) 
x + Rg 
Thus, the reluctances Ryj and R(0), as well as the time constant Ty, contain the equivalent magnetic 
permeability we. 


Having (14) we may easily write the expression for the clutch's inductance operator: 


L(p) = L (0) (18) 


afr) 
1 
Knowing (14) and (18), we may now turn to the determination of the laws governing the rise and fall of 
the magnetic flux in the EMC, 
4, Establishing the Flux when the Clutch is Switched On and Off 


We now write the equilibrium equations for the clutch's electric and magnetic circuits: 


(19) 
U (p) = i(p)r +Wp®(p) 10°, 
Wi 
(20) 
By solving (19) and (20) simultaneously, we obtain 
1 
® (p) =O (22) 
> (P) + PT o 


where U(p) = U = const, Ip = U/r, o = &(0), and 


ok 
2 aif; (P) 


Expressions (21) and (22) describe the increase of the clutch's excitation current and magnetic flux when 
the clutch is switched in by a voltage U = const with a.. air gap 6. 


From (22) we find the law of flux decrease: 
Poff P) = Steady (23) 
where 
®steady=Po, Ty = 0. 
We substitute (22) in ( 23) obtaining, after some transformations, 


2 () 


We emphasize that the values of @) and Tyn which enter into fj(p) are chosen for the permeability te 
corresponding to the given gap 6. 


Expressions (22) and (24) describe the transient response with respect to flux in a general form. Below we 
carry out their transformation to forms convenient for engineering calculations. 


(24) 


5. Armature Motion 


The transient response when a friction EMC is switched on may be decomposed into three intervals: the 
starting interval, in which the gap is 6 = 6, = const; the interval of motion, in which the gap 6 is variable; and 
the interval of steady-state flux (with a tightened rotor), in which the gap is 6 = 6 min = const. 


The relationships obtained above were derived for 5 = const, i.e., for an immovable rotor, The dynamics 
of rotor motion turn out to have a definite effect on the course of the transient response. This effect appears 
both within the interval of motion, due to the appearance of an additional counter-emf, and outside of it, since 
the course of the transient response within the interval of motion determines the boundary conditions for the 
steady-state interval which follows it. 


A host of methods for solving the equations of motion of the rotor (armature) of a piece of electromag- 
netic apparatus are known [11-14]. Of these, as was correctly stated in [1], the most useful is the solution given 
in [13] which describes the motion of an electromagnet's armature. However, the unfounded assumption was 
made in [1] that the starting current ig, = 0. This is not in correspondence with the physics of the phenomena, 
and leads to a noticeable error in the computations for even small counteracting forces, As shown by analysis, 
two cases must be distinguished in considering rotor motion in a clutch: a) @,,<<@», and b) Os < &p. 


In case a) the solution in [1] is correct if, instead of U, the expression used is* 
Use = U— iger. (25) 
In case b) the equation of rotor motion becomes linear, In fact, 


d*x 
mF = A[D, + ae 


where x is the anchor travel, ds; + ©; = @ is the magnetic flux, sr is the starting flux, @; = @ — ®gy is the 
“excess” ("dynamic") flux, and 


A = 2/(25 + 10°Soje). 


*Here, and in what follows, ist is the starting current computed with eddy currents taken into account. 
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By transforming (26) we obtain 


d*x 
mM ie = AQ, [20 (27) 


, The greatest value of &, = Ad is reached at the termination of rotor motion. The mean value is 
;me * A&/2, It is obvious that for $,, > A® it is possible to neglect the second term, and to write (27) in 
the form 


d*x 
By solving (28) jointly with the equation 


Uy =W 
we obtain, as an approximate solution, the function 
x=n't, (29) 
where n’ is a constant coefficient which depends on the clutch parameters and t is time, 
The increase in flux from the value #.; is described by the equation 


©, = Ft (30) 


Line (30) is tangent to the flux curve at the moment of starting, i.e., at the point t = ty. Thus, for the 
construction (or computation) of this line, it suffices to know the time of motion 


3 
& 8 
tmo -\/ (31) 
n 


Thus, in the first approximation, the effect of eddy currents on the process of rotor motion reduces to 
a decrease of the acting voltage Usr at the cost of an absolute and a relative increase in i,,. This effect leads 
to a decrease in the speed of flux increase and to an increase in the time of rotor motion, 


Having (29), we may write 5 (t), He (t) and i(t) inside the interval of motion, and then seek a solution 
for the motion to a second approximation. In practice, the increase in refinement thus obtained is nonessential. 


6. Engineering Computations of the Transient Response 


The laws of magnetic flux increase (22) and decrease (24) were presented in operator form, as derived 
above. A direct transition from the p-domain to the time domain is possible, but not convenient here, since 
the complexity of the time functions thus obtained proscribes their use for practical calculations. 


The methods, described in the literature, of obtaining approximate solutions in analogous problems (6-8, 
15], come to the following. 


A. The first terms in the series expansions for Gj(p) are used [6,8]. 


B. The actual cross section of the solid member is replaced by an equivalent cross section whose shape 
permits a summation of the series for Gj(p) (7). 


C. Graphicoanalytic solutions are employed [15]. 


The disadvantage of method A is the sizeable methodical error at the beginning of the process (19%). 
Method B is useful only for the special case of a rectangular core cross section (we note that in our case, the 
series for Gj(p) are added directly for the actual cross sections), Method C requires the construction of the 
graph of the denominator of (24) for the operator function, and is therefore extremely arduous, Methods A and 
B are also sufficiently difficult. Thus, none of the methods enumerated satisfies the requirements for an engineer- 
ing solution of the problem under consideration. 
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Here, various variants of an engineering solution will be obtained by the use of approximate values of the 


fi(p) = cothz. (32) 


Variation 1, If the exact solution at the beginning of the process is nonessential, method A should be 
used, i.e., a rapidly converging series, of the form given in (6), is to be used, bearing in mind that 


1 
ctha = = (n= 1, 3, 5,...). (33) 
2 i+ 


Substitution of (33) in (22) and (24) leads to the rational-fraction functions 


® (p) = 


bop* + b, 


® ope (P) = 


— by) p* + — by) (my _, — P 


(34) 


(35) 


mop* + +... my my 


When we make the transition to the time domain for the cases of switching on, (34), and switching off 


(35), we obtain, respectively, 
4 Pit 
D() =O [1 + 2 | (36) 
k k 
where the pj are roots of the denominator of (34),5) Ai = —41,since 6(0)=0, 1+ > Ay = 0 and 
1 i=] 
k 
Doff (t) = (37) 


where the pj here are the roots of the denominator of (35), ' B; = B = 5 = 0.81, since, from (35), & o¢¢(0) = 


= by)/ mg] = B Sp. 


mF tp a 


t. sec 


Fig. 6. Variations of the engineering calculations 
of the transient responses with the clutch switched 
off for 6 min = 9.35 millimeters. 


1 


The expressions for the coefficients aj, bj and mj 
are given in Appendix I, while the calculated curve of 
Soff (t)/ So, constructed in accordance with (37), is 
shown in Fig. 6 (curve a). 


Even though k is ordinarily 2 to 4, the solution 
of the k-th degree equation,which is necessary to deter- 
mine the k roots of the denominator of (37),requires a 
definite expenditure of time. In this sense, there is a 
significant amount of convenience in a variation of 
the approximate solution which is best called a linear 


averaging of the process, 


Variation 2, For a small (principal) root p, of 
the denominator of (35), we neglect the influence of 
the rapidly damped components deriving from the co- 
efficient (1 — 6 ), and we set the terms in (35) which 
contain this coefficient equal to zero. Then, the de- 
nominator of (35) is linearized, and its sole root turns 
out to equal 


(38) 
m 


k—1 


“exact” sums: 
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where 


k 
1 


We then introduce the time constant 


T. a;7';. 
i (39) 


This time constant can be defined as a certain equivalent exponent which describes the process of mag- 
netic flux decrease: 


Dore (t) = 7°. (40) 


The curve for the process, computed in accordance with (40), is shown in Fig. 6 (curve b), Physically, 
the linear averaging is equivalent to the assumption that the solid portions have been replaced, for the basic 
portions (81%) of the flux, by short-circuited windings with no inherent leakage. 


Mathematically, it is a variant of the well-known method of approximate calculation of the small roots 
of a polynomial, 


It follows from a comparison of curves a and b of Fig. 6 that the elementary and simple solution by (40) 
entails almost no increase in error, so that such a solution is very convenient for rapid estimates in the design of 
an EMC, 


For the process of switching on, the linear averaging gives: 
1 1 
D(p) = (1— B) Do pT. + (41) 


The time function corresponding to (41) will be 


t 
@(t)=(1—B)D (i—e _™) + pO (i—e (42) 
In the limiting cases we have: 


a) for Te = 0 (no eddy currents) 
t 
(t) =O, (1 —e 7); (43) 


b) for Ty = 0, Te # 0 (supply current) 
t 


@ (1) =(1—B) D+ PO (i—e 7} (44) 


It is possibie to bring (42) to the form 
(e ™ (42a) 


where the second term describes the effect of the eddy currents upon switching-on. The degree of this effect is 
estimated by the ratio Te/ To. 


Variation 3, Relationships (37) and (40) have an increased error (up to 19%) for small t, For ordinary 
friction clutches this fact has little practical effect on the computed release time. However, for powder clutches, 
by virtue of the discontinuity of the flux functions, and for clutches with normally switched-in rotors, possessing 
significant counteracting forces, it is necessary to increase the accuracy of the computations for small t (for 
large p), An approximate value of the sum 1/ f;(p) is given, for large p, by the equality 


= fer. 

é 

mj 

a 
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where nV pT; 
. Substitution of (45) in (24) gives 
off (P) of (4s) 
where 
k 
te= (47) 
1 
Going over to the corresponding time function, we write 
off (¢) = — (48) 
‘ At switching on, if Ty << Te ("supply current"), the process is described by the equation 
(t) = vit Vi. (49) 
’ It should be borne in mind the formulas (48) and (49) are valuable only for small t ; in practice, the limits 


q of their applicability occur at(1 to 0.5), for falling flux and (0 to 0.5), for rising flux. The computed curve 
constructed in accordance with (48) is shown on Fig. 6c. 


It is obvious that the computation of the process of flux fall will be complete if, on one graph, the curves 
for (48) and (40), or (48) and (37) are drawn. The transition from the functional dependence of (48) to that of 
(40), or of (37), occurs at the point at which the curves intersect, a point which always exist, 


. \ For the switching-in process when Ty << Te (forcing by a large additional impedance Ragq >> rt), 

t Formulas (49) and (44) should be used and, for commensurate Ty and Te, Formula (42). When there is no forc- 
ing, ordinarily T, >> Te, and the eddy currents can be neglected in computing the processes occurring at switch- 
ing-in. 


7. Effects of the Air Gaps and of Magnetic Circuit Saturation 


Under our assumptions, both these factors affect only the magnitude of [te and it is therefore convenient 
to consider them simultaneously. Due either to variations in the air gap or to nonlinearities in the magnetiza- 
tion curve, the equivalent permeability will be assumed to take a new value: 


He # 1). (50) 
Then, all the time constants will also change by a factor of 9: 
Tj=nTj, (51) 
Tj = Tp. (52) 
For (35), for example, we then obtain 
Dot (p) = % — bo) + (ma — by) — (53) 


: By introducing the new operator 
q="p, (54) 


we obtain a function which is literally identical to (35), The roots of its denominator will be qj = py, whence 
the new attenuation values will be 


The values of the numerator obviously do not change. Analogous results are obtained by the substitution 
of # & = Ne in Formulas (34), (48), and others. 


| 


We thus arrive at the following ruie: for a given EMC, the form of the law according to which the mag- 
netic flux is established does not depend, for a motionless rotor on the magnitude of the air gap or on the 
magnetic permeability. The effect of these factors reduces to a variation in the scale of the times for the 
transient response curves, a scale change proportional to the relative variation of the resulting magnetic admit- 
tance of the magnetic circuit. 


The rule just given makes it easy to find the magnitudes of the gaps 5 min and 45 max from the conditions 
on the required release of braking times and, for powder clutches, to obtain given dynamic characteristics, 


The conclusion may be drawn from what has been here presented that the effect of the magnitude of » on 
the process is small, When it is necessary to take this effect more accurately into account, a piecewise linear 
magnetization characteristic (in two pieces) should be used, in conjunction with the well-known rules of curve 
fitting. 


8. Taking Leakage into Account 


The phenomenon of leakage corresponds to the appearance of a parallel magnetic admittance G, which 
can be introduced into expressions (14) and (18) and be thus taken into account analytically. 


The admittance Gg is inertialess with respect to the eddy currents and somewhat decreases their effect. 
On the other hand, the additional component of the magnetic flux,corresponding to the leakage admittance,leads 
to a certain decrease of the magnetic potential at the boundaries of the air gap and, consequently, somewhat 
retards the increase of flux. As a result of the interaction of these mutually opposed factors, the effect of leak- 
age on the transient response is small, A more detailed investigation shows that when there are eddy currents, 
the leakage coefficient varies during the transient response and is an additive factor to the transient impedance, 


Due to the small effect of leakage on the transient response, the introduction of a leakage coefficient 
operator into the computational formulas is not worthwhile, since this would significantly complicate the com- 
putations without leading to any essential change in the final results. 


Thus, there should be used in the computations only the ordinary “stationary” leakage coefficient, depend- 
ing on rotor position, and the steady-state value of the magnetic flux, leakage being taken into account,should 
be substituted into the computational formulas, 


The error arising from this, as with the error due to approximate treatment of saturation, enters into the 
total error of the computations. 


9. Experimental Data and Computational Error 


Yee : An experimental verification of the computational results* was 


yj 2 performed on the electromagnetic single-disk clutch (Fig. 7), the mag- 
YY netic circuit of which was implemented by the type of magnetic circuit 
YUE YU in the electromagnetic clutch ENIMS, nh A EM (leading material 
A-100-1, 1955). An oscillogram of the transient response of the flux 
itself (Fig. 6, curve d) was obtained by including a testing winding ,with 
a number of turns W, = 30, in series with an integrating RC circuit lead- 
ing to the dc amplifier at the output of the oscillograph loop. It follows 
from a comparison of the oscillogram with the computed curves that 
the maximum error of computation does not exceed 20% to 25%, 


SUMMARY 


. 1. We obtained general equations for the transient responses in 
the magnetic circuits of electromagnetic clutches, taking into account 

i ie es the inertia of the solid magnetic circuits and the excitation circuits of 
the clutches, 


* The experimental work was carried out in the electrical branch of the Experimental Research Institute for 
Metal-Cutting Machines, 
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2, We developed an engineering method for computing the transient responses, allowing compact com- 
putational formulas to be obtained. 


3. The computations made in accordance with the formulas presented are carried out comparatively 
quickly and give results whose accuracy is satisfactory in practice. 


Appendix I 


The coefficients in expressions (34) and (35), for k = 4, are computed from the formulas: 
for (34) 


b, =(1— (n=0, 1, 2, 3), 
b& = 1, 

ay = boT 

a, =m,_,+5,T, (n= 1, 2, 3, 4), 

as = 1; 


k 

my = (1-8) |] 7;; 

my as — Beer] + [4 — Bors] + 

+ [1 — Bas] 717274 + [1 — Bog] T1727}, 

my = (1 — B) — B + + [1 — B + + 
+ [1— B (ar + + [1 — B (a2 + Ti + 

+ [1—B (@2—a4)] + [1 —B (a: + T27 0), 


4 
ms = —B) + Ba T = 
1 1 


m,= 1. 


For the case when k < 4, the missing time constants Tj and relative impedances @j must be set equal to 


Appendix II. 


‘ Example of the Computations, We shall present the computations for the fall of magnetic flux in the 
f clutch whose magnetic circuit is sketched in Fig. 7, using the following data: Beomp = 12,000 gauss; the mater- 
ial is Armco steel; fay = 2500, and p = 0.15 ohm 6 = 0.35 mm; hs Olen be; = 0.8 cm; beg = 1.7 cm. 


We now determine the reluctance of the pieces (all lengths are taken along the center lines): 
for the rotor and front disk 


for the outer cylinder 


Rog = 0.1832 em” 4, 


for (35) 
} 
zero. 
= In = 0,166 cu-!; 
cl 
bey 
| for the inner cylinder 
| 
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iter- 
cm, 


We now find the relative reluctances: 


2uR, PR oy 


— = 0.733, 


The impedance of the air gap is: 


The equivalent magnetic permeability is: 


The equivalent time constant is 


8 
Te = + = 0.0367 sec 


4 


For = 1.54r;, we recommend 


The parabola parameter is 
1 


(a V Ty + To) = 0.283 sec 2. 


Thus, the process of flux attenuation is described by the curves: 


o (t) (t) —t | 0.0967 
= = 
3.95V t, = 0.81 e 


The construction is given on Fig. 6. 
Received March 31, 1958 
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ON THE DESIGN OF A CONTACT FREQUENCY TRANSFORMER 
FOR ASYNCHRONOUS MOTOR CONTROL 


M. V. Meerov, I. B. Semenov and V. Z. Yarina 


(Moscow ) 


The feasibility and efficacy are demonstrated of using a mechanical contact 
rectifier in the design of a frequency transformer for speed regulation of an ac motor. 
The properties of a current transformer, working in an inverter regimen, are explained, 
and methods of solving the commutation problem are set forth, 


Of all the methods for regulating the speed of an asynchronous motor, the most highly perfected is that of 
regulation by variations of the voltage source's frequency, 


A number of works [1-4] are devoted to the frequency regulation of asynchronous motors. In them are to 
be found investigations of the various aspects of the frequency regulation problem, as well as derivations of the 
most efficient rules for varying the frequency and voltage of the supplying source so as to attain the requisite 
operating characteristics. The great superiority, both for ease of use and for reasons of economy, of frequency 
regulation as compared to other forms of regulation (for example, parametric regulation), has been demonstrated, 


However, despite the clear superiority of this form of regulation, it has not as yet been widely adopted in 
industrial practice, The basic reason for this has been the lack of a frequency transformer which is both reliable 
and economically efficient. 


As a source of variable frequency, a transformer must possess a number of properties, the most basic of 
which are: 1) the ability to control independently the frequency and voltage at its output; 2) the ability to vary 
frequency smoothly within wide limits; 3) a high efficiency rating; and 4) reliability and convenience in use. 


If any of the existing rectifier types is taken as the basis of the frequency transformer, such rectifiers to 
operate in the inverter regimen, then it becomes necessary to assure satisfactory commutation for all frequencies, 
for all transitions from one frequency to another, and for various loadings. It becomes necessary, moreover, to 
assure the virtually independent control of the output voltage and frequency. Moreover, since the desire is to 
use the frequency transformer for ac motor control in various drive designs, it becomes very important to have 
the simple<t possible means for switching the transformer from the inverter regimen to the rectifier regimen 
(for example, during motor braking) and vice versa. 


The present paper is devoted to the exposition of some preliminary results from an investigation into the 
design of frequency transformers. 


As the basis for our transformer design we chose mechanical, or contact, rectifiers (MR). Today, MRs are 
built for currents of 20,000 amps and higher for rectified voltages of 400 to 600 volts. Actual usage of contact 
rectifiers has shown that their contact systems work satisfactorily, without repairs, for eight to ten thousand hours. 
Mechanical rectifiers possess high efficiency ratings — up to 99.5%, The two-way conductivity of mechanical 
rectifier contacts gives, in principle, the capability of switching from a cectifying regimen to an inverting 
regimen and back, reversing both the current and the voltage. 


A number of other properties of MRs, making them efficient for use in the design of a frequency trans~ 
former, wili be considered below. 


Mechanical Frequency Transformers 


A mechanical frequency transformer must be comprised of a rectifier portion, transforming ac voltage at 
industrial frequencies into dc, and an inverter portion, transforming rectified current into ac at the desired fre- 
quency. Obviously, the magnitude of the voltage at the inverter output is determined by the voltage at its in- 
put, i.e., by the output voltage of the rectifier. Therefore, the solution of the problem of independent regula- 
tion of the transformer's output voltage reduces to the solution of the problem of controlling the rectifier's out- 
put voltage (whether the rectifier be of semiconductors, mechanical, or something else). 


The inverter portion of the design is the most complicated and, moreover, is decisive for the transformer 
design, Therefore, as the first step in our work we considered the possibilities of designing the inverter portion 
of a mechanical transformer with the condition that at its input would be applied a de voltage to be regulated 
and, in our experimental designs, we used a dc machine as the source of our dc current. 


As the basis for the design of the inverter portion we chose a contact apparatus, to be used in the mechanic- 
al rectifier, the shaft of which controlled the rotation of the low-power de motor to be regulated. The frequency 
of the inverted voltage was thus determined by the speed of rotation of the contact apparatus. 


Certain analogies between mechanical contacts in series and saturated chokes and rectifiers allowed the 
use, as the basis of the electrical portion of the mechanical inverter, of well-known circuits of rectifier inverters, 


Below, we shall demonstrate the possibility of designing mechanical inverters under conditions where the 
system includes forced commutation and interrelated control of a number of parameters. 


As was pointed out above, in the development of any inverter, 
including mechanical ones, it is necessary to satisfy the following re- 
quirements: a) the commutation of the transformer contacts in all 
work regimens must be sparkless; and, b) the form of the transformer 
output voltage curve must be practically sinusoidal. 


In 


Commutation 


Obtaining sparkless commutation of a mechanical inverter, i.e., 
mA ae obtaining sparkless breaking of periodically operating contacts, is pos- 
Contact A! ' sible if the condition is met that, at the moment of contact switching, 
3 tact g the current flowing through the contacts is close to, or equal to, zero, 
In other words, to obtain sparkless commutation, the design must in- 
Fig. 1 corporate some supplementary device which periodically decreases, 
virtually to zero,the current in the contacts to be broken just prior to 
their switching. 


The use of well-known rectifier circuits as the basis of the mechanical inverter permits a decrease, during 
arcover time, of the current in the contacts to be switched, down to zero and even over into the negative domain. 


In theory, if the commutating device were so adjusted, and contact arcover so devised, that the moment 
of contact switching coincided precisely with the moment of zero current,then the commutation would be sparkless, 


In other words, for sparkless commutation in the ideal case, it is necessary that the condition y = \ be 
met, where y is the commutation angle, corresponding to the time during which the current through the contacts 
decreases to zero from its precommutation value, and \ is the angle of contact overlap (Fig. 1). 


Under actual conditions, there will always be a certain “play” of the contacts in any contact apparatus 
(due to slack in the contact apparatus, thermal changes, nonuniformity in contact wear, etc,), and a certain 
“play” of the commutation angle y , arising from the inaccuracies in the functioning of the contact device. 
Therefore, to obtain sparkless commutation in actual devices, it is necessary to construct a certain spike of low 
current, the length of this spike being such that contact switching would occur during the duration of this spike 
even though there might be “play” of the contacts and oscillations of y. Such a low current spike can be obtained 
by the use of a commutating choke with relay characteristics, connected in series with the contacts. With this, 
the commutating device would have to be so adjusted that, after the current had been decreased, a certain mag- 
netization-switching voltage would be applied to the choke. Then, starting with some small value of current, 
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the choke would begin to switch its magnetic polarity, as a result of which a “spike” would be formed whose 
length (duration) would be determined by the choke's parameters and by the magnetization-switching voltage. 


In this case, reliable commutation would be attained solely 
by so adjusting the contact overlap that contact switching would 
always occur during the low current spike, 


In other words, in contradistinction to rectifier transformers, 
where the duration of the simultaneous firing of the two adjacent 
rectifiers is determined by the magnitude of the angle y and is 
controlled in a natural way, in a mechanical contact apparatus it 
is necessary to maintain artificially a definite relationship between 
the average value of the angle yay and the average magnitude of 
the overlap ay. 


This relationship can be written in the form 


v 
Aav =Yav* 
where v is the duration of the low-current spike. 


Just as in rectifier transformers, the commutation angle y is 
determined by the magnitude of the current I,, to be commutated, 
i.e., by the form of the load, by the design parameters of the de- 
vice, and by the frequency of the output voltage, 


In principle, there are three possible ways of fulfilling the 
condition for reliable commutation, (1). 


1. So to choose the duration v of the spike that, for all varia~ 
tions of the commutation angle y and for fixed ) , the moment of 
contact breaking does not lie outside of its limits (Fig. 2a). It is 
on this principle that the majority of contemporary meciianical 
rectifiers carry out reliable commutation. The disadvantages of 
this solution are the necessity for choosing a spike of long duration, 
which leads to a significant increase in the dimensions of the ex- 
pensive commutating choke, and the limitation of the maximum 
current I; by the choke parameters, 


2. So to regulate the parameters of the commutating device 
that with variations of the load or of the output frequency, the 
average value of the commutation angle yay remains fixed, In 
this case, the spike must be so chosen that its duration somewhat 
exceeds the angles corresponding to the “play” of the contacts and 
the “play” of the commutation angle (Fig. 2b). 


Since the duration v of the spike depends, not only on the 
choke's parameters, but also on the voltage which switches its magnetic polarity, variations in the parameters 
of the commutating device must not lead to variations in this voltage. In this case, the maximum value of Im 
is determined by the heating conditions. . 


3. So to control the magnitude of the contact overlap that, when y varies (Fig. 2c), the following condi- 
tion will hold: 


Fig. 2 


v 
Aav =Yav* (2) 


With this, the problem of the commutating device consists in maintaining a magnetization switching volt- 
age on the choke which will keep the necessary duration v of the spike, 
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Thus, to guarantee sparkless commutation it is necessary, in cases 2 and 3, to introduce continuous con- 
trol of certain system parameters. 


As was stated above, the time ty, for magnetization-switching of the choke and, consequently, the com- 
puted width of the spike v , are determined by the chokes parameters and by the magnitude of the magnetiza- 
tion-switching voltage. The magnetization-switching time can be determined from the following relationship (5): 


2B.wQ-10-* 2B.wQ-10-* 
v= — 7, (3) 


where e is the n.agretization switching voltage on the choke, Bs is the saturation induction of the commutating 
choke, w is the numbe- of turns of the choke, Q is the cross section of the chokes core, and w is the angular 


frequency. 


In the well-known inverter designs, the commutating device (capacitance or master generator) possesses 
the property that the commutating voltage is related, directly or indirectly, to the device's parameters, to the 
load, and to the frequency, 


Variations of the frequency, load, or commutating device parameters will inevitably lead to variations 
of the commutating voltage e and, consequently, to variations in the width v of the spike which, in turn, could 
lead to a rupture of the conditions which make for reliable commutation, 


In order to render such a rupture impossible, it would be desirable to have a commutating device which 
would permit independent control of its voltage during the commutation period. Such a device can be obtained 
by introducing into an ordinary inverter design a special pulse emf generator, Such a generator must satisfy the 
following basic conditions: 1) the form of its emf must be that of pulses, alternating in sign, whose duration 
must not be less than A; 2) during the intervals between pulses, the output voltage must equal zero; 3) it must 
be possible to control the pulse amplitude within wide limits on both the positive and negative sides, 


Methods of introducing such pulse generators will be shown during the consideration of concrete commutat- 
ing device designs. 
Reactive Power and Form of the Voltage Curve 


It was stated in [6] that the rectifier frequency transformers in ordinary use transmit only the active cur- 
rent components to the recipient device, 


Therefore, in cases when the recipient device requires the reactive current components, it is necessary to 
set up special compensators at the transformer‘s output, compensators in the form of capacitive batteries or the 
so-called reactive power machines which, in individual cases, may be more powerful than the motor to be 
regulated, With capacitive compensation, control power is expended only on rotation of the contact apparatus, 
and its magnitude does not depend on the reactive power of the recipient device, In the case where a reactive 
power machine (MFG) has been placed at the transformer's output, variations in magnitude or character of the 
load invariably lead to variations of the magnitude and form of the current curve in this machine, which could 
lead to variations in the power consumed in control. In order to guarantee minimum current in the master fre- 
quency generator and, consequently, minimum power consumed in its rotation, it is necessary to regulate the 
emf of the master frequency generator in such fashion that it coincides, in amplitude and phase, with the first 
harmonic of the voltage on the load. In this case, the only current flowing through the master frequency gener- 
ator is that determined by the difference between the voltage on the load and the sinusoidal voltage of the same 
emf, i.e., only the higher harmonics of the current. 


Inverting Design 


Figure 3 shows a standard single-phase inverting design with a master frequency, where A and B are con- 
tacts of the contact apparatus, DCS is the dc current source, MD is the motor drive of the contact apparatus, 
MFG is the master frequency generator, L¢ is the filtering choke, L, is the commutating choke and L, and Ry 
are, respectively, the inductive and active impedances of the load, 


To guarantee sparkless commutation it is necessary so to choose the angle a of contact- making that con- 
tact arcover occurs during the time when the commutating current from the master frequency generator is 


directed toward the working current in the contacts to be 
switched, This condition holds during half a period and, 
consequently, it is theoretically possible to obtain commuta- 
tion when the angle a varies within the zone (—1 /2 + 

+ 9MFG~ A) < < /2+ A) Fig. 4). 


The magnitude of the voltage applied to the choke 
during commutation is determined by the position of con- 
tact arcover with respect to the master sinusoid, i.e., by 
the contact- angle « and the amplitude of its value, But the 
condition of obtaining the minimum current in the MFG 
engenders the necessity of controlling the angle « and the 
amplitude of the master emf, With this, the magnitude of 
the commutating voltage and, consequently, the width of 
Fig. 3 the spike, are varied which, in general, leads to a rupture 
of the conditions for sparkless commutation, 


In the design it is possible to include, in series with the master frequency generator, a pulse generator 
satisfying the conditions cited above, and so tuned that its pulses coincide with contact arcover (Fig. 4). Such 
a design would allow, by the control of the magnitude and sign of the pulses, the necessary voltage of the com~ 
mutating device to be obtained. With this, the magnetization switching voltage on the choke will not depend 
on the phase or amplitude of the master emf. As a result, there will be reliable sparkless commutation of the 
transformer in any regimen of functioning. 
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With capacitive commutation, the voltage on the capac- 
itor at the moment prior to commutation is switched. This 
voltage is determined by the load parameters, the frequency, 
and by the magnitude of the capacitance, To cbtain the re- 
quisite duration of the spike, it is necessary to be able to con- 
trol this voltage independently of the parameters cited, Such 
control can also be effected by connecting a pulse generator in 
series with the capacitor, 


It follows from what has been said that it is necessary, in 
implementing normal operation of a mechanical frequency 
transformer, to regulate simultaneously several parameters. 


Actually, when there is a variation in the regimen of 
work (speed, torque) of the asynchronous motor to be regulated, 
it is necessary to change, in a definite way, the output voltage 
and frequency of the transformer, There occurs with this a 
variation in the magnitude and character of the load on the 
transformer, 


Corresponding to this, there should be changes in the tun- 
ing of the phase and amplitude of the master generator (in the 
case of an inverter with a master frequency) or in the magnitude 
of the capacitance (in the case of inversion with capacitive 
commutation), i.e., the parameters of the commutating device 
should be changed, Simultaneously with this, so that the con- 
ditions of reliable sparkless commutation may hold, there should 
be a corresponding regulation of the magnitude of the commutat- 
ing voltage, i.e., of the amplitude of the pulse generator's emf, 


It becomes clear, from the consideration of the processes 
in a transformer to work on an asynchronous motor, that varia- 
tion of any of the parameters to be regulated leads to variations 
of the other parameters subject to regulation. The end result is 
a complex control system with several interdependent quantities 
to be regulated. 


The present paper describes just the basic elements of the design of the inverter portion, and presents 


some experimental results. 


The contact apparatus (Fig. 5) was prepared at the IAT AN SSSR in accordance with a design of the Lenin 
Polytechnic Institute at Khar'kov. The permissible current across the contacts was 200 amperes. 


In contradistinction to ordinary designs of contact apparatus, applicable to mechanical rectifiers, worm 
gears and handwheels were used here, which enabled the angle of contact making to be varied during operation 
of the apparatus. The axis of the apparatus was maintained horizontal, Changes in the contact time, i.e., the 
establishment of desired contact arcover, was done at the forward portion of the contact apparatus. This device 
permitted the variation of the relative positions of the immovable and movable contacts. 


The three~phase generator of sign-alternating pulses (type GZI-1) has a maximum voltage of 110 volts, 
a maximum pulse current of 100 amperes, and a frequency of from 30 to 50 pulses/second. The pulse duration 
is from 15 to 20 electric degrees, The generator was prepared by the Lenin Polytechnic Institute at Kharkov 
with a type PN-85 machine as the base, in accordance with the technical conditions imposed by the IAT AN SSSR, 


The oscillogram (Fig. 6) shows the curve of the pulse generator's emf. 


The cc nmutating chokes were prepared at the IAT AN SSSR. The material of their cores was KhVP, As 
the de source there was used a type PN-100 de machine, the voltage of which was regulated. 


In the investigation of designs with master frequency generators, the latter device was implemented by a 
synchronous generator with sharply defined fields of 4.5 kilowatt power. 
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The shaft of the MFG and the shaft of the pulse generator were mechanically linked by means of a clutch, 
which allowed the shafts to be matched at the required angle. There was a rigid connection between the shaft 
of the pulse generator and the shaft of the contact apparatus, This connection was so established that contact 
arcover coincided in time with the application of the pulses, 


The whole aggregate was hooked up to the rotation of the machine to be regulated, a type PN-68 de 
machine. 


A schematic of the circuit, and the timing diagram for the contact operations, are shown in Fig. 7. The 
following notation is used on the figure: PG is the pulse generator, L, are the commutating chokes, L¢ is the 
filtering choke, and DCS is the de source, 


The experimental investigation showed the following: 
1. With the proper choice of system parameters, the voltage and current curves at the inverter output are 


. virtually sinusoidal, 


2. With variations of the load of the asynchronous motor to be regulated, it is necessary so to regulate 
the contact-angle @, the amplitude of the master frequency generator’s emf EyyFG, and the emf of the de source 
Epcs that the following conditions hold: 


_2V 
Encg 


where zy is the total load impedance, x; is the inductive load impedance, R; is the active load impedance, 
and yg; is the angle of phase shift between the load current and voltage. With this, one attains minimum cur- 
rent in the master frequency generator, 


3. The power consumed in the rotation of the aggregate is an insignificant portion of the transformer's 
output power, 


4, Since the master frequency generator had significant inductance, the process of commutating current 
increase was protracted, which significantly decreased the effectiveness of the pulse generator and led to the 
impossibility of obtaining sparkless commutation in all regimens of operation of the asynchronous motor being 
regulated with minimum current in the MFG, 


The usefulness of the design could be enhanced by a proper choice of the parameters of the MFG, speci-~ 
fically, by decreasing its inductance, 


The design of a contact frequency transformer with pulse-capacitance commutation differed from the one 
described above in that the MFG is replaced by a capacitance. 


Completely satisfactory commutation was obtained experimentally in controlling an asynchronous short- 
circuited motor of 7 kilowatt power, both free-running and underload, in the range of 6 to 60 cycles, 


With this, when the speed of the regulated asynchronous motor was varied, the pulse amplitude, the mag- 
nitude of contact arcover and the capacitance were regulated, 


With this design it turned out to be possible to carry outstarting and speed regulation of an asynchronous 
motor of 28 kilowatt power when it was running without load. 


Received January 22, 1958 
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MULTIFREQUENCY PULSE DEVICES FOR TELECONTROL 


F. A. Katkov 
(Kiev) 


The basic principles are presented for electromechanical vibrators with pulsed 
excitation which are used to form the frequency signals for telecontrol and tele- 
metering multifrequency devices (TC-TM), Variations of the TC-TM design are 
given for the case when the objects are distributed along the common communica- 
tion channel, 


Multifrequency devices for telecontrol, whose principles of design were presented in [1], require for their 
operation sources of frequency signals which possess significant stability. 


The requirements on the stability of the frequency signal generator depend on the bandwidth for operating 
the relay, account being taken of bandwidth variations. The admissible instability of the frequency generator 
can be determined from the expression 


— 
bf,= OK, 
where A fo, is the bandwidth for relay operation, 6 fre is the shift of relay operation bandwidth when external 
conditions (temperature, voltage, etc.) vary,and K, is the safety factor, lying between 2 and 3. 


For use in multifrequency TC-TM devices with narrow-band electromechanical frequency relays having 
operation bandwidths of two to three cycles, frequency generators are subject to such high stability requirements 
that only electromechanical vibrators should be employed in their construction. Frequency signals are obtained 
most simply from them by virtue of their pulsed excitation. For this, a condenser is discharged across the input 
winding at the proper moment, giving rise to a frequency pulse at the vibrator output with a gradually damped 
amplitude, For telecontrol purposes, this pulse is given preparatory amplification. 


When connecting cables are to be used, where the noise level is low, the frequency signals obtained from 
the vibrator output may be passed directly into the communication channel without preliminary amplification, 
for subsequent use in telemetering. 


For this case, we determine the distance of signal transmission along the connecting cables, As is known 
[2], the noise level in cables with no wireless links equals —7 to —8 nepers in a 5 kilocycle bandwidth. 


Since the band used by the given TC-TM devices does not exceed 200 cycles, it is possible to choose a 
lower noise level for the computations. However, we will start with the average value of —7,5 nepers. 


To guarantee reliable operations of the TC-TM frequency device, it suffices to have an excess of 2 nepers 
of minimum signal level over average noise level. Consequently, the minimum signal level at the amplifier 
input must be —5.5 nepers. 


By taking into account the exponential character of frequency pulse damping at the vibrator output, we 
may write the expression for the power generated by the vibrator during time t after excitation in the form [3] 


Pp = Pye" ddwot 


‘ 
| 
| 
uf 


where Py is the initial power at the vibrator output, dq is the damping coefficient in the vibrator system, and 
Wo is the angular frequency of resonance. 


To determine the minimum value of power at the vibrator output, we should set t = 1.5top, where top is 
the operating time of the receiving frequency relay. 


If we let fp = 160 cycles, dg = 0.01, Py = 1.5 + 10° watts, and t = 0.3 seconds, then Pp = 7.35 - 1075 watts, 


If we take into account the mutual shunting of the vibrators for two-frequency signals, then Pp = Pg/2 = 
= 3.67 + 107° watts, or —5.1 nepers 


For noncoil-loaded telephone lines of 1.2 mm diameter the kilometric damping, for a frequency of 300 
cycles, is B 39 = 0.031 nepers /kilometer, By bearing in mind that the damping is proportional to the square root 
of the frequency, we obtain, for a frequency of 160 cycles, 8 yg = 0.023 nepers /kilometer. Consequently, in 
this case the transmission distance will be 


5.5-5.1 


0.023 


= 17.4 kilometers 

For higher frequencies, the transmission distance will be correspondingly smaller. In order that the com- 
bining frequencies, created when multifrequency signals enter the amplifier input, do not give rise to additional 
noise, the amplifier should be tuned to the working frequency band which is defined by the expression f,, = 
= 2(f; — Affe), where fy and f; are, respectively, the upper and lower frequencies of the working band, and 
Afre =(2 to 5)Afop is the difference frequency of the working frequencies. 


If the noise level in the channel is high and more than two frequencies are simultaneously transmitted, 
then an amplifier should also be used on the transmitting end. These amplifiers are prepared from transistors 
according to standard circuits. 


When vibrators are used as pulsed frequency generators, it is important to determine the length of the ex- 
citing pulse for which the initial amplitude of the free oscillations will be maximal. To simplify the analysis, 
we shall assume that the pulses are rectangular, of width Tr. 


When the leading edge of the pulse is operative, the equation of the free oscillations will have the form 


1 
X(t)= —Ape * 
where A, is the initial amplitude of the oscillation. 


SiN ol, 


At time t =r the free oscillations corresponding to the trailing edge of the pulse arise. The initial ampli- 
tude of these oscillations will have a sign opposite to that of the initial amplitude of the oscillations arising 
when the leading edge of the pulse impinged. Consequently, 


t—t) 
X(t = Aye * a sin — *). 


1 1 
14 
“de 


X(t) + X(t—t+) = Are "sin Wo(t — t)— sin 


From this expression we can determine the greatest amplitude of the oscillations fort >t. For this, and 
bearing in mind that dg << 1, and +r < Ty = 24 /wo,we set 7%d** 4, Then, the expression for the resulting 
oscillations can be put in the form: : ig 


le 


1 
X(t) + X(t — 2) = —24y sin [e Wy 
and the desired maxinium amplitude of the oscillations will equal 


Am == 2Ay sin 


| 
The resulting oscillations, for t > r , will be 
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‘ We shall now determine the pulse length Tm 
7. a for which the resulting oscillations will be maximal. 
— For this, we differentiate the equation for A, with 
TA respect to r , and equate the result to zero: 
4 = Ay (09 608 = 0. 
Remembering that r< Ty, we get Tr = 2. 
a When the vibrator is excited by pulses from a 
va a Ya m2 condenser discharge, the process of setting up the oscil- 
=_——= lations occurs in a much more complicated fashion, 
since the rise and fall of the pulse do not occur in- 
f- — stantaneously, There is importance also in the phen- 
— 
7 Tt omenon of energy exchange and, hence, additional 
0 L weil dissipation, between the condenser and the vibrator. 
/- It has been established experimentally that the optimal 
capacity of the discharge condenser somewhat depends 
am on the vibrator load, in such fashion that a decrease of 
, load impedance increases the optimal value of the dis- 
4 charge condenser'’s capacitance. This phenomenon 
can be explained by the fact that a decrease of load 
Fig. 1. Curves for determining the number of impedance lowers the time constant of the condenser 
frequencies m in a multifrequency signal as a discharge circuit, since the equivalent discharge im- 
function of the required number of commands pedance is lowered. 


N and the number n of working frequencies. It should be mentioned that the dependence of 


the initial amplitude of vibrator oscillations on the , 
capacitance of the discharge condenser, when close to the optimal value, is not significant, which renders pos- 
sible the use of one discharge condenser for several vibrators close to the tuned frequency. 


In constructing a multifrequency TC-TM system, it is very important to determine the number of fre- 
quencies to be transmitted simultaneously, For this, it is convenient to use the curves shown in Fig. 1. It is 
necessary to bear in mind with this that an increase in the number of frequencies entering into the multifrequency 
signals increases the noise stability of the TC-TM device on the one hand and, on the other hand, increases the 
unit number of contacts and increases the power dissipation at the receiving amplifier's output. 


The increase of noise stability of a multifrequency code over a single-frequency one, with independent 
source deviations, can be determined on the basis of the following considerations, 


We shall assume that the probability of incorrect operation of one of n frequency receivers (frequency 
relays) equals P, Then, for a single-frequency code, the total probability of incorrect operation for a set of N 
two-position objects under control, equals 


Py = 2NP, 
and for a multifrequency code, 
Pm = 2N(KpP 


where Kp is a coefficient which expresses the fact that there is an increase in the probability of an incorrect 
operation of the individual receiver relays due to the decrease in power of the signal components of the multi- 
frequency signal with the limitation on their mean level, and m is the number of frequencies simultaneously 
transmitted, 

From a comparison of these two expressions one may draw the conclusion that the use of a tifrequency 


code, for the same number of controlled objects, varies the noise stability by a factor of (Kp. yy”: 
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Fig. 2. Circuit for a pulsed multifrequency telemeteriiig device. 
1VFC,,.nVFC...nVBD are telemetering position vibrators, ICR 
eee MCR are contacts of the relay controlling object position, RR 
is the release relay, Cy. ~. «Cp are discharge condensers, A is the 
amplifier, 1FC...nFC, 1FD ...nFD, are frequency relays, R 

« oRy are additional impedances, 


inv 


Fig. 3. Circuit for a pulsed multifrequency TC-TM device with telemeter- 
ing based on the energy of the control pulses, CP is the control point, 1TF 

«TF are frequency relays receiving telemetering pulses, 1FV .. .nFV 

are vibrators forming the control pulses, A is the amplifier, TR is the trigger, 
relay, 1CK . . «nCK are the control keys, C;, C; and Cp are condensers, CKT 
is the control knob trigger, IKT is the interrogation knob trigger, 1V...nV 
are vibrators of the TC-TM, 1C....mC are the output control relays, and 
1CR, ..mCR are contacts of the relays for controlling the state of the objects. 


Insofar as Kp © m and P << 1, the increase of noise stability obtained by using a multifrequency code 
quite overshadows the decrease in noise stability of the individual receivers due to the decrease in power of the 
components of the multifrequency signals, On balance, the noise stability of the system is increased, as has 
been verified experimentally. 


The average, over one period, of the power of the multifrequency signal at the input to the set of n fre- 
quency relays, connected in parallel, must equal 


Pap = mnPfr, 
where P¢, is the power required by one frequency relay. 


vec] Tac 6° | 
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The input amplifier must be designed for this power . 


A pulsed multifrequency TC-TM device can be built in two basic variants: with vibrator excitation 
om the energy of discharge condensers or from the energy of control pulses. A circuit for a telemetering de- 
ice constructed in accordance with the first variant is shown in Fig. 2. At the controlled point CP there are 
‘t up vibrators in accordance with the number of telemetering positions and individual control relays and sig- 
nal switching, while at the dispatch point DP are the amplifier A and the frequency relays. 


With a change of state of an object, the contacts of control relay CR are switched into the circuit of the 
.ischarge condensers, C;— Cy, and relay TR, and a telemetering frequency pulse is sent through the channel 
causing operation of the corresponding frequency relay at the DP, Relay TR is released after relay RR operates, 
the delayed operation of which is brought about by the presence of resistor R and condenser C. 


A circuit for a pulsed multifrequency TC-TM device built in accordance with the second variant is shown 
in Fig. 3. Vibrators 1V.., nV at the executive point are used both for reception of the control pulses and for 
the sending of the telemetering pulses. 


When the control key is depressed and the control knob trigger is pressed, a frequency pulse is sent along 
the channel, causing operation, at the executive point EP, of the corresponding frequency relay and switching 
in of the controlled object. With this, the control relay CR switches its contacts into the circuit of the vibrator 
winding and the input winding of the vibrator is switched out of the channel while the output winding, with the 
winding output relay, is switched into the channel, When relay TR loses its excitation, the output of amplifier A 
is switched to frequency relays 1TF ...nTF, and its input is switched into the channel, With this, the oscilla- 
tions of the vibrator do not immediately cease and, on the EP side of the channel, there will be sent a frequency 
signal which, after being amplified, will act on the corresponding frequency relay. 


When the interrogation knob trigger IKT is pressed, the duration of the frequency pulse in the channel is 
not sufficient to cause operation of the output relay (since Cy > C,) and, consequently, it becomes possible to 


- interrogate the state of the objects. 


A description of the device and of the characteristics of the frequency relays which were used in construct- 
ing the TC-TM devices considered here are given in [4]. There is therein also given a description of the com- 
plete circuit, and experimental results, for a pulsed multifrequency TC-TM device for centered objects, the 
experimental model of which was built by the Department of Automation and Remote Control of the Kiev Poly- 
technic Institute. The same group developed a pulsed multifrequency TC-TM device for objects distributed 
along a common communication channel. Tue principles underlying the construction of such devices are given 
below. 


Depending on the concrete requirements, one may create many variations of such types of devices. 


For a comparatively small number of executive points and objects at them, one might recommend a 
direct multifrequency combined set of objects at the individual points. 


if the number of executive points is large, or the number of objects at them is large, then the set of execu- 
tive points should first be isolated, after which the object should be isolated. In order, in such a scheme, to in- 
crease the rapidity of action and to simplify the circuitry, the communication channel should carry a multifre- 
quency signal consisting of a combination of the frequency to isolate the executive point and the frequency to 
isolate the object. Only the frequency relays for isolating the executive point (group selection relays) should be 
directly connected to the channel and, upon their operation, the proper block of frequency relays for object 
selection will be connected to the channel, In certain cases, in order to reduce the power dissipated by the 
group selection relays, it makes sense to connect only one of the group relays to the channel and let its opera- 
tion switch the remaining group relays into the channel, 


Telemetering from dispersed objects is carried out in a somewhat more complicated manner than tele- 
control, Use inay be found for the following three variants: simultaneous telemetering, serial telemetering, 
and telemetering by interrogation. 


In the variant with simultaneous telemetering, signals may, as they arise, be accepted simultaneously 
from several points, Since, in this case, signals may accumulate one on another, they must have different 


cr 


Telemetering 
objects 
EP, EP(n-1) 


Fig. 4. Block schematic of a multifrequency device for serial telemetering. 


cP 


CGB Ac 


Fig. 5. Circuit for a multifrequency device with telemetering by interroga- 
tion, A, and A; are amplifiers, CGB is the control generator block, TFRB is 
the telemetering frequency relay block, CB, and CB, are object contact 
blocks, and DR are deblocking relay contacts, 


frequencies. Moreover, the input amplifier must have sufficient power at its output to operate simultaneously 
the frequency relays equal in number to the maximum possible number of frequency signals at the input. “‘he 
adjustment of the receiving frequency relays must take into account that, with the minimum number of fre quency 


signals at the output and their power increased at the amplifier output, there does not occur overlap of the bands 
of relay operation, 


It is clear that the variant just considered can be recommended when there is a comparatively small num- 
ber of executive points and objects at them, However, their number can be increased if the most important 
executive points are given the more “active” codes, This idea can most simply be implemented by so increas- 
ing the power of the signals from a part of the executive points that the signals from these points, when passing 
through the amplifier, will suppress signals from the remaining points, After acceptance of the signals from the 
“active” executive points, the signals from the remaining points may be accepted, the duration of the signals 
having to be correspondingly greater. Several such levels may be used. 


In the variant with serial telemetering, signals may be accepted at each moment only from one of the 
executive points, where the points closest to the control point are given preference in the telemetering. For this, 
whe there is telemetry from an executive point situated close to the control point, all of the tail of the com- 
munication circuit is cut off. After the signal is reproduced at the control point, the latter sends an acknow- 
ledgement (quitting) signal which causes the channel to be reestablished, allowing telemetering signals from the 
executive point next farther removed from the control point to be received. Until reestablishment of the channel 
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to the more remote points, the telemetering signals must be stored. The high speed of the device makes it 
possible to obtain telemetry signals from the most remote executive points even when there is a significant 
number of them, 


A block schematic for such a device is shown in Fig. 4, where each TGB is a telemetry generator block 
and FRB is the receiving frequency relay block, 


With variation of the state of an object, all of the succeeding portion of the channel is cut off by contacts 
of the position relays IN or OU until a quitting signal is received from the CP, Con-:rol of the channel is car- 
ried out by means of a de source on the side of the later EPs, 


Interrogation of object state in the variants so far considered may be carried out most simply by using a 
telecontrol device for switching in and out the corresponding objects. 


In the variant with telemetering by interrogation, data is initially taken from the executive points by 
telecontrol, and is thereafter obtained by interrogation. Transmission may be carried out, as desired, either by 
dispatcher control or cyclically after a definite time interval. The principles underlying the construction of 
such a cevice are shown in Fig. 5. Relay PR is normally open, and amplifier Ay is connected to the channel, 
At interrogation, relay PR operates and a two-frequency signal is transmitted from the CP, selecting the EP, at 
which the corresponding frequency group relays FGR,, then FGR,, operate. Relay FGR, is held and switches 
block frequency relay BFR into the channel, but contacts of relay GR connect the source to generators GSV and 
GSO, After this the CP sends successive combinations of selection frequencies to the interrogated objects, at 
which are operated, and held, relays VR. These relays switch into the channel the corresponding pairs of tele- 
metering generators GIN and GOU, With this, the given EPs send frequency pulses along the channel, thus de- 
fining the state of the interrogated objects. These pulses are received at the CP and act on the previously pre- 
pared signal circuits. What this preparation amounts to is that, in the time intervals between transmissions of 
selecting combinations, relay PR is released and the apparatus at the CP is switched over to signal reception, 


The deblocking (releasing) of group relay FGRg, after interrogation of the following EP, can be imple- 
mented by the transmission of a special combination frequency or, if a remote source is used, by a brief break- 
ing of the supply channel. 


In the circuit considered, interrogation of object states may be carried out episodically, by depressing of 
the corresponding keys by a dispatcher, or cyclically by means of a specially programmed switch, It is also 
possible to construct the device in such a way that, with a spontaneous change of state in one of the objects, a 
warning signal is sent to the CP, this signal either ringing a bell or automatically starting up the programmed 
switch, 


With a comparatively small number of objects, interrogation of states may be carried out by excitation 
of the individual vibrators by pulse control (Fig. 3). 


TC-TM devices, constructed on the basis of the principles considered here, may find use in the remote 
control of many processes in the various branches of the Soviet Union's economy. 


A further problem is the development of a contactless multifrequency TC-TM device. 
Received November 20, 1957, 
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A PRECISION MULTIPLIER 


L. N. Fitsner 


(Moscow) 


A description is given of a multiplier with a static error not exceeding 0.01 
to 0.02% of the full scale of output voltage, and with the ordinary high speed of 
electronic components. The low value of error is made possible by constructing 
the multiplier of a “rough and an “exact™ part. 


INTRODUCTION 


The necessity for multiplying two voltages arises in the design of nonlinear models of dynamic systems 
(1), in certain electrical computing instruments (harmonic analyzers, correlators, etc.) [2], and in a number of 

cases in the automation of productive processes, A device which carries out this operation is called a multiplier 
(abbreviated MU). 


In most widespread use today are pulse~time multipliers [3], diode multipliers [4], and multipliers with 
carborundum resistors [5,6]. The errors of these multiplier types lie within the limits of 0.2 to 0.5% of the out- 
put scale, and seem to be incapable of essential reduction. 


In many cases, the error just mentioned is inadmissible, thus giving rise to the necessity of more accurate 
MU's, The accuracy of an MU can be increased by using compensatory methods in its circuit. This was the 
path followed by Goldberg [7] and by Morril and Baum [8]. In designing pulse-time MU's, these authors used 
precision clutches, based on the compensation principle, which they had developed. 


However, these circuits were comparatively complex, and their adjustment was difficult. There is another 
way to increase the accuracy of an MU, based on the use of the discrete principle, The first to take this path 
was the American investigator Goldberg [9]. 


We shall consider below the circuit for a precision multiplier which was presented by A. A. Fel'dbaum 
and the present author in 1952 [10], this multiplier differing from Goldberg's both in its underlying principle and 
in its quality. It is simpler than the American circuit, containing fewer elements and being easier to adjust. 


The circuit to be discussed differs from the American circuit in that it does not require the use of very 


fast electromagnetic relays whose contacts function in a regimen of high-frequency vibration and, therefore, 
rapidly wear out, 


In a number of instances, the dynamic properties of the MU design,to be described below,turn out to be 
significantly better than those of Goldberg's design, particularly in its transfer function for large deviations. 


A Multiplier with “Rough” and “Exact™ Parts 


The precision MU to be discussed, the block schematic of which is shown in Fig. 1, may be divided into 
two parts. We shall call them the “rough” and the “exact” portions. 


One of the input voltages x, is applied to the input of the separator block, Here voltage x, is separated 
into two portions, the basic portion X, and the supplementary portion Ax;. There is a number of fixed values of 
voltage present in the separator block, That one of them which is closest in value to voltage x, is taken as the 
basic portion X, which is applied to the input of the rough system. 


The second, supplementary, portion 
Ax; =X —X (1) 


ct, Y-CI,L, is applied to the input of the exact portion, In the rough 
anal! | portion, voltage X, is multiplied, with high accuracy, by the 
second factor xq. 


The name “rough” is given to this portion only because 
Fig. 1. 1) Separator block; 2) “rough” one of its factors, X;, differs by Ax, from the factor x. 


portion; 3) “exact” portion; 4) sum- Voltage Ax, is applied to the exact portion, Here it 
ming block. is first amplified with gain k, then multiplied by factor x, 

using an MU of some standard type, for example, one of 
those mentioned above, and then, finally, attenuated by the factor k. 


At the multiplier output there is a summing block for adding the ovtput voltages of the rough and exact 
portions, In this case, the output voltage will be 


= CXyXq = CXyX_ + (2) 


Here, c is the coefficient of proportionality of the multiplier. For (x,)max = (%_)max = 100 volts and 
(ymax = 100 volts, c will equal 0.01. 


Let us now consider the reasons which allow the error in such a combination-type MU to be lowered, 
The output voltage of the MU is 
y = + (3) 
Here 6 denotes the error in multiplication, This error can be represented as the sum of three components: 
5 =5,+6, +5e, (4) 
where 5, is the error of the summing block, 6, is the error of the rough system, and 6 ¢ is the error of the exact 
system. 


If, for the summing block, one uses an operational amplifier with automatic zeroing [11] and with wire 
zesistors in the feedback path, the error 6, can be ignored. 


The error 6, is also negligible,since multiplication in the rough portion is usually implemented, as will 
be shown below, by switching of the impedance at the input of the summation block by means of relay or 
selector contacts, Thus, the error of the combination-type multiplier is close to the error of its exact portion: 


5 (5) 


The maximum value of the reduced voltage Ax, from the output of the separator block is several times 
smaller than the maximum value of the multiplier's input voltage x,. 


Let 
(xy) max/(Ax,)max = 2- (6) 


At the input of the exact system, Ax, is amplified by a factor of k, The gain k is chosen to be close to n 
(ordinarily, k= n). The amplified voltage kAx, is applied to the input of an ordinary electronic MU with an 
error of 59. The voltage at the output of this MU is 


y* = ke Axx, + 4p. (7) 
A voltage divider at the output of the exact system effects a k-fold decrease of y*. 


As a result of these operations, the input of the MU's summing block receives from the exact portion's 
output the voltage 


y = + . (8) 
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we have that 8, 
. (9) 
Consequently, the error of the combination-type MU is lower by a factor of k than the error of an MU of 
the ordinary type. 


Multiplier with a Relay-Type ‘*Rough"™ Portion 


Figure 2 shows the rough portion of a combination-type multiplier consisting of electromagnetic relays 


Reg to Reg, and high-gain dc amplifiers Ay to As. The relay coils are connected in the anode circuits of the 
amplifiers’ output stages. 


In the circuit given, switching of the contacts of relay Rey is determined by the sign of input voltage xj; 
the positions of the contacts of relays Rey~Reg, connected in series with amplifiers Ay-~Ag with several input 
voltages, are determined by the signs of the sums of these voltages. Contact 1Re, is so wired in that x,, with 
positive sign, is always applied to one of the inputs of each of amplifiers A,-Ag. 


With the resistances at the inputs of amplifiers Ay-Ag all equal, relay Re, operates for |x,| > 50 volts. 
Relay Reg operates with |x, | > 25 volts but, in the case when | x,|> 50 volts, relay Rey operates and its contact 


1Re, switches —50 volts to one of the inputs of each of amplifiers A,-Ags, thanks to which the contacts of relay 
Reg are opened. 
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If x, increases to the point when x,;— 50 > 25 volts, the total voltage on amplifier A, again reaches the 
magnitude at which relay Reg operates, The states of the relays in the circuits of amplifiers As, A, and Ag, as 
in the case just considered, depend on the value of xy and on the states of the contacts of the relays in the ampli- 
fiers’ input circuits, 


With such a circuit, the entire range of variation of x, (+ 100 volts) is divided into different intervals to 
each of which there corresponds a definite combination of the relay contacts. 


For example, let x, = 78 volts, In this case, contact 1Re, will be in state 2, since x, > 0, Contact 1Re,; 
will be in state 2 since xy > 50 volts, and contact 1Reg will also be in state 2 since xy— 50 > 25 volts, Since 
x,~ 50 — 25 < 12,5 volts, contact 1Res drops into state 1, 


Obviously, from what has been said, contact 1Re, drops into state 1, Relay Reg will be in the same state 
as relay Rey. No contacts of relay Reg appear on the portion of the MU circuit shown. Instead of different mag- 
nitudes of voltage (—50 volts, —25 volts, —12.5 volts, —6.25 volts, —3.125 volts) at the amplifier inputs, one 
might have applied equal voltages, but across unequal resistances, 


The second portion of the MU we are considering is shown in Fig. 3, Here, 1 and 2 are amplifiers and 3 
is a multiplier of some standard type. By means of contacts of relays Rey-Reg (Fig. 2), resistors are switched in- 
to the inputs of these blocks. In dependence on the sign of x;, contact 2Re, of relay Re, switches either voltage 
+X_ OF —X, to the input of amplifier 1 in such a way that the product y has the proper sign. 


Voltage xX, can be applied to the amplifier inputs across resistances of 2R, 4R,... , 32R. 


At the input of amplifier 2 a voltage of —100 volts is impressed across resistances of 2R, 4R,... , 32R 
and contacts of the relays mentioned, as well as voltage x, across resistance R, 


In the example given, with x, = 78 volts, voltage x, is applied to amplifier 1 via contacts 2Re, and 2Reg 
and resistances 2R and 4R, Simultaneously, a voltage of —100 volts is applied to amplifier 2 via contacts 3Re, 
and 3Re,. 


Voltage x, is applied to the input of this same amplifier. A resistance of 32R is connected between the 
input and the output of the amplifier. 


Consequently, the voltage at the amplifier's output is 


= 32.2, —100(F +f = 96 Vv. 


For X_ = 100 volts, the voltage at the output of block 3 is 
ug = 0.01 X 96 x 100 = 96 volts 
The voltage at the MU's output is 


y= pate 


The circuit given permits a decrease of the static error of an MU of standard type (block 3) by a factor 
of 32, 


If one of the MU's input voltages, say xy, varies stepwise from —100 volts to +100 volts, the transient 
response time in the combination-type MU will equal the sum of the operation times of relays Re» to Reg. 


A disadvantage of this variation of the circuit for the roigh portion is the voltage jump, small in magni~ 
tude and duration, at the MU's output at the moment when the relay contacts transfer, These jumps can be 
significantly decreased if polarized high-speed relays are used for the switching. 


A relay variation can be made without splashes if, during the transient response (for large deviations), 
the exact portion is "roughened." Such a device was tested and gave good results although, it is true, the circuit 
of the MU was somewhat complicated thereby. In the x, channel, these MU's have frequency characteristics 
which are the same as those for ordinary operational amplifiers. 
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Precision Multiplier Circuit with Reversible Selectors 


The rough portion of an MU circuit can be constructed by means of a two-bank reversible stepping switch, 
The circuit of such an MU is shown in Fig. 4. 


Two voltages, of opposite sign, are applied to the input of amplifier Ay: input voltage x, of the MU and 
the voltage from the brush By; of the switch (selector), This latter voltage is tapped off voltage divider D, via 
one of the commutator segments of the switch, 


Fig. 4 


The output voltage of amplifier Ay (Ax;) equals the difference of its input voltages amplified k times, 
If voltage Ax, attains its limiting value (for example, 100 volts), two single-contact polarized relays Re, and 
Re, operate which, in dependence on the sign of Ax,, lock in the winding of one of the electromagnets of the 
reversible switch (electromagnet SM, or SMg). In this case, the brush B, of the switch is moved to a neighbor- 
ing commutator segment, thanks to which voltage Ax, is decreased and becomes close to zero, 


If, as a result of a change in x,, voltage Ax, again attains the value of 100 volts, the switching described 
will again occur in the same order. 


Simultaneously with B;, brush B, steps along the commutator segments on the second bank of the switch. 
These segments are connected to voltage divider D,. At the inputs of this divider are impressed voltages +X, 
and —xX,, equal in magnitude to the second factor in the MU. Two voltages are impressed on output amplifier 
Az: one from brush B, and one from the MB block, which is an electronic multiplier of the usual type which 
multiplies voltages Ax, and xg. To obtain normal operation of the MB, it is generally necessary to supply at its 
input both the positive and the negative values of the factors Ax; and xg. In the circuit given, inversion of these 
voltages is effected by amplifiers Ay and Ay. 


As an example, let us consider a circuit in which the two-bank switch contains 40 commutator segments 
in each bank, We shall assume that the maximum value for the input and output voltages of the MU and MB 
equal 100 volts, We shall also assume that voltage v, = 100 volts, For this case, we shall determine the ratio 
of the impedances switched into the circuit, 


For the selected values of the parameters, the voltage between neighboring switch commutator segments 
equals 5 volts. Consequently, switching must occur for values of x; equal to 5, 10,15,..., 95 volts. Ampli- 
fier Ay must have a gain of k = 20, 


At the amplifier output we can set up a voltage divider consisting of resistors Ry and Rj so that, for suffi- 
ciently large values of the MU's input impedance, resistor Rs need not be too large, Amplifier A, must have 
equal gains for each of its inputs; the ratio Ryy/ Ry = 1/ 20, and ratio Ry/ R, = 1. 


As in the circuit considered earlier, in this circuit a small deviation of the relay operation voltage from 
the nominal value does not lead to an increased multiplication error if the scales of the input and output vol- 
tages of the ordinary MU slightly exceed 100 volts. 


For the same end, the operation voltages of relays Rey and Re are set to be slightly less than the scale 


/ 
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of the input voltage to the MB, With this, one also succeeds in eliminating the jumps in the MU'S output 
voltage and avoids the concomitant additional error in the MU, with the condition that the contacts between the 
commutator segments aud brushes B, and By are not disturbed. In the circuit given as an example, the error of 
the MU design is approximately 20 times lower than for an MU of the ordinary type. 


In the given MU circuit one may, instead of the reversible selector switch, use the following system. 


An MU Circuit with a Relay-~Type "Rough" Portion and a Thyrite "Exact" Portion 


The rough portion of the MU variant to be considered (Fig. 5) includes a row of symmetric triggers 
Try-Tr,, with different operating voltages. The number n of such triggers is determined by the admissible error 
of the MU, If, for example, the exact portion of the MU has an error of 0.5% of the output scale and the error 
of the combined MU must not exceed 0.05%, then the circuit must include 10 triggers with operating voltages 
of 0, 10, 20,... , 90 volts. 


An electromagnetic relay is connected in the anode circuit of the vacuum tube in each of the triggers. 
These relays switch the impedances at the inputs of amplifiers A; and Ag. The exact portion of the MU (enclosed 
by the dotted line) includes two thyrites T; and T;, diodes D, and Dz, and resistors, Its description is given in 
[5]. In order to decrease the number of relays and triggers, the inputs of the latter have only positive voltages 
applied to them, 


This operation is performed by trigger Try, which operates when x, = 0. For xy < 0, contact 1Re, of this 
trigger's relay switches the input of the MU to the output of amplifier A,, thus changing the sign of x. To as- 
sure that the output voltage has the correct sign, analogous switching is also carried out for input x, by means 
of contacts 2Re, and 3Re,. 


As in the previously considered circuits, the present circuit may, for a large number of problems, be con- 
sidered practically inertialess for the x, channel. If polarized relays, for example, type RP, are used for switch- 
ing, one may succeed in significantly increasing the rapidity of action of the x, channel as well, In processing, 
for example, a step function of arbitrary amplitude, the duration of the transient response will not exceed 3 to 
5 milliseconds, 


If it is necessary to obtain several products in which one of the factors is common (x4Xp, X4X3, « - +%4%n)> 
then the rough portion may be implemented as a common element for all the MU"s, It is convenient to use this 
method to obtain, for example, a power series ayx + agxX + asx’ +... + anx", with which one can frequently 
succeed, in computing devices, in reproducing a nonlinear function. 
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SUMMARY 


The circuits we have been considering for multipliers with “rough” and “exact™ portions have greater 
accuracy and speed of action than do existing circuits. They are simple to adjust, and the elements comprising 
them are simple. 


Experimental investigation of a number of circuits of multipliers constructed in accordance with the 
principles here considered were carried out as early as 1952, and gave good results, The static errors of such 
MU's did not exceed 0.01 to 0.02% of the output voltage scale. 


Subsequently these components were manufactured as parts of nonlinear electronic models, 


S. P, Onufriyuk and L, V. Medvedev took part in the adjustment and experimental investigation of the 
sample multipliers, 
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THE ANALYSIS AND SYNTHESIS OF CERTAIN DISCRETE CONTACTLESS CIRCUITS 


B. I. Rameev and Yu. A. Shreider 
( Penza and Moscow) 


The article treats of the code circuits, based on diodes, which are used in the 
Strela computer. 


Algebraic methods for analyzing and synthesizing such circuits are considered. 
The possible ways of constructing these circuits economically are discussed. Use is 
made of a,logical symbolism close to the ordinary Boolean algebra. 


1. Discrete Circuits 


It is quite well known that logical (Boolean) algebra has been successfully used in the analysis and synthesis 
of circuits based on electromagnetic relays. The possibility of such an application was first discussed by V. I. 
Shestakov [1] for the so-called parallel-series circuits (P-circuits), Subsequent works contained the systematic 
application of the apparatus of Boolean algebra to the study of the most general circuits constructed of electro- 
magnetic relays, Recently, the role of electromagnetic relays has been played, in many devices, by various 
contactless elements.* 


In the present paper there is developed an apparatus which permits the analysis and synthesis of diode cir- 
cuits to be carried out. This apparatus is a variant of the ordinary Boolean algebra, and not less simple than 
this latter. 


First, we must define what we shall mean by a discrete circuit. 


Let there be given some electrical circuit Z in which certain points are distinguished A, B,C,..., and 
let there exist for each of these points a succession of moments of time t,(A), tg(A), .. . t¢(B), 4(B),..., 
such that the electric potential at the given points at these moments of time assume only two values: a(t) and 
(ty), a(t,) and A(t,), . . . , and B(ty), etc. These points and moments of time are said to be characteristics 
of the given circuit. 


In what follows, we shall denote one of the two values of potential as 1 and the other as 0. Instead of the 
term “potential” we shall frequently use the term “voltage.” 


Such a definition of a discrete circuit is somewhat idealized. What it actually comes down to is this: 
that for studying the operation of a host of circuits it suffices to know only that the voltages at certain points of 
the circuits at certain times are, or are not, greater, or less, than certain levels. For this, a knowledge of the 
voltages at other points and at other moments of time is nonessential. Thus, the behavior of a discrete circuit 
is determined by the voltage distribution at the characteristic points of the circuit at the characteristic moments 
of time, 


We give here several examples. 


1) A digital bus in a serial computer, Here, each bus transmits binary numbers in the form of a temporal 
sequence of pulses, The characteristic point here is any point of the bus, and the characteristic moments of time 
are determined by the common system of machine synchronization, 


* The questions of synthesis and analysis of electronic circuits were discussed most completely in (2). 
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For example, let the series of voltage values on the bus at times ty, tg, ts, tg, ts and tg equal (in our con- 
ventional symbolism) 1, 0,1, 0,9, 1. This means that the bus is transmitting the binary number 101001. 


2) An electromagnetic relay (Fig. 1). The characteristic points A and B are shown on the sketch. The 
characteristic time corresponding to point B is defined by the condition that t(B) = t(A)+r, where Tr is the 
contact make~break time. 


Since the relay state is described in terms of make or break, it can be completely characterized by the 
voltage at point B being equal either to zero or to +E, 


Unless something is said to the contrary, we shall assume in what follows that we are considering only 
static circuits, i.e., we shall assume that the state of the basic circuit Z does not depend on time. 


A system of characteristic points of circuit Z will be called a complete system of inputs of Z if the vol- 
tages at all these points can be assigned independently and if the voltages at all the remaining characteristic 
points of the circuit are completely determined by the values of the voltages at the inpifts, For example, in 
the circuit of Fig. 2, points A and B form a complete system of inputs, 


For each characteristic point A of circuit Z one may form the proposition a: “The voltage at point A 
equals 1." The negation of this proposition is the proposition a: “The voltage at point A equals 0."* 


2. Structural and Characteristic Formulas 


Here we may pose the general problem of analyzing relay-type circuits, Let there be given a circuit Z 
and a combination of voltages at the inputs. It is necessary to determine what the voltage will be at any pre- 
viously designated characteristic point of the circuit. 


We shall denote the characteristic points by the letters A, B,C, ... , and will use the letters a, b,c,..., 
to denote the propositions: “the voltage at point A equals 1," “the voltage at point B equals 1". . . Then, by 
the definition of input, every proposition of the form “the voltage 2* characteristic point F equals 1" may be 
written in the form of a formula of the propositional calculus in terms of the propositions a, b,c, ... , for the 


corresponding inputs. 
A system of formulas of the propositional calculus 
f =f G@,6,c,...), 


g=g(a,b,c,...), (1) 
h=h(a,b,c,...), 


written for every characteristic point of the circuit Z being considered will be called the system of characteristic 
formulas of circuit Z. 


The problem of analysis of any discrete system can be thus formulated: given a circuit Z, to write its 
system of characteristic formulas. 


*For a time-dependent circuit these propositions are generally meaningless, and should be replaced by proposi- 
tions of the form a(t): “the voltage at point A at time t(A) equals 1." 
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The basic problem of synthesis of discrete systems can be formulated as follows: given a system of 
formulas of the propositional calculus which express the dependence of certain propositions f, g,h,..., on 
other propositions (1), to realize a circuit Z with inputs A, B,C, ... , and characteristic points F,G, H,..., 
for which system (1) would be part of the system of characteristic formulas. 


The example of a circuit built of electromagnetic relays shows that different circuits can be associated 
with one and the same set of characteristic formulas. 


We are thus led to the necessity of characterizing each circuit by some structural formula in such a man- 
ner that different circuits will correspond to different structural formulas, For this, we shall consider as identic 
those circuits which are comprised of identical elements identically connected, i.e., circuits which can be 
symbolically represented by one and the same sketch. 


We shall assume here that we know how to set each discrete circuit, lying in some chosen class of circuits, 
in correspondence with a system of structural formulas in the sense used above, Then, the basic problem of 
analysis is formulated in this way: for a given system of structural formulas, to determine the system of char- 
acteristic formulas of the circuit, The basic problem of synthesis is formulated as: given a system of character- 
istic formulas, to find all the systems of structural formulas corresponding to it. 


Of great interest are problems of the type: for a given system of characteristic formulas, to find structural 
formulas possessing certain extremal properties, for example, minimum number of elements in the circuit, mini- 
mum duration of the transient response when the voltage at the inputs is abruptly switched, optimum conductance. 
Extremal problems of this type are important in the applications of the theory, so that it is quite important to be 
able, by comparatively simple 1. ns, to determine from a system of structural formulas what the properties of 
the corresponding circuit are. 


It is shown in the works of V. I. Shestakov and of other authors that each circuit consisting of contact pairs 
connected in parallel or series (class P circuits, in the usual terminology) can be placed in correspondence with 
some formula of the propositional calculus, With this, the number of literal symbols entering into the circuit's 
structural formula equals the number of contact pairs. In this case, the characteristic formula of the network 
coincides with its structural formula, The essential result achieved here is that any formula of the propositional 
calculus can be realized as a structural (and, consequently, as a characteristic) formula of some relay-contact 
circuit of class P. This result showed the possibility, in theory, of synthesizing a circuit implementing any arbi- 
trary logical function, From this there came the rules for transforming relay-contact circuits with the aim of 
decreasing the number of elements. 


3. Diode Networks 


In this work we shall consider networks whose basic elements are diodes. Circuits consisting of such diode 
nets can be used in the most general cases.* This follows from the result, obtained below, that any character- 
istic formula can be realized by means of such networks. 


At this point we shall define what is to be under- 


pied r ae ; stood by a network built of diode circuits, A network 
TR: s * of this type consists of busses, connected by diode ele - 
is ments, and connecting voltage sources to them, In the 


sequel we shall assume that the busses have zero im- 
pedance and that the diodes are ideal. To each input 
terminal there corresponds a pair of input busses de- 
noted by (p, p). If there is a voltage on input p there 

a b will be a positive voltage +E on the first of these busses 
and zero voltage on the second, In the contrary case, 
there will be zero voltage on the first bus and +E on 
the second. In fact, both input busses are connected to 
the corresponding input terminal via a buffer element (switch), which may be implemented by a contact relay 
or an electronic tube, The input busses must satisfy the following requirements: 1) with no voltage the circuit 

is grounded; and, 2) no internal sources of voltage are connected to them, 


Fig. 3 


*In particular, the arithmetic unit of the Strela computer is based on such circuits. 
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The remaining busses of the network are called working or output busses. The output terminals can be 
connected directly to any of the working busses. 


To each point A on any of the working busses we may associate the proposition "A has voltage on it." 
Since the voltage on any bus can assume only two values, +E and 0, the proposition just stated (which we shall 
denote by the letter A), can be expressed by a formula of the propositional calculus in terms of propositions 
P.q4,..., Of the type: “input p has voltage on it.” 


The working busses are connected to the input busses by means of diodes which are always reversed with 
respect to the working busses. 


Let us consider the examples of networks shown in Fig. 3. 


In example a of Fig. 3, there will be high voltage on bus A only in the case when there is high voltage 
on both inputs P and Q, Indeed, in the other cases, if at least one of the P and Q busses had low potential on it, 
the voltage on the output would also be low due to the voltage drop across the resistance. Thus, the character- 
istic function of this network can be described by the equation A = pq. 


In example b of Fig. 3, the voltage on the output will be high if there is either voltage on input P or no 
voltage on input Q, Indeed, in this case there will be voltage at the output only when there is voltage on one 
of the horizontal busses T or S. But current cannot flow between busses T and S and, consequently, bus T will 
have voltage on it only when there is high voltage on input P, and bus S will have voltage on it only when there 
is voltage on input ©. Thus, the characteristic function of this network has the form: 


A=p+q. 


Various examples show that the elementary logical functions are realized by means of diode networks, It 
is not difficult to convince oneself of the validity of the following assertion. 


Theorem 1. For any formula of the propositional calculus* which contains propositions p,q,r,..., 
there exists a diode network each of whose inputs corresponds to one of the propositions p, q, r, such that its 
characteristic formula will be the initial formula. 


Proof. Every formula of the propositional calculus can be represented in the form of a sum of products of 
the basic propositions p, q,r, ... , or of their negations in such fashion that each product will contain either 
the proposition p (correspondingly for q,r, . . . ) or the proposition p and, moreover, will contain it just once.** 


We set up the network with inputs p,q,r,... 
that have working busses which realize the correspond- 
a ing products, and we connect these busses to the out~- 
ee LES puts via diodes, It is easily seen that the characteristic 

formula of this network is precisely the corresponding 
‘dl aii -? sum of constituent elements, As an example, Fig. 4 
, | 4 shows the network corresponding to the formula 


‘ 
a R 


Fig. 4 par + Par + 
It is apparent from the proof of Theorem 1 that 
for networks of the type considered, with n inputs, any characteristic formula can be realized with no more than 
(n + 1)2""* diodes, 


In order to implement a network with several outputs, with previously specified characteristic formulas 
for each output, it suffices to implement the system of networks, each of which has one output and which realizes 
one of the proposed characteristic formulas, aad then to joint the input busses and buffer elements of all these 
networks. From this follows Theorem 2. Let there be given n formulas which depend on the propositions p ,q,f e+ 
Then there exists a network (diode net) whose outputs have the original formulas as their characteristic formulas. 


* The propositional calculus is sometimes called the algebra of logic, or Boolean algebra. The formulas of the 
propositional calculus are sometimes called Boolean functions. 

* *In other words, each formula of the propositional calculus is decomposed into sums of the constituent units 
(cf. 
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4, Propositional Calculus and Operators 


The results of the previous section showed the possibility, in principle, of realizing any logical function 
by means of diode nets. These results, however, give no means of obtaining networks which are the most desir- 
able, either in terms of the number of elements required or in terms of their physical properties. 


Fig. 5 Fig. 6 


For this purpose we shall now develop a formal 
apparatus which will permit the structural formulas of 
any diode net to be determined. 


L& In contradistinction to the propositional calculus 


| as Ordinarily applied, we shall be dealing, in what 
owe ro follows, with quantities of a dual nature: with proposi- 
| ay | tions and with operators. This bifurcation corresponds 
R a ? to the fact that the networks of the type considered 
Fig. 7 consist of two forms of busses: input and working. For 


this, each input bus P breaks up into two input busses 
P and P in such fashion that if there is voltage on input P then bus P will have high voltage and P will have low 
voltage; in the contrary case, the voltage distribution on busses P and P will be reversed. 


We introduce the following concept. Bus S will be called subordinate to bus T if, in the corresponding 
network it is possible to distinguish a sequence of busses, Ry, Rg, . . . » Rn, such that the connection of bus T 
with bus Ry is implemented by a diode which conducts current from bus T to bus Ry, busses Ry and S are con- 
nected by a diode conducting current from Rp to S; finally, any two busses Ry and Ry +4 are connected by a 
diode conducting current in the direction from Ry to Ry 44. 


For example, in Fig. 5, bus G is subordinate to bus A and also to busses B and F, while bus H is subordinate 
to busses D and E, 


To each working bus § there is associated the proposition “bus S has high voltage." We will denote this 
proposition by the letters. To each input P there is associated the operator (proposition), “P has high voltage,” 
which we shall denote by the letter p. In the sequel we shall find it convenient to assurne that the operator p 
is associated with input bus P and operator p is associated with input bus P, Each voltage source is associated 
with the tautologically true proposition E. 


In a general form, the rules by which one can form the proposition S$, corresponding to some working bus 
S,consist of the following. 


1) We consider all working busses T, R,. . . , which are connected to bus S by diodes which conduct 
current toward bus S, and we form the sum of the propositions T+R+...,. In particular, if a voltage source 
is connected to bus S, then the tautology E should be added to the sum, 


2) We next consider all input busses P, Q, R, each connected to one of the busses to which S is subordinate, 
and we apply successively the operator propositions p,q, ..., to the proposition T+ Ri+...,. As the result 
we obtain the expression 


S=pq...(T+R+...) 
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3) Operators will bear upper and lower indices, The upper index on an operator p is one less than the 
number of diodes through which current must pass in order to flow between the working bus under consideration 
and input bus P, The lower index gives the ordinal number of the diode connected input bus P with the working 
bus. Each vertical bus will be numbered independently of the other, 


Consider the example of Fig. 6. In accordance with the rules formulated, we have the sequence 


A=Kk |-T -+- B. 

R=- E, 

T= pre, 

s,(D-\ C)= (qi? ry py). 
Thus, 
A = srpqE + rsqE + s(qpr + rpE). 
The operator calculus developed here allows each single-output network of the type considered to be un- 
ambiguously put in correspondence with a structural formula. One may also easily convince oneself that the 


converse is true: each formula of the operator calculus can be interpreted as the structural formula of some 
network (diode net) with a single output. 


We are now in a position to formulate the following rules which permit the characteristic formula of a 
network to be derived from its structural formula. 


1) By using the distributive law, multiply out any parentheses in the structural formula, 
2) Delete every monomial which does not contain E, 

-3) Delete E from every monomial in which it appears, 
4) Delete all indices from the operators, 


5) The expression thus obtained is to be interpreted as a formula of the propositional calculus, the 
operators being considered as propositions. 


Thus, the characteristic formula of the network shown in Fig. 6 has the form: 
A = srqpE + TsqpE + s(qpr + rpE) = srpq + stpq + rps. 


A rule can be stated for calculating the total number of diodes from the system of structural formulas. 
This rule reduces to the following: 


1) Count the total number of terms v in each structural formula. 


2) For each operator p, choose the maximum value of the lower index, and sum these values for all the 
operators, The result is denoted by p. 


3) For each term in each structural formula, choose the maximum value of the upper index, and sum 
these values for all the terms, The result is denoted by A. 


The total number of diodes in the circuits equals the sum ). +P +v. 
Example. Consider the system of structural formulas for some network: 
| pyr PE. 
Letv =8,A = 2, 


Figure 7 shows the corresponding network. 


5. Simplification of the Structural Formulas 


We now consider transformations of the networks for which the system of characteristic formulas remains 
invariant. 


— 


4, Propositional Calculus and Operators 


The results of the previous section showed the possibility, in principle, of realizing any logical function 
by means of diode nets, These results, however, give no means of obtaining networks which are the most desir- 
able, either in terms of the number of elements required or in terms of their physical properties. 


= 


Fig, 6 


For this purpose we shall now develop a formal 


pa ~ apparatus which will permit the structural formulas of 
any diode net to be determined, 
¥ ez b, In contradistinction to the propositional calculus 
ron as ordinarily applied, we shall be dealing, in what 
Oo > follows, with quantities of a dual nature: with proposi- 
| f tions and with operators. This bifurcation corresponds 


to the fact that the networks of the type considered 
consist of two forms of busses: input and working. For 
this, each input bus P breaks up into two input busses 
P and P in such fashion that if there is voltage on input P then bus P will have high voltagg and P will have low 
voltage; in the contrary case, the voltage distribution on busses P and P will be reversed. 


We introduce the following concept. Bus S will be called subordinate to bus T if, in the corresponding 
network it is possible to distinguish a sequence of busses, Ry, Rg, .. . , Rn, such that the connection of bus T 
with bus R, is implemented by a diode which conducts current from bus T to bus Ry, busses Ry and S are con- 
nected by a diode conducting current from Rp to S; finally, any two busses Ry, and Ry +4 are connected by a 
diode conducting current in the direction from Ry to Ry 44. 


For example, in Fig. 5, bus G is subordinate to bus A and also to busses B and F, while bus H is subordinate 
to busses D and E. 


Fig, 7 


To each working bus S there is associated the proposition “bus S has high voltage." We will denote this 
proposition by the letters. To each input P there is associated the operator (proposition), “P has high voltage,” 
which we shall denote by the letter p. In the sequel we shall find it convenient to assume that the operator p 
is associated with input bus P and operator p is associated with input bus P, Each voltage source is associated 
with the tautologically true proposition E. 


In a general form, the rules by which one can form the proposition S, corresponding to some working bus 
S,consist of the following. 


1) We consider all working busses T,R, . .. , which are connected to bus S by diodes which conduct 
current toward bus S, and we form the sum of the propositions T+R+...,. In particular, if a voltage source 
is connected to bus S, then the tautology E should be added to the sum, 


2) We next consider all input busses P, Q, R, each connected to one of the busses to which S is subordinate, 
and we apply successively the operator propositions p,q, ... , to the proposition T+ Ri+...,. As the result 
we obtain the expression 


S=pq..(T+R+...). 
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3) Operators will bear upper and lower indices, The upper index on an operator p is one less than the 
number of diodes through which current must pass in order to flow between the working bus under consideration 
and input bus P, The lower index gives the ordinal number of the diode connected input bus P with the working 
bus. Each vertical bus will be numbered independently of the other. 


Consider the example of Fig. 6. In accordance with the rules formulated, we have the sequence 


|-T -+- B. 
R= E, 
T= rng pre, 


Be C) = (qi pry! ry 
Thus, 


A = srpgE + rsqE + s(qpr + rpE). 
The operator calculus developed here allows each single-output network of the type considered to be un- 
ambiguously put in correspondence with a structural formula, One may also easily convince oneself that the 


converse is true; each formula of the operator calculus can be interpreted as the structural formula of some 
network (diode net) with a single output. 


We are now in a position to formulate the following rules which permit the characteristic formula of a 
network to be derived from its structural formula. 


1) By using the distributive law, multiply out any parentheses in the structural formula, 
2) Delete every monomial which does not contain E. 

3) Delete E from every monomial in which it appears. 

4) Delete all indices from the operators. 


5) The expression thus obtained is to be interpreted as a formula of the propositional calculus, the 
operators being considered as propositions, 


‘ Thus, the characteristic formula of the network shown in Fig. 6 has the form: 
A = stqpE + TsqpE + s(qpr + rpE) = srpq + stpq + rps. 


A rule can be stated for calculating the total number of diodes from the system of structural formulas. 
This rule reduces to the following: 


1) Count the total number of terms v in each structural formula. 


2) For each operator p, choose the maximum value of the lower index, and sum these values for all the 
operators. The result is denoted by py. 


3) For each term in each structural formula, choose the maximum value of the upper index, and sum 
these values for all the terms. The result is denoted by A. 


The total number of diodes in the circuits equals the sum A +H +. 
Example. Consider the system of structural formulas for some network: 


AN | palm 
| 


Letv =8,A =2, 
Figure 7 shows the corresponding network. 


5. Simplification of the Structural Formulas 


We now consider transformations of the networks for which the system of characteristic formulas remains 
invariant. 


| 


ing 


For this, we formulate the following rules for operating with propositions and operators, 
1) Addition of propositions is commutative, so that 
A+B=B+A, 


2) Propositions can be multiplied by operators, the result being a proposition and the following relation- 
ship being valid: 


= qpA. 
3) The sum of the tautology with any proposition is again a tautology: 
E+A=E, 
4) The distributive law is valid: 
p(A + B) = pA + pB, 


It is clear that with such transformations of the structural formulas there is no change in the characteristic 


formulas, and that with each such transformation one obtains a structural formula which corresponds to some 
actual network, 


Fig 10 Fig. 11 


Let us now consider in more detail transformation of networks with one output, We note immediately that 
in this case the logical apparatus can be simplified by dropping the indices on the operator propositions in the 
formulas, In fact, the necessity for having upper indices was due only to the fact that some input bus could have 
been connected to two working busses by means of one diode. In such a case it was necessary that these working 


busses be separated by diodes, and the purpose of the upper index was just to allow these separating diodes to be 
taken into account, 


As an example, we take Fig. 8. Such a method for connecting busses turns out to be possible when there 
is a common factor in several terms. However, it is clear from this same Fig. 8 that by transforming this net- 
work to the form shown in Fig. 9 we can save two diodes by eliminating the separating diodes. 


The example just considered shows that in the general case it is possible to avoid the use of isolating 
(separating) diodes. But in this case it is no longer necessary to retain either the upper or the lower indices, 
since each diode can be used only for connecting one factor in the structural formulas to an input bus. Here, 
the rule for computing the number of diodes realizing the given structural formulas can be formulated as follows: 
the total number of diodes is equal to the number of terms in the structural formulas plus the number of all the 


operators entering into these formulas. With this, each operator is counted as many times as it enters into the 
formulas. 
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Example. The number of diodes realizing the structural formula: 
A = pqE + p(qrE + qrE), 
equals 4+ 7 = 11. 
This formula can be transformed to the form 
A = pgE + parE + parE. 
In this case, the number of diodes equals 3+ 8 = 11. The corresponding networks are shown in Fig. 10, 
where Fig. 10a corresponds to the first of these formulas and Fig. 10b corresponds to the second. 


Thus, ...nplification of structural formulas should be effected by factoring out common factors in the 
terms, In so doing, however, the total number of terms should not be too greatly increased, since such an in- 
crease can lead to a deterioration of the frequency characteristics of the circuit, to the appearance of false cir- 
cuits, and to disruptions of circuit operations due to the finite magnitude of the diodes’ back resistance. 


The same rules for simplification remain valid also for cixcuits with several outputs. In addition, however, 
certain special simplifications are possible if several outputs have common factors or terms. 


As an example, consider a one~bit binary adder with three inputs P, Q and R, and two outputs: © is the 
sum and II is the carry, It is easily shown that the characteristic formulas for such an adder have the forms: 
Part par + par, [| = par + par + par + par. (2) 
We can replace Formulas (2) by equivalent formulas of the logic of algebra 
>} = + ar) + + 97), ]] + (ra + 79). (3) 
By introducing the corresponding propositions one can, from Formulas (3), construct the system of structural 


formulas: 
A=rqE, B=rqE, G=p(A+B), 
C=rqE, D=rqk, H=pC+D), 
E=rqE, F=rqE, K= F), 
L=rqE, 


The circuit of Fig. 11, consisting of 27 diodes, corresponds to these structural formulas. 
This circuit can be simplified by using common working busses for the outputs £ and I. 
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CHAIN MATRICES AND SUBSTITUTION CIRCUITS 
FOR ELECTROMAGNETIC~-~TYPE ELECTROMECHANICAL TRANSFORMERS 


A. Ya. Mel'nichuk 


(Moscow ) 


Different variations of four-terminal substitution circuits for electromagnetic- 
type electromechanical transformers are considered. The possible methods of esti- 
mating the influence of “negative rigidity” on transformer operation are discussed. 


Electromagnetic-type electromechanical transformers are widely used in the apparatus of automation and 
remote control, In particular, polarized transformers with linear characteristics are used in tuned relays, electro- 
mechanical filters [1,2,3], tuning-fork frequency regulators [4,5], vibration gyroscopes, vibration transducers, etc. 


We shall consider the processes occurring both in the electrical and in the mechanical portions of the 
transformer, since the transformer is most effectively considered as an electromagnetic four-terminal network, 
described by the corresponding matrices [2, 6-9]. 


In [8] there was proposed a matrix system for different types of transformers which was based on the analogy 
of “force-voltage." However, it was shown in [7] that in order for the principle of duality® to hold, it is neces- 
sary to use the analogy of "force-current” in setting up the matrices of the electromagnetic transformer, More- 
over, in the matrix system used in [8], phenomena attendent on the transformer processes were not taken into 
account, specifically, the occurrence of a negative component in the expression for the mechanical rigidity of 
the electromagnetic transformer (described, for example, in [1]). 


Below we shall compare various methods of matrix and schematic presentation of polarized electromag- 
netic transformers, and shall give the possible methods for estimating the influence of “negative rigidity” on 
their operation. 


1, Chain Matrices of an Electromagnetic Transformer 


Figure 1 gives the schematic representation of an electromechanical four-terminal network which cor- 
responds to an electromagnetic transformer with a movable core, The directions chosen for the voltage U, cur- 
rent I, force F and velocity V are nominal, Four-terminal network a transforms electrical energy into mechani~ 
cal energy, while network b transforms mechanical energy into electrical. A series (chain-like) connection of 
such four-terminal networks is possible, both on the electrical and on the mechanical side (for example, in 
electromechanical frequency filters). 


The system of equations which describes the processes in an ideal electromagnetic polarized transformer 
working with small sinusoidal oscillations (account not being taken of magnetic leakage and losses in the elec- 
trical portion) has the following form: 


U = + pV, 


*It is necessary that the duality principle hold in order that the electromechanical transformer might be re- 
Placed by the purely electrical passive circuit equivalent in form [7]. 
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Here, Zyp is the mechanical resistance of the moving portion of the transformer with the electrical circuit 
open, 8 is the coefficient of the electromechanical connection, Lf is the inductance of the braked transformer 
(cf. Appendix). 


The force F engendered by the total magnetic flux @ of the transformer is defined by the expression 


(2) 
where d is the constant polarizing magnetic flux, ¢ is the variable magnetic flux generated by the coil current i 
(ordinarily, @ >> y), lg is the magnetic permeability, and S is the cross section area of the effective air gap. 


With no alternating current in the transformer coil, the change in force can be expressed in the form (cf, 
Appendix): 


dF = 
(3) 
where A, is the constant, and 6 the variable, component of the air gap width. 


One of the effects of the polarizing field is the engendering of an additional mechanical force (3), acting 
upon the armature in the direction away from its equilibrium position, This effect leads to the appearance of a 
negative component of mechanical rigidity of the transformer, this component equalling 


(4) 


It follows that the mechanical impedance will consist of two parts: 


where Zr, is the inherent mechanical impedance, and Z7, is the impedance of the "negative rigidity.” 


After substitution of expression (5) for Z,, in (1) and some regrouping of terms, we obtain a new system 
of equations which, in matrix notation, has the following form: 


oh, —™ j 
] hy B jolh V 
P 
The first factor in the right member of (6) is the chain matrix of the transformer. The rationale for this 


name is that, for a series connection of transformers (in a “chain"), the chain matrix of the whole series equals 
the product of the chain matrices of the individual transformers. 


(6) 


b 
Fig. 1 
| 
dF o? 
(5) 
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2. Circuit for a Substitution Transformer 


In substitution circuits it is effective to connect the corresponding electrical and mechanical elements on 
the electrical and mechanical sides of the transformer, Between these elements is to be found the specific four- 
terminal network which relates the electrical and mechanical quantities. 


Fig. 2. Substitution circuit of a transformer with a four- 
terminal network for an electromechanical connection 


depending on frequency. 


Fig. 3. Substitution circuit of a transformer with a four- 
terminal network for an electrical connection not depend- 


ing on frequency. 


The substitution circuit can be obtained by decomposing the chain matrix of (6) into its constituent 
matrices. One possible decomposition gives the following matrix product: 


ol 


According to the theory of four-terminal networks, the circuit of Fig. 2 corresponds to expression (7), The 
circuit given has the disadvantage that the transformation of the mechanical element to an electrical quantity, 
and vice versa, depends on frequency (since the factor w occurs in the matrix of the contenenems network of 


the electromechanical connection). 


Another method of decomposing the chain matrix in (6) gives the following matrix product: 
1 jol,, B 1 0 


x x 1 x (8) 
0 1 Jo 1 0 
Figure 3 corresponds to expression (8), Its advantage is that the mechanical quantity is transformed to an 
electrical quantity, and vice versa, independently of frequency, since the matrix of the four-terminal network 
of the electromechanical connection does not contain w. Moreover, the circuit in Fig. 3 allows the effect of 
"negative rigidity” to be estimated graphically. 


The quantity which is the reciprocal of the rigidity K of the mechanical system is transformed, on the 
electrical side, to an inductance 


1 


where w_ is the number of turns in the transformer's coil and I is the current necessary to give rise to the polariz- 
ing flux &. 


| 
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g x x (7) 
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Fig. 4 


If K is sufficiently large, then Lf, >> Lx. Consequently, the admittance jwCy, = j/ wL}, in comparison 
with the admittance 1/ jwLy, can be neglected; in this case, the capacitance C, will not appear in the curcuit 
of Fig. 3. 


In many cases the effect of “negative rigidity" must be taken into account, In particular, it can be used 
for the fine tuning of electromechanical filters and tuned relays, used in remote control frequency devices [1,3] 
and in the adjusting of tuning-fork frequency regulators [5], etc. 


As an example, we may take the substitution circuit for an electromechanical band-pass filter with flexibly 
connected vibrators which was developed for frequency telecontrol devices. Fig. 4a shows the equivalent electro- 
mechanical circuit, and Fig. 4b shows the electrical circuit which replaces such a filter. In the given case, the 
circuit parameters had the following values: 


30? a 
Lx, = & = = L,=0,9 


0.031 h, 


(m_ is the equivalent mass of the vibrator, N = newton), 


With this, the natural frequency of the tuned vibrator, corresponding to the upper boundary of the filter's 
pass band, equals 


— 


1 1 
j= [,0 =~ 1180 cps. 
With the vibrators in the magnetic field of the transformer filter with a small air gap (0.2 mm), the effect 
of “negative rigidity” (element C, in the circuits of Fig. 4) appeared. The magnitude of C; in the case given 
equals 


Cy = 0 02x10 Ff 


wo? b 


Due to the effect of C; the natural frequency of the vibrators must be reduced to the value 


which actually occurs, 


When account is taken of the voltage drop across the coil's ohmic impedance and of the inductance leak- 
age, as well as the active component of the input current, related to the loss in the iron (core), the transformer 
substitution circuit becomes somewhat more complicated (Fig. 5). Here are appended the element L},, the in- 
ductance leakage, the winding impedance Rj, and also the loss impedance Rj in the iron, All the previous dis- 
cussion applies also to this circuit, starting with terminals 1 and 2. 
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Fig. 5. Transformer substitution circuit with dissipation 
in the electrical portion accounted for. 


CONCLUSIONS 


Chain matrices and substitution circuits for electromagnetig circuits can have different forms, depending 
on the type of analogy employed and on the method of decomposing the matrices. The duality principle holds 
if the analogy “force-current” is used, For this, the most convenient substitution circuit is that in which the 
matrix of the four-terminal electromechanical connection does not contain a factor of frequency. In this cir- 
cuit, the transformation of electrical element to mechanical element, and the converse transformation, occurs 
independently of frequency. Moreover, this circuit allows the effect of “negative rigidity” on transformer opera- 
tion to be estimated graphically. This latter effect can be used for the regulation of frequency tuning in reson- 
ant systems, 


The author wishes to express his gratitude to P. I, Zubkov for his valuable comments on the questions 
considered here. 


Appendix 
1. The total differential force acting on the moving portion of the electromagnetic transformer equals 


OF OF OF a® oF a® 


The first term of this expression defines the variable force created by the coil's alternating current i, The 
value of the second term is found from the definition of its factors for i = 0: 


From the law of total current 


it follows that (for Ay >> 5) 


If we take the above expressions into account, we obtain 


4 


2. L} is the coefficient of proportionality between the increment of effective flux linkage and the incre~ 
ment of current. Since, for a braked transformer, we have the relationships 


= dy =9, 9 = then 


Thus, Lj, is the inductance of the braked transformer without account being taken of the flux leakage. 
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3. The value of the coefficient 6 is determined from the expression for the emf e directed to the trans- 
former's opened coil (for i = 0) and due to the oscillation of the armature with speed y: 


This expression can also be obtained from the dependence between current and force for a braked trans- 
former: 


_ oF ,,_ a® wo. 
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ON THE POSSIBILITY OF CERTAIN TYPES OF OSCILLATIONS 
IN SAMPLED-DATA CONTROL SYSTEMS 


Shao Da-Chuan 
(Peking) 


1. Posing of the Problem 


We shall consider an automatic control system containing a linear portion and a pulse (sampled-data) 
element (Fig. 1), As shown in Fig, 2, the system is characterized by the repetition period T,, the pulse ampli- 
tude ksign x(jT; )and length Tj = g|x(jTs)|.* 


yt), 20 


Pulse 
element 


U 


Fig. 1 Fig. 2 


We shall assume that the transfer function of the linear portion of the system has the form: 


G(s)= = 


where Red; < 0 and the gain is taken equal to unity; in other cases it may be taken equal to a constant k, We 
now consider the possibility of tl.ere being symmetric oscillations with period 2T,;. The presence of the integrat- 
ing term makes the existence of asymmetric auto-oscillations impossible, We choose the origin (t = 0) in such 
fashion that the input is positive (Fig. 3), and then 


It suffices to consider the process during the half-period Ts. 


2. Periodic Solution 


We first transform the equation of the linear portion 
pQpx = 


*If g|x(iT,)| > Ts, then Tj = Ts. 


| 
| 
7 | 
| 
1 = gx(0). (1) | 


where p = d/dt, to the normal form [1]: 


a (3) 
Then 
Figure 3 shows the form of y(t) when there are oscillations. With this, 
y =k for0<t< Ty, y =0 for Ty<t< Ts. 
Moreover, €; is continuous and € ;(0) = —€ ; (Ts). 
Taking these conditions into account, we present the solution of Eq. (3) in the form 
— O<t<T;: 
Ver, 
MMm<t<tT, 
k i+e t (5) 
| 
where 
kT 
£0) 
Ag 
k k 
tag 
and 
mets 
i— 
Fig. 3 +2 O) 4 (6) 
From (1) and (6) we obtain 
n 
4 1 1 i—e 


If Eq. (7) has a real and positive solution T,>0, then this means that the form of oscillation considered is 
possible, while the contrary case means that such an oscillat.on is impossible. 


We obtain the complete solution by substituting T, in (5). It is possible to estimate the relative magnitude 
of Ty since x(0) = T;/g. For T; > Ts one should set Ty = Ts. With this, 


rT, 


i=? 


3. Stability of the Oscillations 


We shall assume that there is an oscillation, We now consider a small disturbance €4 of the variable x. 
If this perturbation tends to zero as t—> oo, then the oscillation is stable ("in the small"), If we denote the 
periodic solution for €; by €f and the perburbation by 6£;, we then have 


Since £4 and €j are solutions of Eq, (3), the equations in the variations will take the form: 


=0, (i = 2,...,), (9) 
= (0) = const, 
BE, = BE, ..., (10) 


Equations (10) are valid in the interval 0 < t < Tj, where 


=1 


where 
or a n n 
T= © + 88, ()] = GF; 0) + =T, + 
i=1 i=1 i=1 
where 
= Py GE, (0). (12) 
Fort = Tj 


E,(T,) = + 88, (T)). 


But in the interval Tj, < t < T;, y = 0 and 
= (7) + 88, = + 88, (0), 


For t = T, we get 

E, = + 38, (0), 
A(T 


Since all the 6 €; are very small, the terms of Eq. (13) which depend on Tj can be computed approximately. 
Thus, for example, 
= (Fp » 


but in the interval 0 < t< Tj, a 


| 
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= & + BT yy = + Ty). (14) 


Similarly, we get 
(15) 
Substituting (14) and (15) in (13) we obtain 
(T,) (T,) + + (0), 
or 
(7) — 86, (0) = 
(16) 


BE (T,)— BF, 
(i) 


We now consider the process in the interval T, < t < 2T,;. After t = T, the pulse length equals Tj = 
= g|x(T;)|= —g(x)(T;), since x*(T,) and, consequently, x(T,) are negative. Letting Tj = T; +5 T(g), we get 


n 
= — > (T,). 
i=1 


If we wish to retain Formula (12), i.e., to set 


yy) = G,3E; (T,), 


then we must have a minus sign in the right members of Eq. (16). Moreover, with this,y = —k and the quantity k 
in the right members of Eq. (16) must be replaced by —k. As is obvious, Formulas (16) do not change. In 


general, 
BE, (nT,) — [(n — 1) T,] = ke 


8T = + G,3&, [(n — 1) T,}. (17) 
=) 


Equations (17) are difference equations with constant coefficients. Let the solution of (17) have the form 
(nT) = Mje", 8T ny) = No". 
Substituting this in (17) we obtain 


M, a— kNe tT (18) 


N= ¥' 
(19) 
From (18) and (19), taking (11) into account, we get 
{ 


1 


(20) 
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The stability conditions reduce to the requirement that all the roots p; of (20) satisfy the condition 


| pil<1. 
We now consider the case when gx(0) > T,. With this, Ty = T; = const, Here, 


(0) = GE, (0) = gz* (0) + 0) >7,. 


i=] i=] 


Since all the 6€; are very small, the pulse lengths, as before, equal T, and we may write 


(T,) = + 88, (T,) = & (T,) + ©), 
= & (T,) + (T,) = (T,) + 88, 


+ G88, |>7,, 
i=] 


i.e,, the length of the second pulse, as before, equals T,. In addition to 6 € ;(0), be (oye also decreases and, 
therefore, after a sufficiently long time, x(t) = x*(t) +[1/g)(G45 € 3(0)]. If we assume that the perturbations 

are random, we may consider that some perturbations increase the deviations and some decrease them, There- 
fore, after a k-fold action of the perturbations, DG,6 € 4, (0) will also be small, from whence we may determine 


that the oscillation is stable. 


4. Decomposition of Parameter Space Into Regions 


: We consider the simplest problem in which 


Q(s) =s+a. 
With this, T, satisfies the equation 


Equation (20) is written in the form 


i 
| (22) 
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We now introduce relative quantities, and the notation: 
O@=aT,, a=T,/T,, Y=gkT,, m= l= a8, 


and then (21) is replaced by 


(23) 


By considering / as a parameter, we can construct the curve showing the dependence of @/y on m (Fig. 4), 
On Fig. 4, curve A corresponds to the case when e! =m, or &@ =1, The roots of Eq. (22) are the roots of the 
equation 


The condition for stable oscillations hasthe form: f(0)< 1, f(+ 1)> 0. 


By taking (23) into account one can show that f(+ 1) is always larger than zero and that the limits of 
stable oscillations are determined from the equality 


(0) =e [1 + — |= 1. 


If we cancel y/@ from (23) we obtain 
lm? — 2 (e! — 1) (4 + I) m + [2(e' — 1) (1 =0. (24) 


If we are given 1 , we can compute m from (24), We thus obtain curve B (Fig. 4). We can decompose 
both the (7/6 , m) and the (y , 8 ) planes into two subplanes each. Subplane I corresponds, in each case, to 
stable oscillations, and subplane II to unstable oscillations. We might mention that if g is very small, the system 
corresponds to Tsypkin's first type of sampled-data system in which the limit of stability is defined by the curve 
1 =0 on Fig, 4. 


LITERATURE CITED 
A. 1. Lur’e, Certain Nonlinear Problems in Automatic Control Theory [in Russian] (Gostekhizdat,1951). 
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AUTO-OSCILLATIONS IN A SINGLE-LOOP AUTOMATIC CONTROL SYSTEM 
CONTAINING TWO SYMMETRIC RELAYS* 


Tu Syui-Yan" and Tei Lui-Vy 


(Peking) 


By the use of Ya. Z, Tsypkin's method [1], the exact equations are derived for 
the determination of the symmetric states in a single-loop circuit which contains 

two symmetric relays. It is shown that, if one assumes that the relays are separated 
by harmonic filters, then the exact method leads to the same results as the approxi- 
mate method based on harmonic balance, which is used by E. M, Nadzhafov (2). 


Ya. Z, Tsypkin [1], in considering the exact determination of the symmetric auto-oscillations in a single- 
loop system containing one relay, showed that if one assures that the linear part of the system is an ideal filter, 
then the exact methods for determining the periodic regimens (but not the conditions for their stability) lead to 

the same results as the harmonic balance method. 


In the present work, analogous results are obtained for a single-loop system containing two symmetric 
ideal relays. For this, the exact solution is derived by the method of Ya, Z, Tsypkin [1], but the harmonic 


balance method for the system with two nonlinearities is used in the form given by E. M. Nadzhafov [2).** 


1, 


The Exact Solution 
We shall consider the single-loop system (Fig. 1) described by the equations 


yi =f (21), Yo =f 
t2a=Wilp)y, t1=— 2, (1) 


where f(x,) and (x2) are ideal relay characteristics (Fig. 2) and W,(p) and W,(p) are the transfer functions of 
the linear portions of the system (Fig. 1). 


First % 


Fig. 1 


* This work was carried out under the direction of M, A. Aizerman, 
**The exact solution for the problem of symmetric oscillations of the system containing two ideal relays was 
obtained by Yu. I, Neimark [3], using other methods, 
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“ne We are interested in the symmetric oscillations and, following 
dite Ya. Z. Tsypkin’s method, we assume that the form of the oscillations of 
wid the coordinates y, and y is as shown in Fig. 3. Then the functions y,(t) 
and y,(t) can be expressed by Fourier series 


sin (2m — 1) wt, 


v(t) =< sin (2m — 1) w(t — 1), 


Fig. 2 


and the functions x(t) and xs(t) by the series 


> [Wal — gin — 1) ot + 


+ arg W; [i (2m — 1) @}}, (3) 
x(t) = —%(t)= 
sin {(2m — 1) (¢ — t) + 
: + arg W, [i (2m— 1) 
-* > | W, [i (2m — 1) | cos — 1) wt + arg W, [i (2m —1) 
m=1 (4) 


= > |W [i (2m—1) «|| 0s {(2m—1) (t—t) + arg We [i (2m—1)o]}° 


™m=1 


1 


Fig. 3 


In order that series (2) and (3) satisfy Eq. (1), the following conditions must hold: 
1, 2, (0) =0, 22 (t) = 0, 
2.21(0)>0, >0. 

It must be verified, moreover, that 


3. 21(t)>0 for 


(t)>0 for (7-=). 


The first and second conditions lead to the following equalities and inequalities: 


(2) 
| m=) 
/* 
nts \ 
' f \ / 
| 
84 


8 
ns of 


Gm — ell sin — 1) +) + arg Walt (2m —1) = 0, (5) 


4 
=»! 


1) ©} sin ((2m — 1) wt + arg W; fi (2m — 1) wh) = 0, (6) 
m=1 
eo < 
Dy | Walé (2m — 1) | cos ((2m — 1) (— +) + arg Ws [i (2m — 1) <0, 


m=1 (8) 


The solutions w» and Ty, of Eqs, (5) and (6), which satisfy inequalities (7) and (8), define the functions 
x(t) and x(t), which satisfy conditions 1 and 2, If, with this, conditions 3 hold, then the system under con- 


sideration will have a closed trajectory in phase space; if this is a stable limiting cycle, then the sought auto- 
oscillation of frequency w, can arise in the system. 


Im/, tml, 


¢& 
Fig. 4 


If we introduce the following “relay system characteristics” we can present the equalities and inequalities 
in (5)-(8) in a form more convenient for solution: 


= — 24, (0)— i210), — (x) — ina (9), 


By introducing the real and imaginary frequency characteristics Uj(w) and Vj(w) of the linear portions of 
the system 


W;; (iw) = U;(o) + iV;() (j= 1,2), (9) 
we may write conditions (5)}(8) in the form 


Im 1; (@, *)= 4 3 — [— Us (@,,) sin @,,t + V2 (@,,) cos @,,t] = 0, 
m=1 
Re J, (@, = 4 [Us (@,,) C08 @,,t + Ve Sin w,,t] <0, (11) 
m=1 
Im 5) = (mq) + Vi (Cp) £08 = 0, (12) 
m=1 
Re Iz (@, t) = [U1 (@,,) @,,t — V3 (w,,) Sin w,,t] <0, (13) 


m=1 


where w,, = (2m— 1)w. 


Being given various value of T , one can construct the graphs of 1,(w,r ) and I,(w,r ) (Fig. 4), The points 
in which the curves I,(w,r ) intersect the negative real semiaxis correspond to a series of values of rf and w; the 
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yx(t) 
) 
£2 — 1) cos ((2m — 1) + arg Ws [i (2m — 1) > 0. 
| 
| 
= | 


analogous points for the curves I,(w,T) correspond to 
another T ,w series. If we now construct the curves rT = 
= f;(w) and t = f,(w) on the f ,w plane, we obtain the 
desired Ty, wy as the coordinates of their point of inter- 
section (Fig. 5). 


2. The Approximate Solution, and its 
Comparison with the Exact One 
ff (0) We assume here that both linear portions of the 
r system are ideal filters which pass only the first harmonics 
Fig. 5 of the periodic functions y,(t) and y,(t). 


It is then necessary, in series (3), to consider only the first terms of the sums in the determination of x, 
and xs. Conditions (5) and (7) take the form: 


A | W2 (iw) | sin [— wt + arg W; (iw)| = 0, (14) 


40 
= | Wa (je) | cos [— wt + arg Ws (iw)] <0, (15) 


and Conditions (6) and (8), respectively, 


| W; (ie) | sin [wt + arg W, (iw)] = 0, (16) 


(ie) | cos [ar + arg (ie)] > 0. (17) 


It follows from (14) and (15) that 
—t + arg W; (iw) (18) 
It follows from (16) and (17) that 
wt + arg W; (iw) = 0. (19) 
By summing (18) and (19) we find that 
arg W, (iw) + arg W2 (iw) = x. (20) 
Equation (20) can be used for determining the frequency wy by auto-oscillation, However, we find it 
more convenient to replace it by the equivalent system of two equations: 


sin [arg (iw) + arg (iw)] = 0, 
cos [arg W, (iw) + arg W2(iw)| <0 


U; (a) V2 (@) + U2(@) Vi (@) = 0, 
U; (@)|U2 (@) — Vi (@) V2 <0. (21) 


We will now seek the frequency of auto-oscillation by the harmonic balance method. As is known [2], 
for the determination of w»,in this case,we have the equation 
ImI (iw) = 0 (22) 
and the condition 
Rel (iw) < 0. (22") 
In our case, the reciprocal of the frequency characteristic is defined by the expression 


I = 


‘ics 


where 
W (iw) = W, (iw) Wa (iw) = U (@) + iV (@), 


where we have taken U(w) = Uy(w)U,(w) — Vy(w) (w) and V(w) = Uy(w)V,(w) + U,(w) Vy (w). 


Consequently, 
U (@) . 


Then, (22) and (22") take the form: 


V (w) = U; (@) Vz (@) + U2 (@) Vi (@) = 0, 
U (a) = U; (@) U2 (@) — (@) Ve <0. (23) 


We now determine the amplitudes A, and Ag of the periodic functions x(t) and x,(t) by the method of 
harmonic balance. We find A, from the conditions of amplitude balance 


I(iw) = -R(A, w), 
where 
R(A,w) = Ry (Aq) Rp (Ag). 


To determine A, we set w = Wy and eliminate A, by means of the equation describing the passage of the 
sinusoidal signal of frequency w, and amplitude A, through the portion of the circuit made of of “first relay- 
first filter": 


Ay =|W, (i@e) | Ri (A) Aj. (24) 


In our case, the equivalent transfer function of the first relay elements is Ry(Ay) = 4 mA, and Eq. (24) 
can be used for the determination of A, 


Ag = | Ws (iene) |= - (25) 
We determine A, analogously: 
4 
A; = | Ws (ie) | =. (26) 


A comparison of (23) with (21) convinces us that the equations for determining the auto-oscillation fre- 
quency obtained by the exact method (if the additional assumption is made that each of the linear portions of 
the system is an ideal filter) coincide with the equations obtained by the harmonic balance method. 


A comparison of the amplitudes A, and Az, determined in accordance with (25) and (26), with the ampli- 
tudes of the first harmonics of x, and x, which can be obtained from (3) also convinces us that, in the case we 
have been considering, they also coincide. 


From what has been presented we may also draw the converse inference: if at least one of the linear por~ 
tions passes at least one of the higher harmonics produced by the relays, then the exact and approximate equa~ 
tions will not coincide, 


We note in conclusion that the conclusions drawn may easily be extended to single-loop systems contain- 
ing an arbitrary number of symmetric relays which are separated by harmonic filters. 


Received May 17, 1958 
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CARBORUNDUM THERMOCOMPENSAT ORS 


R. I. Bresker and N. I. Voronin 


(Leningrad) 


High-ohmic and low-ohmic thermocompensating electrical impedances have 
been made of carborundum which allow heating up to 300°C and power dissipation 


up to 24 watts, 


There have been developed, in the Leningrad Institute for Refractory Materials, low-ohmic (from 5 to 
20 ohms) and high-ohmic (from 400 to 4000 ohms) thermocompensating impedances which, in the temperature 
range from —60 to +300°C, possess negative temperature coefficients equal, for 20°C, to 0.4 + 0.6% per °C, 


These impedances can dissipate power of up to 25 watts, 


R, ohm 
50 


100 200 


Fig. 1. Impedance variation as a function of 
ambient temperature, 


ohm 


Fig. 3, Impedance variation as a function of 
work-piece exposure time at 300°, Curve 1 is 

for an initial impedance of 1000 ohms, and 

Curve 2 is for an initial impedance of 1800 ohms, 


0-00 100 0 


Fig. 2. Temperature coefficient variation as a 
function of ambient temperature. 


Zamp 


Fig. 4. Volt-ampere characteristics of low- 
ohmic work-pieces, Curve 1 is for a piece with 
10 ohm impedance, and curve 2 is for a 25-ohm 
piece, 
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These thermocompensators are prepared from 
2 black carborundum, elemental silicon, titanium di- 
odixe and kaolin. The carborundum, titanium dioxide 
and silicon are used in invariable proportions — 40:30;10, 
The kaolin content varies from 10% to 40%, The 
moistened mixture, after stirring, is loaded into a hy- 
draulic broaching press, from which it is extruded 
through a spout 15 to 17 mm in diameter. The extruded 
material is cut into lengths of 65 to 70 mm and, after 


of eS 


Fig. 5. Volt-ampere characteristics of high- being dried, is baked in a coke charge in furnaces with 
ohmic work pieces, Curve 1 is for a 500-ohm oil heating with a temperature of 1200 to 1300°C for a 
piece, and curve 2 is for a 4000-ohm piece. period of 3 hours, Contact strata are deposited on the 


ends of the roasted pieces by spraying on aluminum 
leads in a Schoope apparatus, 


Figures 1 and 2 show the dependence of the electrical impedance R and the temperature coefficient of 
impedance @ on the ambient temperature, The measurements were made in the temperature range of —60°C 
to +300°C, In this interval the variation of impedance with temperature was exponential, The thermocompen- 
sators of carborundum changed their impedance in the first hour of operation, It is necessary to maintain the 
work-pieces at a temperature of 300° for some 50 to 60 hours in order for them to attain a state at which varia- 
tion of the electrical characteristics of the compensators ceases, Fig. 3 shows the curves of impedance varia- 
tions as functions of time of exposure to 300° in a furnace, 


The behavior of the thermocompensators was studied in a medium with 98% relative humidity, Exposure 
to such a medium for three to five days resulted in an increase of impedance by 10 to 30%, To protect the work- 
pieces from the action of humidity, they were sheathed in an enamel consisting of 70% lead oxide, 25% boric 
anhydride, and 5% pyrolusite. The enamel was sufficiently heat-resistant up to temperatures of 300-350°, The 
reliability of the enameled work-pieces was increased by a factor of two or three. 


Use of thermocompensators will be accompanied by significant heating of them, Currents up to 2 amps 
may flow through the low-ohmic impedances, It is very important to know to what temperature the compensator 
surface will heat when it is loaded with a power of 25 watts. Work-pieces of 10 and 16 mm diameter, both low- 
ohmic and high-ohmic, were investigated. The dissipating surface in the first case was 16 cm and, in the 
second case, 25 cm*, With 25 watt loads, currents up to 0.25 amps flowed through the high-ohmic pieces and 
currents of 2 amps flowed through the low-ohmic ones. In the steady state, surface temperatures for both types 
did not exceed 150°C, The compensators retained their stable electrical characteristics in use. Periodic or pro- 
longed loading with power up to 25 watts did not lead to any noticeable variation in the electrical characteristics 
of the previously aged work-pieces, 


Figure 4 gives the volt-ampere characteristics for low-ohmic work-pieces, while Fig. 5 provides them for 
high- ohmic pieces, 


The thermocompensators developed can be used successfully in those cases where, for implementing 
temperature compensation, impedances are required which are capable of withstanding large current loads and 
which possess impedance temperature coefficients equal in absolute magnitude to the impedance temperature 
coefficient of copper. 


Received May 4, 1958 
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CHRONICLE 
SEMINAR ON THE TECHNICAL APPLICATIONS OF MATHEMATICAL LOGIC 
(1957-1958)* 


During 1957 and 1958 the seminar on the technical applications of mathematical logic continued its 
meetings, devoting them almost entirely to the discussion of new original work by the seminar participants. 


During this time there were 16 sessions, of which 14 were devoted to original work and one to an instructive 
presentation of the translations of Soviet literature on the seminar subject matter in the USA. The work of the 
seminar participants was principally devoted to methods of analysis and synthesis of contact, and contactless, 
relay circuits, Foreign scientists took part in the work of the seminar, 


On September 20, 1957, in his paper, "On the fundamental concepts of relay-contact circuit theory,” 
devoted to general methodological questions in relay-contact circuit theory, G. N. Povarov expressed the opinion 
that in existing relay-contact circuit theory and, in particular, in the terminology recommended in 1953 by the 
KTT AN SSSR, a two-fold ambivalence of concepts can be observed: the structural point of view is mingled 
with the functional, and concepts from electrical engineering are mingled with information theoretical ones. 
According to the speaker's point of view, the functioning of a relay-contact circuit can be reduced to the proces- 
ses of changing the circuit's state — to its own type of absolutely inflexible nonhomogeneous Markoff chains. 

The concepts related to the changes in state of the circuit as a whole have, not a structural, but a purely func- 
tional character, and are applicable also to contactless relay circuits and to other switching systems, Such, for 
example, are the questions in classifying sequential relay-contact circuits or questions of determining the mini- 
mum number of internal circuit states necessary for storing, in sequential circuits, the necessary information as 
to previous states, Even questions of coding the circuit's state by combinations of states of its elements, and 
questions of interconnection of elements are amenable to treatment with structural concepts. 


Further, in the speaker's opinion, relay-contact circuits are essentially devices for the processing (trans- 
formation) of information, for which engineering questions are only means, not ends, While admitting the pos- 
sibility of treating relay-contact circuit theory as a branch of nonlinear electrical engineering, the speaker still 
asserted that a more general relay system theory can only be a portion of the theory of discrete transformers of 
information. In particular, pure'y functional questions in the theory of relay-contact circuits must be treated 
just in the domain of the theory of discrete information transformers, The speaker also made several remarks 
anent existing terminology. 


The paper gave rise to lively discussion, Present for the paper was the American scientist Koni, then in 
the Soviet Union as a tourist. 


V.N. Grebenshchimov (post-graduate physics student at the Moscow State University), in the paper 
"Methods of synthesizing multiterminal contact circuits of the general type," (October 25 and November 1, 1957), 
presented an algorithm for circuit synthesis using a method he had previously presented. The method consists of 
a successive removal of redundant contacts from the sought multiterminal network and in the determination of 
the contact admittance of the remaining portions of the circuit by solving certain systems of Boolean equations 
(DAN SSSR, 1958, 119, 2). The algorithm consists in writing, in some general form and using a special symbol- 
ism, all the solutions to each system of equations,and in then applying three types of selection rules, Use of the 


* Proceedings of the seminar prior to 1957 are reported in Automation and Remote Control (USSR) 18, 10 (1957). 
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first type of rule eliminates those superfluous solutions which, by different paths, would lead to identical circuits, 
The application of the second type of rule leads to the elimination of those solutions which would lead to cir- 
cuits which could be transformed into each other by renumbering the variables in the functions or by reposition- 
ing the contacts in the circuit without disrupting the contact admittance realized by the circuit. The use of 

the third type of rule leads to the elimination of all the remaining solutions but one. Changing the third type 
of selection rule allows one to obtain circuits which satisfy various requirements. 


Methods of analyzing multiterminal contact circuits were treated in the papers of V. N. Grebenshchikov, 
"A method of analyzing multiterminal contact circuits, and its mechanization,” (Feb. 6, 1958), "A method of 
analyzing multiterminal contact circuits with rectifying elements" (March 6, 1958), and "A method of analyz- 
ing multiterminal circuits consisting of contact multipoles" (March 20, 1958). 


The first of these papers presented an analytical method, suggested by the author, for multiterminal con- 
tact circuits, the method consisting of successive determinations of the contact admittances of chains of sub- 
circuits of the given circuit, each of which subcircuits is formed from the previous chain either by the operation 
of joining one contact to two terminals of the circuit or by the operation of introducing one free terminal. To 
determine the contact ad.nittance of each subcircuit, one uses a system of Boolean equations which relates the 
admittance of the new subcircuit with the contact admittance of the previous chain. Also presented was a 
method of mechanizing the foregoing method, which consisted of coding the Boolean functions in the binary 
number system (one number to each term of a Boolean function) and a reduction of the algebraic operations for 
transforming contact admittances to a set of machine operations on numbers, 


In the second paper, the method just presented was extended to circuits with rectifying elements by the 
use of functions of direct and reverse admittance for each terminal pair and by adding the operation of joining 
a rectifier to.two terminals of the circuit. 


In the third paper this method was extended to circuits of contact multipoles with given contact admit~- 
tances between the terminals of these multipoles, For this there was introduced the operation of connecting to 
two terminals of the circuit an element which posses a constant contact admittance of 1 in both directions, 


M. L, Tsetlin (post-graduate student in physics at the Moscow State University), in the paper "A matrix 
method for synthesizing relay-contact and electronic circuits," (December 19 and 26, 1957), presented a method 
suggested by himself for the synthesis of primitive and nonprimitive relay-contact and electronic relay circuits 
(DAN SSSR 117, 6 (1957) and 118, 3 (1958); Elektrosvyaz' 4 (1958)]. 


The author calls a relay circuit primitive if the set X; of quantities xf, . . . , xf *S which determine the 
states of the circuit's input busses completely determines the states of the circuit's output busses, i.e., if the 
following relationships hold: 


fi afi 1,2,..., p48). 


The author calls a circuit nonprimitive if it is obtained from a primitive circuit by the introduction of a 
feedback loop with unit time delay: 


= (i= 4, 2,....5), 


i.e., by the identification of the states of the input busses a. . at time t + 1 with the states of the output busses 
at time t. 


By use of a matrix method certain general properties of such circuits were established, and formulas were 
deduced which permit the structure of the circuit to be synthesized to be determined directly from the matrices, 


In the paper there was given a detailed presentation of a method for synthesizing nonprimitive circuits 
constructed of electromagnetic relays and triggers. The method is also easily generalized to circuits constructed 
from other elements. 


In the paper of M, L, Tsetlin and L, M, Shekhtman (NIISChETMASh Gosplans SSSR), "Ferrite~transistor 
circuits and methods of synthesizing them" (April 3 and 10, 1958), circuits were described which contained 
ferrite cores and transistors, principles were presented for the design of universal elementary cells which realize 
the function xy, and a methodology was presented for the synthesis of circuits constructed from these cells. 


The methodology for circuit synthesis consists of the introduction of a special symbolism for the element- 
ary cells and in the rules for the superposition of the elementary functions realized by the cells, It was shown 
that an arbitrary function can be realized in not more than two cycles (stages) of a device, and the necessary 
formulas were obtained. A technique for minimizing transformations was described. By the use of special 
formulas a method was presented for constructing nonprimitive circuits of ferrite-transistor elements directly 
from the matrices of these circuits, It was noted that the synthetic method presented by the authors was used 
in the design of digital computers. 


In the paper of V. 1, Shestakov, "A card method for synthesizing sequential relay systems,” (September 27, 
1957), the author gave his method for implementing the vector algebraic method of synthesizing two-position 
relay systems directly from special punched cards [Automation and Remote Control (USSR) 19, 6 (1958))? 


The vector-algebraic method of synthesizing two-position relay systems, presented in previously published 
works of the author [Automation and Remote Control 15, 4 (1954); Doklady AN SSSR 98, 4 (1954); and 99, 6 
(1954)) always permits the Boolean functions describing the conditions of operation of each relay of the system 
to be synthesized to be obtained from given sequences of operation of all the relays of this system, plus given 
sequences of external stimuli acting on this system, However, with a large number of relays in the system to 
be synthesized and a large number of independent stimuli acting on this system, the computation of the Boolean 
functions describing the operation of the system's relays requires a large expenditure of work and time, This 
gives rise to the requirement for alleviating and accelerating this computation by mechanizing at least part of 
it. The simplest method of effecting this mechanization is the punched-card method presented. 


V. I. Shestakov, in his papers, "Simulating the operations of D, A, Bochvar's calculus by three-position 
switches," (November 22, 1957), and "On one arithmetic interpretation of D, A. Bochvar's three-valued calculus,” 
(November 26, 1957), considered questions of simulating the operations of the three-valued propositional cal- 
culus by means of three-position relay-switch circuits, 


In the first of these papers the three-valued calculus of D. A. Bochvar was extended to a functionally 
complete propositional calculus and a method of simulating the operations of this extended calculus by means 
of three-position switches was considered, In this method, as the simulating parameter and the coefficient of 


voltage transmission, k = E,/E,, was chosen, where E, and E, are the voltages at the switch's input and output, 
respectively. The correspondence between the values of the three-valued proposition P and the coefficient k 
was established by the equation k = sign P, where the function sign P, “value of the proposition P,” introduced 
by the author, is defined by the following conditions; 


1 if P is true, 
sign P= 4 0 if P is vacuous 
—1 if P is false 


This method of simulating was treated in detail in the paper, “Simulating the operations of propositional 
calculus by means of relay-contact circuits," which was included in the collection of papers on logic published 
by the Philosophy Institute of the AN SSSR, 


The second paper of V. 1, Shestakoy was devoted to the author's new arithmetic interpretation of the func- 
tionally complete propositional calculus, The previous interpretation was based on the equation given above, 
k = sign P, The new interpretation is based on the equation k = sign I holding between the current I flowing 
through the relay winding switched to the output of the three-position switch, and the coefficient k, By eliminat- 
ing the coefficient k from these two equations, we obtain the equation sign P = sign I, which lies at the base of 
the new arithmetic interpretation of the three-valued propositional calculus, This equation permits each opera- 
tion of the functionally complete three-valued propositional calculus to be represented by a corresponding 
arithmetic operation, Each of these latter operations is realized by some relay~switch circuit. 


In the paper of Yu. A. Shreider (Moscow), "On the analysis and synthesis of certain electronic circuits,” 
(December 12, 1957), there was described the work carried out jointly by the author and B, I, Rameey, In this 
work a somewhat modified mathematical apparatus of the propositional calculus was used in the analysis and 
synthesis of networks constructed of diodes. Rules were formulated for finding the total number of diodes con- 
tained in a given circuit from the structural formulas of the circuit. 


*See C, B, translation, 
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Taking part in the seminar session of October 12, 1957 were Rumanian scientists led by Academician 
Gr. K. Moisil, They presented two papers on work on relay circuit theory in the Rumanian National Republic, 


In the paper of P.Constantinescu, “On circuits with the structural formula x,;@ xg®... @Xn," there were 
presented results which were contained in the paper previously published by the author, “Asupra reducerii 


numarului de contacte prin introducerea circuitelor in punte ," (Analele Universitatii "C. 1, Parhon," Bucuresti, 
No, 11, 1956, Ser, Stiin, Nat). 


In the paper of G, Ioanin, "On the synthesis of circuits with relays and stepping switches,” there were given 
methods for synthesizing circuits constructed of stepping switches and relays, these methods having been de- 


veloped by the author in cooperation with Gr. K, Moisil, As noted by the speaker, he tended to construct algo- 
rithmic methods which would be amenable to subsequent mechanization, 


The paper of M, V. Surgucheva (Physical Institute, AN LatvSSR), "A survey of digital computers from the 
point of view of the operations executed by them,” (May 15, 1958), contained a general classification of the 


types of operations executed in electronic digital computers, and estimates of these operations as applied to the 
practical construction and use of one of these machines. 


At the same seminar session (May 15, 1958), G. N. Povarov reported on the translations in the USA of the 
works of Soviet authors on relay circuit theory. The speaker gave a list of such translations made within the 
last few years by the translating firm of Morris D, Friedman, Inc. (Needham Heights, Mass., USA). 


V. N. Grebenshchikov 


CONFERENCE ON THE THEORY AND APPLICATION 
OF DISCRETE AUTOMATIC SYSTEMS 


From September 22 to 26, 1958 in the Moscow Dome, there was held a conference on the theory and ap- 
plication of discrete automatic systems, organized by the National Committee of the Soviet Union on Automatic 
Control in conjunction with the Institute for Automation and Remote Control of the AN SSSR. More than 600 
scientific workers and engineers from 265 organizations in 24 cities took part in the conference. 


The plenary session of the conference was opened by V. A. Trapeznikov, chairman of the National Com- 
mittee of the Soviet Union on Automatic Control, director of the Institute for Automation and Remote Control 
of the AN SSSR, and corresponding member of the AN SSSR. 


In his introductory remarks,V. A. Trapeznikoy stated that the basic task of the conference was a broad dis- 


cussion of the state of the theory and practical application of discrete automatic systems, and the prospects for 
their further development, 


One of the principal directions in the development of technological processes in our country is the broad 


development of automata, Automata cannot be successfully developed without being based on scientific achieve- 
ments, 


Recently, the science of automatic control has been concentrated in a new direction, related to the de- 
velopment of the theory and practice of discrete automatic systems. The use of such systems allows flexible 
control and adjustment of technological processes to be carried out, and allows great economies to be effected 
in many cases. V. A. Trapeznikov further stated that the present conference should be considered as a pre- 
liminary to the International Congress on Automatic Regulation and Control to be held in Moscow in 1960. 


The plenary session continued with a paper by Ya. Z. Tsypkin (Moscow) on the subject, "Discrete auto- 
matic systems and prospects for their development.” He noted that contemporary automation is characterized 
by complex problems, arising not only from the complexity of the objects of regulation, but from the increased 
requirements on intensification, economy, and effectiveness of the regulated processes, 


The solution of these problems led to the creation of a new class of automatic systems, possessing better 
qualities than the ordinary automatic control systems. In this class fall the so-called self-adjusting systems. 


These systems carry out analyses of the previous state of the process, observe or simulate the present state, 


| | 


ate, 


and a prediction of the future state. Execution of these operations requires a large amount of computation, 
both arithmetical and logical, The advent of digital computers had an essential effect on the development of 
self-adjusting systems, Instead of the ordinary setting of the given and desired values of the regulated quanti- 
ties, such systems carry out searches for optimal states (minimum control time, expense, power, maximum 
productivity, efficiency, etc.) 


Ya. Z. Tsypkin dealt in detail with the classification of discrete automatic systems and their character- 
istic idiosyncrasies. Using as the basis for his classificatory scheme the method of transforming and transmit- 
ting information (the method of quantizing physical quantities), the speaker distinguished three basic classes of 
discrete automatic systems: relay, pulsed, and pulse-relay or digital automatic systems, In relay systems, the 


physical quantities are characteristically quantized by level,while the quantization is by time in pulsed (sampled- 


data) systems, In systems using digital computing devices, quantification by level and by time occurs simul- 
taneously. 


Of particular interest in Ya. Z, Tsypkin's paper was the theoretical discussion of the principles of quanti- 
fication and the effects engendered by quantification, such as, for example, the observation of high-frequency 
periodic signals by means of narrow-band apparatus, stabilization, and improvement of dynamic properties of 
systems with lags, etc, Considered also were the theoretical bases of reproducing a continuous signal from its 
values at discrete moments of time, as well as the limitations inherent in time quantification. Quantification 
by level is a nonlinear process. However, if, instead of signals, one considers their distribution densities, then 
many of the results applicable to the operation of time quantification become applicable to this case as well, 


After a brief survey of the results obtained in discrete system theory [methods of analyzing and synthesiz- 
ing relay and pulsed (sampled-data) systems, approximate methods for investigating periodic regimens in ex- 
tremal systems], the speaker formulated a number of problems whose solutions are important for the further de~ 
velopment of the theory, Among these problems are; the analysis of self-adjusting systems described by a set 
of difference or differential equations plus equations describing logical operations; the analysis and comparison 
of various methods of searching; investigation of methods for coping with noise in self-adjusting systems, Solu~- 
tion of these problems will require the use of results from information theory, the theory of games and dynamic 
programming theory, decisive for the various optimization problems, 


The further work of the conference was carried out in the following sections: 

1, Section on the theory of sampled-data and digital control systems, headed by G, S. Pospelov, 
2. Section on self-adjusting (extremal) systems, directed by A. A. Fel'dbaum, 

3. Section on optimal system theory, headed by A. Ya. Lerner. 

4, Section on sampled-data and digital elements and devices, headed by V, A. Kotel*nikov. 

5. Section on finite automata theory, headed by M, A. Gavrilov, 


The section on the theory of sampled-data and digital control systems contained 11 papers, The speakers 
covered a broad range of questions. 


The paper of G, P, Tartakovskii (Moscow) contained a generalization of the basic theory of sampled-data 
systems with constant parameters to systems with variable parameters. For the various cases of parameter varia- 
tion in the linear and pulsed portions, the author introduced the concept of the transmission function depending 
on time, by means of which the reaction of the system to arbitrary regular input stimuli is found, Also general- 
ized were the stability criteria of systems with constant parameters to systems with variable parameters, Using 
a pulsed radio range finder as an example, the author showed the methodology of investigating stationary and 
nonstationary discrete processes, 


The paper of V. P. Perov (Leningrad) was devoted to the statistical synthesis of sampled-data systems with 
variable pulse repetition periods, the variation occurring either in accordance with a previously given rule, or 
randomly, The speaker considered the possibility of constructing an optimal system where the criterion would 
be minimum dispersion of the system error. A generalization of the method of synthesizing optimal systems to 
the case where the variation of pulse repetition period had an arbitrary character allowed the author to draw cer- 
tain conclusions as to the possibility of constructing self-adjusting optimal systems, 
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In the paper of Fan’ Chun-Wui (Moscow), a method was presented for the analysis of sampled-data systems 
in which the pulse sequence is set up by a pulse element with complicated modulation: by frequency and by 
amplitude with weighted coefficients. Based on this method were the considerations of applications of compli- 
cated modulation for solving problems of interpolation, for implementing systems which are optimally fast, for 
increasing noise imperviousness, the quality of signal reproduction, and the degree of system stability. 


As one particular case of a sampled-data system with complicated modulation there was analyzed a 
system containing several pulse elements with different repetition periods. 


In his paper, F, M, Kilin (Leningrad) considered exact and approximate methods for analyzing and syn- 
thesizing sampled-data systems with variable parameters, in which the variation is stepwise. The author adduced 
results of an investigation of the statistical dynamics of such sampled-data systems with gating of random signals 
by a system of coincidence amplifiers, and he also presented criteria for optimal systems based on the conditions 
of noise passage for the case when the noise is produced in the intervals between control pulses, 


In his second paper, Ya, Z. Tsypkin gave a general formulation of the problem of synthesizing the struc- 
tures of pulsed devices designed to compensate for the effect of object delays on the dynamic characteristics of 
the servo system, He obtained the conditions for the possibility of complete and of partial lag compensation, 
and considered the example of synthesizing a compensated optimal system. 


A. A, Krasovskii (Moscow) devoted his paper to the synthesis of self-adjusting systems with digital correct- 
ing devices, The author showed how to choose the transfer function of the discrete correcting device of a linear 
control system which would guarantee optimal processes in the closed system for given characteristics of the ob- 
ject and of the executive portion of the regulator. He then presented a method for adding subsidiary correcting 
links to the self-adjusting system which would guarantee the automatic adjustment of the discrete correcting 
device when there were changes in the characteristics of the object or the executive organ. The method of 
synthesis presupposes that the probability characteristics of the noise are such that the correlation time is suf- 
ficiently small, With this, the optimal distribution of the weights of the closed-system coefficients is a linear 
one, 


In his paper, E. A. Krogius (Leningrad) gave a frequency method for designing linear sampled-data control 


systems with random noise. The author obtained simple relationships between the parameters of the correcting 
device and those of the basic system parts, The method is based on the use of logarithmic frequency character- 
istics, In the paper were given conditions under which a sampled-data control system can be considered as a 
continuous one. 


In his paper, E. M. Esipovich (Moscow) introduced the concept of the so-called conditional transfer func- 
tion which, in contradistinction to the ordinary transfer function, depends on the form of the input stimulus, The 
author showed that the properties of the conditional transfer functions are such that the theory of sampled-data 
control can be extended to cover them, 


Certain features peculiar to conditional transfer functions were illustrated in the example of a discrete 
system with a linear extrapolating device operating under noisy conditions, 


I, V. Pyshkin (Moscow) devoted his paper to the analysis of auto-oscillations in systems with pulse-width 
modulation, In the paper were given the necessary conditions for the existence of symmetric auto-oscillations, 
expressed in terms of the frequency characteristics of the continuous portions of the system, and a method of in- 
vesting the stability of these auto-oscillations was given. For the simplest forms of auto~oscillations, the 
speaker gave nomograms which permit one to determine, directly from the form of the frequency characteristics 
of the linear portion, whether the conditions for existence of auto-oscillations hold. 


I. I, Perel’man (Moscow) gave, in his paper, the design principles for automatic control systems called, 
by the author, systems with self-established programs, Such systems guarantee the required responses to the sig~ 
nal components whose durations are either greater, or smaller, than the net object lag, Systems with self- 
established programs are designed to work with stimuli which repeat in time (are cyclical), As an example, 
the author considered a system for stabilizing the thickness of hot-rolled steel sheet in a continuous rolling mill. 
The paper contained an investigation of the stability of a system with a self-established program, and provided 
data on the choice of the optimal coefficients when the input to the system contained a cyclical stimulus accom= 
panied by random noise. 
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In the paper of B, A. Sigov (Kiev) there was provided an investigation of position, velocity, and integral 
systems of digital program control, A comparative investigation of the different structures was given, using as 
an example the control system for copying machines, The author explained a new structural scheme for a 
system of digital program control (DPC) which is called “integral,” In an integral DPC the integral requiring 
translating is written, but the translating itself is done automatically. This leads to a definite shortening of 


program length, 


Thirteen papers were presented in the section on self-adjusting (extremal)systems, these papers being 
devoted to self-adjusting system theory, design principles for extremal systems, and also to descriptions of ap- 
paratus tested in the laboratory or in prototype use, 


The paper of Yu. V. Dolgolenko (Leningrad) was devoted to the derivation of precise relationships between 
the parameters of relay extremal systems and the frequency and amplitude of auto-oscillations, The author 
compared the exact method with the harmonic balance method, previously used for the same purpose by I, S, 
Morosanov, 


As a result of this comparison, Yu. V. Dolgolenko drew the conclusion that it is possible to employ the 
harmonic balance method in synthesizing systems of the class considered, since the errors introduced by this 
method for high frequencies are small, 


In his paper, I. S. Morosanov (Moscow) considered an approximate frequency method for calculating the 
steady-state periodic regimens in extremal control systems with independent scans, Expressions were given for 
the equivalent complex gain for various relay and relay-pulse extremal systems, The speaker showed that the 
necessary conditions for the existence of auto-oscillations may be determined graphically from the equivalent 
frequency characteristics of the object. 


Yu. I, Ostrovskii (M.scow) devoted his paper to a method of designing extremal systems when noise is 
taken into account, The m, ‘thod was illustrated by the example of the design of a pneumatic extremum-regu~ 
lator, The essential point of the method is the following, Initially, the system is considered without noise at 
its output and with a given sinusoidal stimulus at the input, In such a system, the harmonic balance method al- 
lows one to find the errors accompanying this input excitation, Then, separately, one finds the probability of 
an incorrect step, or the frequency of incorrect steps, at the object output when noise is present, The system 
parameters are then so chosen that, for a given error, the attendant frequency of incorrect steps will be a mini- 
mum, The author provided experimental verification of the results of his theoretical investigation. 


The paper of V, V. Kazakevich (Moscow) considered questions of stability and quality of extremal con- 
trol systems, and also the principles of an extremal controller with a commutator implementing a forced reversal, 
and methods of introducing dummy extremum indicators, 


V. V. Kazakevich questioned the feasibility of extremal controllers which are developed without external 
disturbances being taken into consideration. If the speed of variation of the external excitations is commensu~- 
rate with the speed of processing of the stimuli to be regulated, the system may turn out to be unstable and, in 
order to avoid incorrect system operation and its passage from an extremum to a limiting position, it would be 
necessary to carry out periodically a forced reversal, The speaker suggested that, in those cases where the char- 
acteristic either has a bounded limb or no extremal indicator expressed at all, the characteristic be transformed, 
by the introduction of additional connections, so that the extremal indicator would become perceptible. Among 
the methods of improving the quality of pulsed-data systems there is the use of a variable processing step, by 
means of “steregushchikh" contacts, which do not permit extrema to occur abruptly. 


The paper of V. M. Kuntsevich (Kiev) contained a frequency method for synthesizing extremal systems of 
the stepping type, the method being convenient for engineering computations, Considered in the paper was the 
block schematic of a system, developed at the Institute for Electrical Technology of the AN SSSR, and data were 
Biven on the introduction of this system into industrial use for the maintenance at the maximum level of the 
efficiency coefficient of the steam boiler of a generating station. Work on the installed regulator is now at the 
stage of experimental testing. The first trials of the regulator working with an object, carried out jointly by the 
author and Yu, M. Bulavitskii, gave encouraging results, 


L. N. Fitsner (Moscow), in his paper, described two forms of automatic optimizers: the simplest relay type 
and the proportional action type with small additional perturbations introduced into the system during the scanning 
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process. The author showed that, in an inertial system operating with a relay optimizer, three fundamental 
forms of motion are possible, and he found the relationships between the object and optimizer parameters for 
which the scan process is stable and the amplitude of hunting is minimum, 


The paper of E, K. Krug (Moscow) treated of extremal systems for the regulation of objects with respect to 
two and more parameters. The method for finding extrema is based on the determination of partial derivatives 
at discrete moments of time, followed by an incrementing of the input quantities proportionally to these deriva- 
tives. The paper contained recommendations for the choice of amplitude, phase and frequency of the harmonic 
components of the input quantities, coeffi sients of proportionality between the increments of the input quantities 
and the values of the derivatives obtained as functions of the static characteristics and inertia of the controlled 
object, and also contained a short description of a regulator which provides an extremum in accordance with the 
method cited. 


In the presentation of his second paper, V. V. Kazakevich gave the results of an investigation of the non- 
linear processes in extremal regulators, The Galerkin method underlay the investigation into the steady-state 
and transient processes in the nonlinear extremal control system. In getting the form of the periodic solutions 
(including those different from the harmonics) for given boundary conditions, the author obtained a complete 
description of the processes. 


The paper of A, Fel'dbaum (Moscow) was devoted to the synthesis of optimal, or near-optimal, extremal 
control systems, in the sense of minimum duration of the transient response and maximum accuracy of the ex- 
tremum determined , taking noise into account, The average scan time of such a system is composed of a num- 
ber of discrete steps necessary for reaching an extremum from given initial conditions, plus a number of trial 
repetitions necessary for choosing the correct direction for the next step, given the probability of an incorrect 
step. Investigated also were several variations of the step system and systems with extrapolation for normally 
and equally distributed noise. One result of the investigation was that the trial repetitions are either unnecessary 
or very little necessary, since increases of the average time due to incorrect operations were, in the cases con- 
sidered, always compensated for almost completely by the time not taken up by the trials which weren't made, 


In a communication he presented, Sun Tszian* (Moscow) stated that the statistical method of analyzing 
automatic scanning processes may play an essential role in the investigation of discrete automatic scanning sys- 
. tems, particularly for those systems which contain complex logical portions, The methods in use up to now, 
which are based on spectral theory, cannot give a complete probability representation of the convergence, and 
of other characteristics, of the scanning process. Further development of the statistical method of investigation 
which was discussed in A. A, Fel'dbaum's paper may lead in the direction of a solution to the following problems: 
to determine the noise stability of a scanning system as a whole; to determine the probability that a system, ex~ 
posed to noise of a definite type, will arrive at an extremum after a determinate number of steps; to determine 
the optimal scanning strategy whereby the probability of arriving at an extremum will be a maximum for a given 
type and level of noise, The author obtained satisfactory solutions for several of the problems just formulated, 


The paper of R, 1, Stakhoyskii (Moscow) contained a comparison of several scanning methods for automatic 
optimizers, The author considered the Gauss-Seidel method of alternate changes, the gradient method, and the 
method of steepest descent. The comparison was made from the point of view of minimum average scan time 
and increase of accuracy in the extrema determined, It was proven that, for a system which moves by steps, 
the best scanning algorithm should be taken as a combination of the methods of steepest descent and gradients, 


The paper of R, V. Kornilovy and N, T, Khristoforoy (Moscow) was devoted to the description of an elec- 
tronic extremal regulator which retains ("memorizes") the extrema, together with an analysis of its operation. 
Laboratory testing of the system showed that its operation was stable with objects having significant inertia, and 
also with intense external perturbations of a monotonic character present, The device may be easily coupled 
with industrial electronic regulators. 


M, G, Eskin (Moscow), in his paper, considered an extremal control system for turbine boring with discrete 
devices, The optimum regimen of the turbine boring is established by so choosing the axial force on the gouge 
that maximum boring speed is obtained. The absence of any direct means for measuring the face parameters 
leads to the necessity of indirect measurement of these parameters. To do this one must take into account the 
dynamic state of the column of boring tubing and the columnar friction against the well shaft (random noise). 
The relative slowness of the departure from extrema suggests a two-stage control system, An ordinary control 
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system maintains a given established load on the gouge. This "given" magnitude of load is established by an 
extremal regulator, The paper considered relay and proportional systems of extremal control with shutdown in 
extremal zones and self-starting with departures from extremal zones. In the conclusion, the paper contained 
a description of the block schematic of a circuit with discrete elements, as well as the results of an investiga- 
tion into devices using a special dynamic model of a turbine borer. 


In the communication of N. V. Grishko (Moscow) there was formulated a criterion for the quality of ex- 
tremal systems in given statistical conditions which are characterized by disturbances of the extremum states 
and by high-frequency and low-frequency noise. Results were given of an experimental investigation of 
sampled-data systems, optimal in the sense of the given criterion, with decaying extrema. 


The work of the section on optimal system theory was devoted to the consideration of questions of con- 
structing automatic control systems which are optimal in respect to their rapidity of action. 


The paper of A. G, Butkovskii and S, M, Domanitskii (Moscow) presented the results of work on the crea- 
tion of optimal control systems for objects with delays in the case when the aim of the system is to work off 
initial mismatches, The authors succeeded in constructing a control system with a very simple structure, The 
system was used for the automatic matching of the flying blade of a rolling mill with the motion of the rolled 
strip. As shown by investigation, the system's accuracy was completely satisfactory for practical purposes, 


The paper of G. A, Nadzhafova (Moscow) showed ways for determining the limiting rapidity of action of 
a control system for objects with limitations up to the second order, The results were obtained by plotting the 
limits of the isochronous domain. This rendered it possible to limit the domain of object states from which ob- 
jects could be transferred to given states with defined boundaries. 


O, G, Varshavshii (Moscow), in his paper, presented the results of the first steps in the work of construct- 
ing a general industrial controller for objects with delays, where the control process is close to an optimal one, 
An analytical investigation and a comparative study of the dynamic properties of ordinary industrial controllers 
and a controller of the suggested structure showed the essential superiority of the latter with respect to the speed 
with which the control process was accomplished, 


The communication of A, A. Pavlov (Moscow) was devoted to a new method of analyzing optimal systems, 
illustrated by the example of the synthesis of a control system for an object with third-order boundaries, The 
special feature of the author's approach to solving the problem is his use of plane movable switching surfaces 
for the formation of the optimal control laws, 


There were 15 papers in the section on sampled-data and digital elements and devices. Considered in 
the papers were questions of the automation of complex technological processes entailing a large amount of 
computation, and also questions of creating simplified special-purpose compyters, devices for discontinuous con- 
trol, and elements for computers for a broad class of automatic devices, 


In his paper, A. G. Moskalev (Moscow) dealt with the problem of the automatic control of the working 
regimen of an energy system both by frequency and by active power so as to provide the economically best 
working regimen for maintaining the frequency at a given level,and also for providing high reliability both in 
normal and in emergency regimens, The author spoke of two control systems which were implemented and 
tested in actual working conditions during 1956 and 1957 in the energy systems of the Kherson and Khar"kov 
Councils of National Economy. Automatic control of the working regimen was implemented on the majority of 
generating stations, with a nominal power between 80% and 85% of all the nominal power of the energy systems. 
The experiments provided good results. 


The paper of O, M, Kryzhanovskii and M, Yu, Beiderman (Kiev) was devoted to automatic control systems 
for mines’ multicable hoists which provide both accident-free and maximally productive machine operation, 
The paper gave particular attention to the development of indicators for deep mines and to measuring elements 
for the control systems, Both devices were implemented as sampled-data devices, and operated with high accu- 
racy and reliability. 


The paper of S, M, Khlytchiev and G, A. Aleksandrov (Moscow) dealt with questions concerning the auto- 
mation of radio receivers in a long-line short-wave radio communications system, The authors showed that it 
is economically most efficient to automate completely the receiving centers, The paper contained block schem- 
atics of automation devices based on the use of electronic controlling machines, 
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The paper of Yu. I, Kaznacheev, N, D, Kolesnikova and Yu. N. Prozorovskii ( Moscow) was devoted to 
velocity electrovacuum instruments, transforming continuous data to discrete, and vice versa, for broad-band 
message transmission systems with pulse~code modulation, The authors described coding and regenerating de- 
vices constructed on the basis of electron-beam tubes, 


B, L, Ermilov (Leningrad) gave a paper on digital analogs of analog computers, these digital computers 
being sufficiently compact, maintaining the high accuracy of digital computers while combining this with the 
speed of analog devices. The special features of these digital analogs were shown by the problem of controlling 
a milling machine. 


In his paper, A. K. Zavolokin (Moscow) presented a classification of transformers of continuous quantities 
into digital equivalents. On the basis of this classification, the author presented a method for obtaining digital 
equivalents which characterizes, with sufficient definiteness, such important considerations for applications as 
rapidity of action, noise-resistance, the cyclical nature of the operation, and the unambiguity of null-point 
definition, The paper treated in detail groups of transformers of serial counters, counter transformers and digit- 
by-digit coding transformers. 


The paper of M. A. Shnaidman (Moscow) was devoted to the investigation of a pulsed-data control system 
by means of an analog computer, The problem of creating means for simulating pulsed-data systems leads, 
essentially, to the construc:ion of pulsed-data blocks and their connecting with de analog devices, to be manu- 
factured by industry, which contain both linear and nonlinear elements, The paper considered circuits for blocks 
which could provide any form of pulse modulation, these blocks being based on standard operational amplifiers 
which allowed very high-quality characteristics to be obtained from them. 


In the paper of E, L, Itskovich (Leningrad) there was considered a discontinuous-type universal correcting 
device for implementing a given control of an aggregate's output parameters by varying the regiinen by which 
its operations are stabilized by means of stimuli to the regulator transducers of the stabilizing circuit. The 
author provided a method for designing correcting devices so as to insure stability and quantitative control of 
the aggregate parameters under consideration and for creating the conditions of autonomy for other of its param- 
eters. The author cited positive results from the application of the correcting devices developed in a circuit 
for regulating the process of roasting cement in a rotating kiln. 


The paper of V. A. Oleinikov (Leningrad) was devoted to the operation of an electromachine amplifier 
(EMA) with transverse poles in a sampled-data mode of operation, The author showed that the introduction of 
the sampled-data mode of operation, due to the vibrational linearization of the nonlinearities, significantly im- 
proved the free-running characteristics and the external characteristics of the EMA. The hysteresis loops of the 
characteristics cited were decreased by a factor between 2 and 3, and this improved the linearity of the character 
istics themselves, , 


The paper contained recommendations on the choice of frequency and pulse length for guaranteeing that 
the EMA would function within given limits of accuracy. 


The paper of B, M. Kagan, V. M. Dolkart, A. I. Voitelev and B, P, Brudnii (Moscow) described a complex 
universal discrete-continuous machine for the investigation of the processes in control systems controlled by 
digital machines. The setup consisted of a whole bank of machines. 


In the paper of Z.B,Vartanov (Moscow) there was considered the circuit for a discrete machine for con- 
trolling a blooming mill. The apparatus consists of a digital servo system of the position type in which a static 
program is introduced by means of a memory device, 


The paper of V. A. Brik (Moscow) investigated the possibilities of constructing simple computers to carry 
out programs for controlling machines for piece-working, where the profiles of the pieces consist of straight 
lines and circular arcs. The paper gave special attention to methods for most simply shaping circular arcs, The 
author showed that the operation of multiplication was not necessary for this sort of program. This significantly 
simplifies the computer, shortens computation time, and reduces computational error to a minimum, 


The paper of G. G. Men'shikov (Leningrad) was devoted to a planned digital computer for computing lineaf 
combinations of functions, where the methods used did not include multiplication, The machine for solving the 
problems formulated would not only lack the operation of multiplication, but would also shorten the total numbef 
of operations, The author provided results of experimental testing of this method. 
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The paper of Ya, E, Gel'fand (Leningrad) considered a device for discontinuous control, with dependent 
intervals, for cement-producing aggregates. 


The very slow flow of the processes and the large lags in the object, sometimes as much as tens of minutes, 
made it preferable to introduce a discontinuous mode of operation into the controller, The best qualitative con- 
trol is obtained by introducing a variable control interval as a function of the deviation of the controlled quantity. 


The paper contained a description of the developed controlling device, currently being produced in small quanti- 
ties. 


In his paper, 1, O, Oskolkov (Moscow) presented a method for the multipoint inspection of parameter de- 
viations from given levels by the use of electron-tube indicators for the simultaneous measurement of the mag- 
nitudes of the deviations of a large number of parameters, The author considered the possibilities of implement- 
ing two-position and three-position control, His paper was based on material obtained from testing of a develop- 
ed system. 


The paper of V. M. Eigenbrot (Moscow) was devoted to a discrete three-position regulator of twenty param- 
eters by one regulator where the objects have large time constants and are self-aligning, The author answered 
questions concerning the practical building of the circuit, the choice and disposition of the equipment, the 
method of adjusting the regulator, and he provided results of industrial testing. 


In the section of finite automata theory there were considered several questions dealing with the theory 
of automatic systems which consist of finite numbers of elements, each of which can have a finite number of 
states, 


M, A, Gavrilov (Moscow), in his paper, formulated the problem of constructing complex discrete remote 
control systems and showed that, by their functions, these systems fall in the class of multiprogram finite auto- 
mata with discrete feedback loops and with the selection of the program depending both on external stimuli and 
on the results of the operations of the system itself. The author also showed that further development of discrete 
remote control systems leads to self-organizing systems, 


The paper of P, P. Parkhomenko (Moscow) presented the underlying principles and the mechanisms of a 
special logical machine for the analysis of relay devices. The machine carries out structural analyses of planned 


or already constructed relay systems by the method of decomposing the circuits into constituent units, and de- 
rives from these decompositions the sequences of operations of the circuit elements, the results being taken off 
either manually or automatically, 


The paper of Yu. Ya. Bazilevskii (Moscow) considered finite determinate automata with constant structures 
which represented arbitrary occurrences of several arguments. It was shown by the author that from the calculus 
of temporal logical functions one may form a system of equations and may determine the properties of the 
system of automata states. The author considered the relationships between abstract automata and the families 
of subautomata derivable from them, gave a method for decomposing automata into nets with independent nuclei 
and _ presented several forms of equivalent transformations and a method for minimizing the matrix of states, 


The second part of the paper dealt with a n.ethod for solving generalized temporal logical equations in implicit 
form. 


In the paper of G, K, Berends and A, A, Talia (Moscow) on pneumatic element systems of finite automata, 
there was given a description of logical elements “not,” “and,” and “or,” of elements which realize the logical 
functions of "implication" and “equivalence” and also of delay and repetition elements, Based on these pneu- 
matic elements, devices of industrial importance were developed: regulators with automatic crossover, semi- 
Proportional regulators, and extremal regulators with several regulating organs. 


The communication of V, N. Roginskii was devoted to a description of a machine for the synthesis of relay 
device structure, This machine is completely automated and deduces the structure of a device from a switching 


table for it, given in graphic form, the machine simultaneously determinirg the required number of additional 
elements, 


The majority of the papers presented at the conference gave rise to lively discussions, While the theoreti- 
cal level of the papers presented was generally high, the degree of technological development was inadequate, 
Representatives of institutional branches spoke of their developments in the domains considered, and of the re~ 
sults of testing in use. 
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At the concluding session, the section leaders briefly addressed the entire conference, and also formulated 
the essential problems for further development in the theory and application of discrete automatic systems. 


I. S. Morosanoy and P, I. Chinaev 


INFORMATION 


FIRST IFAC CONGRESS 
Plans for the Scientific Program of the First Congress 


of the International Federation for Automatic Control (IFAC) in Moscow 


The plans were worked out and presented by the National Committee of the USSR for Automatic Control, 
and approved by the Executive Committee of IFAC at the meeting in Zurich in March, 1958, 


1. General Characteristics of the Congress 


The first international congress of IFAC was called by a resolution of the General Assembly of IFAC which 
was passed in Paris on September 12, 1957, Tentatively, the total number of congress participants was set at 
1500. 


The congress must include three types of scientific and engineering activities in the field of automation, 
Branch 1. Theory and methods of automatic regulation and control. 


Branch 2. Technological means for automation (switching and computing devices) and the problems re- 
lated to their creation on new scientific bases. 


Branch 3, Industrial applications of automatic regulation and control (including the use of computing 
devices), 


The watchwords of the congress should be “Theory — practical application, technological means ~ maxi~- 
mum reliability, automation of production — high effectiveness.” 


It is assumed that the congress will last 10 days and will include: a) the reading and discussion of scientific 
Papers; b) excursions to scientific institutions and also to various undertakings equipped with automatic devices; 
and, c) visits to the cultural centers of the Soviet Union, 


It is assumed that a special program will be set up for the excursions and visits to the Moscow cultural 
centers for participants and guests of the congress. 


2. Scientific Program of the Congress 


Because, as mentioned, there are three major forms of scientific and engineering activity in the field of 
automation, the scientific program of the congress will be given in three sections, The program contains an 
enumeration of the scientific problems and the objects of investigation, Both types of enumeration are presented 
with a broad brush and do not necessarily exhaust the themes of the congress. 


It is assumed that in the first section there will be considered general and particular theoretical problems 
in the investigation and design of automatic regulation and control systems, as well as methods of solving these 
problems, Any general or special mathematical posing of a problem and any (exact or approximate) method of 
solving it will be admissible. 


The range of investigations in the theoretical domain will include both abstract problems, which are posed 
with the aim of finding new methods of solving them, and problems with concrete physical content, which will 
specify the results obtained which are of scientific or practical value and their suitability for use in engineering 
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calculations, The objects of investigation may be open-looped or closed-looped systems, with one or several 
inputs and outputs, with constant parameters, deterministic or stochastic, as well as systems which contain 
humans or computers as links. 


In the second section, the objects of investigation are almost always elementary components of automa- 
tion systems: transducers, transforming devices, controllers and executive organs, For the entire set of these 
devices it is assumed that there will be scientific papers on the theory and design of transducers as measuring 
instruments, on the theory and design of devices which transform one physical quantity into another and, finally, 
the theory and design of controllers and executive mechanisms. 


The general problem of this theory and design reduces to the determination of exact data on these ele- 
ments, and also the determination of their static and dynamic characteristics, This same problem also pertains 
to any set of elements of technological means which form a single system for automatic inspection, regulation, 
and control, or of computing devices. Also relevant here is the most important problem of technological means, 
the problem of maximum reliability. The presentation of solutions of these problems can be illustrated by the 
results of experimental study of the technological means, and also by generalizations of their trials in use, 


The same section must include descriptions of the principles of concocting technological means, illumi- 
nated by trials of their design, preparation, and questions of economy. Here also should be presented new phy- 
sical principles on which may be based the development and design of new instruments, automation systems, 
and inspecting and computing devices. 


An important place in the third section is filled by theoretical and experimental methods for studying the 
static and dynamic properties of objects and technological processes, and also methods for the analysis of dis- 
turbing stimuli, 


All these data, together with the mathematical descriptions of the objects and technological processes to 
be automated, as well as with generalizations from the empirical testing of these objects, will be used in for- 
mulating the scientific bases of automatic production, On this basis there may be presented the design principles 
of automatic inspection and control systems, automated complexes, and also the bases of the automatic control 
of production flow, Also included in this section should be new regulation and control principles. 


The final problem is that of improving the conditions, and increasing the productivity, of work, as well 
as the economic efficiency of automation, The objects of investigation here are either the processes to be auto- 
mated, grouped by some inherent physical or other principle of their action, or the automated objects. 


In every section the problem of the historical investigation of the area of — as well as the prob- 
lem of creating a single international terminology, will be discussed. 


The subject matter of a paper can be any scientific or technological question given on our program, or 
close to it, Below we give a list of the scientific problems and objects relating to them in all three sections of 


the program.* 
Section 1, Theory, General and Particular Problems in Automatic Regulation and Control Theory. 
1) Approximations to, and simplifications of, design methods for controllers and for control systems. 


2) Methods for determining optimal adjustments of automatic controllers, 


3) Methods for finding the design characteristics of a system based on the use of experimental data, 
4) Analyses of concrete control systems and generalization of their results, 
5) Applied methods of synthesizing optimal systems, 


6) Nomograms, tables, graphs, slide rules, and other subsidiary aids for the analysis and synthesis of con- 
trol systeins. 


*For example, problem 1, “Approximations to, and simplifications of, design methods for controllers and for 
control systems,” is relevant to all the objects from 1 to 10 listed under "Objects of investigation" in the given 
Section, etc, 
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1) The stability problem, 
8) Periodic motion problems, 


9) The problem of estimating the quality of the transient response, as well as constructing it, The theory 
of sliding motion, 


10) The problems of structure synthesis and parameter selection, 

11) The problem of optimal regulation and control, 

12) The problem of programmed regulation and control. 

13) The problem of invariance and the theory of combined systems. 
14) The statistical theory of automatic regulation and control systems. 
15) Programroing theory. 

16) The theory of operation and logic diagrams, 

17) The theory of the transmission of remote control signals, 

18) Simulation and machine synthesis, 

19) The statics of controlled systems and the computation of their steady states. 
20) Terminology and history. 


Objects. of Investigation. 


1) Continuous systems (analog systems), 

2) Systems with discontinuous characteristics, 

3) Pulsed (sampled-data) systems. 

4) Regulation and control systems with computing devices. 

5) Self- adjusting (self-organizing) systems, 

6) Systems of extremal control. 

7) Servo systems. | 

8) Prediction systems, 

9) Systems with programmed control, | 
10) Systems with combined forms of operation, | 


Section 2, Technological Means, The Problems and Theory of the Technological Means of Automatic 
Control, 


1) Theory and computational methods for the static and dynamic characteristics of the elements of auto- 
matic regulation and control systems. 


2) Methods for the experimental investigation of the properties of elements and systems of automatic in- 
spection, regulation, and control, 


3) The analysis of disturbing stimuli. 
4) Trials in use, 


5) Classification, principles, and the theoretical bases of designs for elements and systems of automatic 
inspection, regulation, and control. 


6) The design of technological means, 
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1) The technology of technological means. 
8) The economics of technological means. 
9) New physical principles. 

10) Terminology and history. 


Objects of Investigation, 

1) Transducers and transformers. 

2) Controllers and their elements, 

3) Executive organs. 

4) Systems of automatic inspection, regulation, and control. 
5) Systems of remote regulation and control, 

6) Controlling machines, 


Section 3. The Automation of Production, Problems of Control Theory, 


1) Methods of investigating and designing the static and dynamic characteristics of the objects of regula- 
tion and control of technological processes. 


2) Methods for the experimental investigation of the properties of objects of regulation and control, and 
of technological processes. 


3) Analysis of disturbing stimuli. 
4) The classification and mathematical description of the objects of control and of technological processes, 
5) Trials in use. 


6) Principles and theoretical basis of designing automatic regulation and control systems. 
7) Automatic control of production flow. 


8) The development of automatic complexes (automatic factories, plants, etc), 
9) New principles of regulation and control, 
10) The productivity and economic effectiveness of automation. 
11) Terminology and history. 
Objects of Investigation, 
1) Electrical machines: generators, motors, transformers. 


2) Prime movers: turbines, motors, steam machines. 
3) Thermotechnological machines: furnaces, steam boilers, reactors. 
4) Chemical processes. 
5) Metallurgical processes: reducing processes, steel working, nonferrous metal production, alloys. 
6) Electrometallurgical processes: electrolysis, plating. 
1) Hydrometallurgical processes. 
8) Assembly-line material flow: rolling mills, paper-making machines, textile finishing processes. 
9) Forging and construction, 
10) Transport: railroads, cranes, rail-less wheeled transport. 


egula- 


11) Flying and sea-going apparatus: airplanes, ships, boats, rockets, torpedoes, 
12) Mechanical processes for material processing: mills, presses, hammers. 


13) Processes for handling freely-flowing material: concrete, plastic mixtures, ore concentration, 
neutralization and filling of ores. 


14) Mining and extracting operations: oil drilling, water, gas, oil and gas piping, water piping. 


3. Procedures for Preparation and Selection of Papers 


It is proposed to establish the following procedures for the selection and preparation of papers for the con- 
gress. 


1) In each country, the National Committee must make the selection of papers submitted for presentation 
at the congress. 


2) The papers must contain scientific (theoretical or experimental) or practical results of the work of one 
or a group of authors in any of the enumerated domains, or in a new domain related to those given in Section 2 


above, as well as subjects related to the development of methods of designing automatic systems, or experience 
in their industrial use. 


3) The National Committees will send the papers selected by them for presentation at the congress to the 
Executive Committee of IFAC. 


4) All papers submitted by the National Committees for presentation at the congress and approved by the 
Executive Committee of IFAC will be published in the Proceedings of the Congress in the English and Russian 
languages, 


5) Individual scientists and engineers of all countries may take part in the congress and present papers 
without submitting them to their National Committees, 


However, special permission of the Executive Committee of IFAC is necessary for all such papers. 
6) It is desirable that the size of typed page be not greater than half a printed page, including illustrations. 
1) Papers are to be given either in Russian or in English, 


4. Procedures for Preparing and Conducting the Congress 


1) The congress opens June 23, 1960. Duration of the congress will be ten days. 


?) The deadline for the submission of papers by the National Committee of the Soviet Union is October 31, 
1959, 


3) Papers prepared for printing and for authors’ corrections by December 31, 1959, 

4) Printed by March 31, 1960, 

5) Distribution of papers to the authors and participants in the congress on April 30, 1960, 
6) Over-all congress schedule. 


The first day will be devoted to the plenary session, which will include: a) the opening of the congress; 
b) a paper on the subject, ''The science of automatic control and the development of technology.” 


The following days, from the second to the ninth, inclusive, will be given over to the work of the Sections, 


On the final day there will be a plenary session devoted to reports of the Section leaders, current Federa-~ 
tion matters, the passing of resolutions, and the presenting of various addresses. 


1) The dividing up of the congress into Sections must be carried out after the papers have been obtained, 
and will be dictated by the scientific content of these papers. 
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Instructions for Congress Participants 


The National Committee of the USSR on Automatic Control is now accepting applications for the present- 
ing of papers at the First Congress of IFAC. 


All those desirous of presenting papers should submit their applications by February 15, 1959 to the 
National Committee: Moscow, Kalanchevskaia 15a, Institute for Automation and Remote Control, AN SSSR. 
The applications must contain the title of the paper and its abstract (without illustrations), taking no more than 


three typewritten pages. 
The complete text of the papers must be submitted to the National Committee no later than May 15,1959,* 


The papers must be prepared in accordance with the rules for preparation of papers for publication in the 
journal, “Automation and Remote Control" (USSR) as given in the issue of No, 1, 1959, If desired, additions to 
the papers, suitable for discussion purposes, may be prepared by the participants and sent to the National Com- 
mittee of the USSR no later than May 15, 1960, 


The National Committee of the USSR for Automatic Control requests all organizations and authors to pre- 
sent their papers, if at all possible, both in Russian and in English. 


All papers approved by the National Committee will be published in the Proceedings of the Congress; all 
papers approved by the Executive Committee of IFAC will be read at the congress as bases for discussions on the 
given questions, 


In order that the discussion have maximum scientific value, these papers will be printed and distributed 
to the delegates* * to the Congress within two months of closing. 


The Presidium of the National Committee 
of the USSR for Automatic Control 


* The deadline for submission of papers was dictated by the necessity of translating them into English, reviewing 
them, and submitting them to the Executive Committee of IFAC. 
* * The formation of the delegation from the USSR to the congress will be implemented subsequently. 
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